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Abstract It has been assumed untill very recently that all 
long-range correlations are screened in three-dimensional 
melts of linear homopolymers on distances beyond the cor- 
relation length ^ characterizing the decay of the density fluc- 
tuations. Summarizing simulation results obtained by means 
of a variant of the bond-fluctuation model with finite mono- 
mer excluded volume interactions and topology violating lo- 
cal and global Monte Carlo moves, we show that due to an 
interplay of the chain connectivity and the incompressibility 
constraint, both static and dynamical correlations arise on 
distances r ^ ^ . These correlations are scale-free and, sur- 
prisingly, do not depend explicitly on the compressibility of 
the solution. Both monodisperse and (essentially) Flory-dis- 
tributed equilibrium polymers are considered. 

Keywords Polymer melts • Monte Carlo simulations • 
Time-dependent properties 
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1 Introduction 

1.1 General context 

Solutions and melts of macromolecular polymer chains are 
disordered condensed-matter systems fP\ of great complex- 
ity and richness of both their physical and chemical prop- 
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Being of great industrial importance 



erties 

and playing a central role in biology and biophysics [8 9], 
they represent one relatively well-understood fundamental 
example of the vast realm of so-called "soft matter" sys- 
tems mOI comprising also, e.g., colloids III], liquid crys- 
tals [II and self-assembled surfactant systems [lOlQ The 
notion "polymer" is not limited to the hydrocarbon macro- 
molecules of organic polymer chemistry a la H. Staudinger 
or W.H. Carothers |I5 I6| but refers also to biopolymers 
such as DNA or corn starch I.I0. .8J and various self-assembled 
essentially linear chain-like supramolecular structures such 
as, e.g., actin filaments f9] or wormlike giant micelles formed 
by some surfactants 1 10, I7|. 

Obviously, dense polymeric systems are quite compli- 
cated, but since their large scale properties are dominated 
by the interactions of many polymers, each of these interac- 
tions should only have a small (both static and dynamical) 
effect. A sound theoretical starting point is thus to add up 
these small effects independently and to correct then self- 
consistently for the deviations due to correlations between 
the interactions [4 ||j Obviously, the number of chains a ref- 
erence chain interacts with, and thus the success of such a 
mean-field approach, depends an the spatial dimension d of 
the problem considered. Let A'^ be the number of monomeric 
repeat units per polymer chain (with A' ^ 100), 7?n ^ N^ 
their typical size, p^ w ^ l^n their self-density (also called 

' A soft matter system may be defined as a fluid in which large 
groups of the elementary molecules have been permanently or tran- 
siently connected together, e.g. by reversibly bridging oil droplets in 
water by telechelic polymers or similar systems of autoassociating 
polymers 112113114 ] . so that the permutation freedom of the liquid 
state is lost for the time window probed experimentally |10|. The 
thermal fluctuations which dominate the liquid state must thus coex- 
ist with constraints reminiscent of the solid state. Reflecting this solid 
state characteristics the dynamic shear-modulus /i(?) thus often re- 
mains finite up to macroscopic times t and the associated shear viscos- 
ity 7] ~ /Q°°df/x(f) may become huge. Most soft matter systems behave 
as liquids for very long times, i.e. /x(f) — > for f — > <» and 77 remains 
finite. 

^ As seen in Sec. |2] it is, e.g., necessary to "renormalize" |31 the 
local bond length / to the "effective bond length" b |4| or the second 
virial coefficient V2 of the monomers to the "effective bulk modulus" v 
fTSl to take into account the coupling of a reference chain to the bath. 




Fig. 1 Two-dimensional sketch of one chain of an idealized poly- 
mer melt represented by self-avoiding walks on a simple cubic lat- 
tice (d = 3). As notations we use r, for the position of a monomer ;, 
/j = £,_,_i — £, for its bond vector, / for the root-mean-squared bond vec- 
tor, £,„„ =£,„—£„ for the end-to-end vector of the subchain between 
the monomers n and m = n + s and r = \r„,„\ for its length. Taking 
advantage of the translational invariance along long chain contours, 
subchain properties, such as the 2p-th moment of £,„„, are generally 
sampled here by averaging over all possible pairs of monomers (n,m). 
The bond vectors between monomers are taken from a set of allowed 
lattice vectors. The compressibility of the melt may be imposed by 
means of a L agran ge multiplier e conjugated to the local density fluc- 
tuations (Sec. l2.3t . Two monomers of a chain at a distance r interact 
through an effective po tentia l v{r) as predicted by Edwards ||4j and 
further discussed in Sec. 12.41 



number Mmon = Vp, the volume V =L^ of the c?-dimensional 
simulation box and the temperature T are imposed. (Boltz- 
mann's constant ^b is set to unity throughout the paper.) A 
bond vector Z, = r,_^ j — r, connecting the monomers / and 
J + 1 is not necessarily restricted to a step of length I = a 
(with a being the lattice constant) to the next-nearest lattice 
site. As in the well-known bond-fluctuation model (BFM) 
123 "29 "3 01 presented in detail in Sec. [3] one may in gen- 
eral take bonds from a suitably chosen set of allowed lattice 
vectors to achieve a better representation of the continuum 
space symmetries of real polymers ITTlBlL We focus on sys- 
tems of high monomer number density p w a~'^ . In princi- 
ple, although not in practice as seen below, one could sample 
the configuration space of such a lattice polymer melt using 
a Monte Carlo (MC) algorithm 126112711 with local hopping 
move attempts to the 2d next-nearest neighbor sites on the 
lattice, as shown for the monomer k in panel (a). The specific 
local MC move illustrated is an example for the more gen- 
eral local and global rearrangements of the monomers 1271 
29 , 32 , 33 , 34 , 35 , 36 , 37 1 (respecting the set of allowed bond 
vectors) one may implement as discussed in Sec. 13. 31 Due to 
their simplicity and computational efficiency, such coarse- 
grained lattice models have become valuable tools for the 
numerical verification of theoretical concepts and predic- 
tions ll2ll2lf3QI. especially focusing on scaling properties in 
terms of the chain length A^ or the arc-length s = m — n < 
N -1 of subchains 13|. 



"overlap density") and V their inverse fractal dimension fT9|. 
Assuming the monomer number density p to be A^-indepen- 
dent, a mean-field theoretical approach may be hoped to be 
successful if 1 3 1 
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For to leading order Gaussian chains (V = 1/2) this implies 
that typically p/p^ ^ Ar«'/2-i = ^/^ ^ 1 chains interact in 
d = 3. Hence, dense three-dimensional (3D) polymeric liq- 
uids should be much easier to understand than normal liq- 
uids, say benzene and let alone water, where each molecule 
has only a few, say ten, directly interacting neighbors 12011 . 



1 .2 Coarse-grained lattice models for polymer melts 

Focusing on the generic statistical properties of flexible and 
neutral homopolymer melts in d = 3|j one of the simplest 
idealizations consists, as illustrated in Fig.[T] in replacing the 
intricate chemical chain structure by self-avoiding walks on 
a periodic lattice [,2,3,21 J. Let us assume canonical ensem- 
ble statistics ll22ll23lEtlE51l26l l27 1 where the total monomer 

^ Although we focus on 3D bulks the general rf-dependence is often 
indicated since this may help to clarify the intrinsic structure of the re- 
lations. The reader is invited to replace dby d = 3. For similar reasons 
we often make explicit the inverse fractal dimension v. It should be 
replaced by its value v = 1/2. 



1.3 Flory's ideality hypothesis for polymer melts 

One of the most central concepts of polymer physics is that 
in a polymer melt all long-range static and dynamical cor- 
relations are generally assumed to be negligible beyond the 
screening length ^ (defined properly in Sec. 12.3 1 and tested 
numerically in Sec. 14.4b characterizing the decay of the den- 
sity fluctuations |4|. The polymer chains are thus expected 
to adopt random- walk statistics as already stated by Flory in 
the 1940th Il2l [38]l39]l and worked out more systematically in 
terms of a perturbation theory by Edwards two decades later 
[i4i i40ii41ii42L431 . The (normalized) probability distribution 
G{r,s) of the end-to-end distance r of a subchain of a suffi- 
ciently large arc-length i should thus be a Gaussian, 
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= (. 
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for Is ^ r ^ ^ and N > s ^ g with b being the "effective 
bond length" and g the arc-length spanning the correlation 
length ^ . Note that in general the effective bond length b 
differs for various reasons from the root-mean-squared bond 
length I between monomers and ultimately must be fitted as 
will be shown in Sec. 15.31 One immediate consequence of 
Eq. (|2]l is of course that the second moment of G{r,s) must 
increase linearly, 

Rl = {{r„,=n+s-rnf) = b^s for ^ » g, (3) 

and the same holds for the mean-squared total chain end- 
to-end distance/?^ = ((r^ —LiY) with s =N — 1 KiN, i.e. 



V = 1/2 for the inverse fractal dimension. Since the bond- 
bond correlation function Pi{s) measures the curvature of 
the mean-squared subchain size R^, 



(l 






1 d^ 2 



(4) 



it follows from Eq. Q that Pi (s) = for ^ ^ gQ This sug- 
gests to test for possible deviations from R^ ^ shy comput- 
ing directly the bond-bond correlation function Pi (s) as we 
shall do in Sec. l5.4l That Eq. (|4| holds follows from the gen- 



eral relation for displacement correlations Eq. ( 11731) demon- 
strated in Appendix lArl 



1 .4 Rouse model hypothesis for polymer melts 

Based on Flory's ideality hypothesis for the chain conforma- 
tions, the celebrated reptation model suggested by Edwards 
and de Gennes mm provides a widely accepted theoretical 
description for the dynamical properties of entangled poly- 
mer melts explaining qualitatively the observed increase of 
the melt viscosity T] with respect to the chain length and 
other related observables 145 L The central idea is that chains 
cannot cross each other and must due to this topological 
constraint "reptate" along an effective "primative path" set 
by their own contourQ Using local topology conserving MC 
moves various authors have thus attempted in the past to ver- 
ify the reptation predictions by lattice MC models as the one 
sketched in Fig, ft] 1 50 !l5T][52l|53l|54l|5l]|56 57] . Similar tests 
have also been performed by means of molecular dynamics 
(MD) m, Brownian dynamics (BD) EJl and MC simu- 
lations of various off-lattice bead-spring models Il58ll59ll60l 

IMMllMllSllZI 

'^ This lack of memory of the bond orientations may be taken as 
the defining statement of Flory's ideality hypothesis. Using the central 
limit theorem 111] one immediately gets back to Eq. ^. 

^ It is also common to define the bond-bond (or angular) correlation 
function as the first Legendre polynomial P] (.s) = {e„ ■ e,„) using the 
normalized tangent vector e„ = ?„/|/„|. The difference between Pi {s) 
and Pi {s) is irrelevant from the scaling point of view and negligible 
quantitatively for all computational models we have studied. Higher or- 
der angular correlations may be described in terms of general Legendre 
polynomials P, 1 44 1 sampled as a function of the arc-length s between 
the bond pair or their spatial distance r. The polynomials Pair) may 
be used to characterize the interchain angular correlations, a question 
out of the scope of the present study. 

^ In its original formulation the reptation model is a single chain 
mean-field theory, the constraints due to other chains being accounted 
for by the effective topological tube |4|. Due to density fluctuations 
generated by the motion of the chains and the finite compressibility 
of the elastic mesh, different chains must be coupled, however. As 
suggested by A.N. Semenov et al. |46 47 1, the relaxational dynamics 
of long reptating chains (N ^ N^ in d = 3) is ultimately "activated", 
i.e. the longest relaxation time and the viscosity inc rease e xponentially 
with chain length N and not as simple power laws I48I49I . 

^ The computational activities have focused recently on the latter 
off-lattice models, one essential advantage being that the experimen- 
tally relevant shear-stress correlation function /x(f) can be computed 
directly without using additional input from theory. 



We remind that reptational motion along the primative 
path is assumed to be the dominant relaxational process only 
for chains much longer than the (postulated) entanglement 
length A'e and for times t larger than the entanglement time 
Te ^ N^ |3|. In the opposite limit the dynamics is thought 
to be described by the Rouse model proposed in the 1950th, 
i.e. by a simple (position) Langevin equation for a single 
chain with all interactions from the bath being dumped into 
a friction term and random uncorrelated forces |3|. Thus the 
reptation model "sits on top" of the Rouse model and should 
the latter model prove not to be sufficient this must be rele- 
vant for long chains in the entangled regimeQ 

Interestingly, it is relatively easy to construct for the lat- 
tice models sketched in Fig.[T]local monomer moves switch- 
ing off the topological constraints without changing the ex- 
cluded volume and other interactions and conserving thus 
all static properties. This has allowed to demonstrate that 
topological constraints are paramount for the relaxation pro- 
cess Il54ii55j . Since the local MC hopping attempts corre- 
spond to the random white forces of the Rouse model and 
since all other long-range correlations of the bath are sup- 
posed to be negligible beyond ^ , one generally assumes that 
such a topology non-conserving lattice model should be de- 
scribed by the Rouse model for chain lengths N ^ g. Since 
the forces acting on the chain center-of-mass (CM) tjq{t) 
are uncorrelated, the probability distribution of the CM dis- 
placements must be a Gaussian just as Eq. (|2]l with the CM 
displacement r^{t) — £n(0) replacing r, the time t replacing 
the arc-length .$ and 2dD^ replacing b^. The self-diffusion 
coefficient Dn and the longest Rouse relaxation time T^ of 
chains of length A^ should scale as ||4j 

Z)n « b'^W/N, Tn « 4/Z)n « N'^^'+'/W (5) 

with W being a convenient local monomer mobility defined 
properly in Sec. l6.3l l5ri. In analogy to Eq. (O, Dn may be 
obtained by fitting the mean-square displacement (MSD) of 
the chain CM 

%(f) = ((/:n(0 -z:n(0))'> «2JZ)Nf. (6) 

The chain relaxation time Tn corresponds to the motion of 
the CM over the typical chain size and may be defined by 
setting /%(? = Tn) = R^ \5\\. Note that Eq. © should hold 
for all times t (taken apart very short times corresponding to 
displacements of order of £, or the lattice constant a) and not 
only for / ^ Tn as for reptating chains. In analogy to Eq. (|4| 
one may characterize the CM displacements by a four-point 
correlation function in time, the "velocity correlation func- 
tion" (VCF) 

m(0)\ \d^h^{t) 



Cn(0 = 



uit) 



dt dt 



2 dt^ 



5rforf>5f, (7) 



As a consequence of the "hydrodynamic screening hypothesis" for 
dense polymer solutions |4|, all Zimm-like long-range correlations are 
generally regarded to be negligible beyond the static screening length t, 
166.67.6 81. The monomer momentum is not lost, of course, but trans- 
mitted to the bath and for large times and distances the melt is still 
described by the Navier-Stokes equation 1 1 1. Although these hydrody- 
namic effects are beyond the scope of this paper focusing on lattice MC 
simulations, a short hint on scale-free correctio ns to the hydrodynamic 
screening hypothesis is given at the end of Sec. 17.21 



which measures directly the curvature of h-f^{t) with respect 
to time with u{t) = r^{t + 5t) — r-^{t) being the CM dis- 
placement for a time window 8t. Since in a MC algorithm 
there is no velocity in the sense of an Euler-Lagrange equa- 
tion, "velocity" refers to the displacement per time incre- 
ment dt. As shown in Appendix lAl Eq. (|7]) states that Cn(0 
becomes independent of dt for t^ 5t. Since Eq. (|6]l implies 
that Cn(0 must vanish for all times, this begs for a direct nu- 
merical test as will be discussed in Sec.|6] Interestingly, the 
above statements do not only hold for the displacement field 
of the total chain CM r-^{t) but also for the displacements 
of the CM r\^(t) of arbitrary subchains of arc-length i for 
times f <C Ts w s^^^^JW with T^ being the Rouse relaxation 
time of the subchainO We can thus state more generally that 
according to the Rouse model one expects 



function is predicted to decay in (effectively) d dimensions 
as IMlSll^llTl 

1 /(0\'a Coo 1 



/>!(.) = 



2ft) 



ft;\' 

71 J 



bdp'. 



with (0 = dv 



(11) 



for g <^ s <^ N with p being the (i-dimensional monomer 

densityI3 

At variance to Eq. ^ the displacement correlation func- 
tions C-^{t) and Cs(f ) of chains and subchains are observed 
not to vanish. Instead an algebraic long-time power-law tail 



Q(0 « - 



b^W 1 



(m)^'" with ft; = 



2 + d 



= 5/4 (12) 



Q{t) =OfoTt<t:T,andg<^s<N-l 



(8) 



with Cs(f) being the VCF associated to the subchain CM 
displacements u(t) = r^(t + 5t) - r^it). 



1 .5 Aim of this study and key results 

Summarizing recent theoretical and numerical work made 
by the Strasbourg polymer theory group but also by other 
authors I3|69]|70ll2l]|721|73]|Ell2l]|7Sl|77l^^ we address 
and question in this review the validity of the general screen- 
ing assumption for dense polymer solutions and melts both 
concerning Flory's ideality hypothesis for the chain confor- 



b'ip^'"' """ 2+1/v 

is found 'md = 3 for alls <N andt <^Tsf90\. Due to Eq. 
this finding implies for .s = A'^ - 1 that the CM MSD of chains 
/jn(0 takes an additional algebraic contribution which dom- 
inates for short times such that 

%(?) - t^ with j3 = 2 - ft) = 3/4 for f < f* - Af° (13) 

in qualitative agreement with the anomalous power-law ex- 
ponent j3 ~ 0.8 obtained in various numerical Il51ll55ll56ll57] 
|59'i61 62 98 99 100 101] and experimental studies |77,78J 
|79i . Our theoretical and numerical result, Eq. ( 1121) . may thus 
help to clarify a long-lasting debated] 



1.6 Outline 



We begin this review by summarizing in Sec.|2]a few well- 



mations i0l[8l][82l|83]|84ll83[86ll87ll88]|8l and the Rouse mods^^own theoretical concepts from polymer theory and outline 



assumption for the dynamics of systems with irrelevant or 
switched off topology 190 "911 1921193] . We argue that due to 
the overall incompressibility of the polymer melt both the 
static and the dynamical fluctuations of chains and subchains 
are coupled on scales beyond the screening length ^ . In li = 
3 dimensions this leads to weak but measurable scale-free 
deviations which do not depend explicitly on the (low but 
finite) isothermal compressibility K"t of the solution. 

The static correlations are made manifest by the find- 
ing 1 80 84 85 86 1 that the intrachain bond-bond correlation 
does not vanish as implied by Flory's hypothesis, but rather 
decays analytically ind = 3 dimensions as 



^iW 



^WithCp: 



Co.c,/8 for g<^s<^N 



(9) 



the amplitude cp being given by the dimensionless stiffness 
coefficient Coo = (b/l)'^ |H| and the so-called "swelling coef- 
ficient" Cs = ^/TA/n^/pb^ [84|. Due to Eq. © this is con- 
sistent with a typical subchain size 

R^ =b^s {I ~cJ^/^) for g<t:s<^N, (10) 

i.e. R^/s approaches the asymptotic limit b^ monotonously 
from below and the chain conformations are thus weakly 
swollen as shown in Sec. 15.3.21 Note that generalizing Eq. (|9]l 
as demonstrated in Appendix IB.4l the bond-bond correlation 

' For times t > T^ a. subchain feels the connectivity with the rest of 
the chain and must follow the typical monomers motion for f ^ Tn- 



the perturbation calculations which lead, e.g., to Eq. (|9]l. 

Since the predicted deviations Eq. ^ and Eq. (fT2l) are 
rather small in d = 3 dimensions, the key challenge from 
the computational side is to obtain high precision data for 
well-equilibrated polymer melts containing sufficiently long 
chains to avoid chain end effects. This review focuses on 
numerical results obtained by means of the BFM algorithm 
with finite excluded volume interactions [86| and topology 
violating local an d g lobal Monte Carlo moves 1. 84 J as de- 
scribed in Sec. l3FT 0ur numerical data for monodisperse 
melts are crosschecked Il81ll83 i87"881 using systems of an- 
nealed size-distribution, so-called "equilibrium polymers" 
(EP) n7l ll02lll03lf!04yi05lll06lfT07]ll08lll09lfnol[TTTl . where 
chains break and recombine constantly and for this reason 
the sampling becomes much faster than for monodisperse 
quenched chains (Sec. l3.5l l. 

Static properties are discussed in Sec.|4]and Sec.|5j In the 
first section we investigate some thermodynamic properties 

'" For d <2 th e rang e of validity of Eq. dlU is further restricted as 
discussed in Sec. 17. 2.11 

' ' We emphasize that Eq. dllt and Eq. <12t are both scale-free and do 
not depend explicitly on the compressibility Kj of the solution. Inter- 
estingly, both key findings depend on the subchain length s considered 
and not on the total chain size A' as long as s ^N. 

'^ To check for caveats due to our lattice model we have verified our 
results by MD simulation of a standard bead-spring m odel with and 
without topological constraints [80.82.841851911921931 . In order not 
to overburden the present summary, this story must be told elsewhere. 



of dense BFM solutions such as the chemical potential /iN in 
EP systems in Sec. 14.61 1 88 1. We verify in Sec. l4.4l the range 
of validity of the (static) "Random Phase Approximation" 
(RPA) 14 1 describing the coupling of the total density fluctu- 
ations, encoded by the total structure factor S{q) at wavevec- 
tor q, with the degrees of freedom of tagged chains, encoded 
by the intrachain coherent single chain form factor F(^) (5 1. 
Section |5] contains the central part of our work related to 
conformational intrachain properties where we demonstrate 
Eq. (O. The deviations from Flory's ideality hypothesis are 
also of experimental relevance, as shown in Sec. 15.71 since 
the measured form factor F(^) should differ from the corre- 
sponding Gaussian chain form factor Fo{q) as II82II83I 



1 



1 



F{q) Foiq) 



g|^forl//?N«?«l/^ 



(14) 



Our work on dynamical correlations in Rouse-like sys- 
tems without topological constraints |90 91 1 is summarized 
in Sec.|6l A simple perturbation theory argument leading to 
Eq. (1121) will be given in Sec. 16.61 This argument uses the 
■'dynamical Random Phase Approximation" (dRPA) 



113 ','91l describing the coupling of the degrees of freedom of 
the bath (encoded by the dynamical structure factor S{q,t)) 
to the degrees of freedom of a tagged test chain (encoded by 
the dynamical form factor F{q,t)). We shall verify explicitly 
this important relation. 

We conclude the paper in Sec. |7] where we comment on 
related computational work focusing on dimensionality ef- 
fects (Sec. 17.2. lb |94 95l [TT4ll96]| . collective interchain cor- 
relations (Sec. 17.2.2b I18II115II and long-range viscoelastic 
hydrodynamic interactions (Sec. 17.2.3]) I.92..93 I. Several the- 
oretical issues are relegated to the Appendix. 



2 Some theoretical considerations 

2.1 Introduction 

Let us go back to the simple generic lattice model sketched 
in Fig.[I] Simplifying further, we begin in Sec. l2.2l bv switch- 
ing off all short- and long-range interactions between chains 
and monomers, the only remaining interaction being the con- 
nectivity of the monomers along the chain contours. Transla- 
tional invariance along these contours is assumed. No partic- 
ular meaning is attached to the orientation of the monomer 
index i]}^ To make this toy model more interesting let us on 
the other side introduce two additional features: local chain 
rigidity (Sec. l2.2.TI) and polydispersity (Sec. 12.2.21) . Remind- 
ing some standard properties |l3]|4l|5]|7][8l we will then intro- 
duce for later reference the intrachain single chain form fac- 
tor F{q) (Sec. l2.2."4] l and the total monomer structure factor 
5'(^) (Sec. 12.2.41) . Since the chains do not interact we can 

'^ This implies / -H- — / symmetry with respect to the monomer in- 
dex i which can be read as a time variable t as seen in Appendix IaI 
Due to this reversibility the bond-bond correlation function can be ex- 
pressed in terms of the n on-Gaussian "colored forces" [116] acting on 
the monomers, Eq. (176). 



consider each chain independently. The interaction between 
chains and monomers will be switched on again in Sec. 12.31 
by means of a Lagrange multiplier £ limiting the density 
fluctuations (Sec. l2.3.2] l. Note that assuming Flory's ideality 
hypothesis most intrachain properties discussed in Sec. 12.21 
should also hold rigorously in incompressible melts with the 
effective bond length b being the only fit parameter. The ef- 
fective entropic correlation hole forces for chains and sub- 
chains arising due to the incompressibility constraint are in- 
troduced in Sec. 12.3.41 The generic scaling of the deviations 
from Flory's ideality hypothesis is motivated in Sec. 12.3.61 
before we turn to the systematic perturbation calculation in 
Sec. 12.41 Following Muthukumar and Edwards 114311 we ar- 
gue that the reference length bo of the Gaussian reference 
chain of the calculation (Sec. IB .21) should be renormalized 
(Sec. 12.4.51) . The bond-bond correlation function Pi {s) for 
asymptotically long chains is presented in Sec. 12.4.61 before 
we turn finally in Sec. l2.4.'7l to finite chain size effects. 



2.2 Connectivity constraint 

2.2.1 Local rigidity 

Keeping only the chain connectivity we switch off all other 
monomer and chain interactions. Let us apply a local stiff- 
ness energy proportional to the cosine of the bond angle 6, 
cos(0) = e„ ■ e„+i, with e„ — Ln/\Li\ being the unit tangent 
vector. The stiffness energy is assumed to be not too large to 
avoid lattice artifacts f52l. Due to the multiplicative loss of 
any information transferred recursively along the chain con- 
tour the bond-bond correlation function Pi{s) must decay 
exponentially with arc-length i IJ), 



Pi{s) Kiexp{-s/sp) 



(15) 



with Sp = Ip/l being the curvilinear persistence length, i.e. 
P\{s) w for i ^ Sp. Using Eq. (|4]l it follows of course 
that /?s « b^s for subchains of arc-length s ^ Sp. Assum- 
ing Eq. dTSl) this implies b^ = fi{2sp — 1) for the effective 
bond length b since more generally it is known that 1291 



= /2f2f;PiW-lj=2ZZp-i 



(16) 



The latter relation is consistent with the more common defi- 
nition |4 



b^ = 2da^ 



;lim^ 



(17) 



using the total chain mean-squared end-to-end distance iU. 
We have introduced here a convenient monomeric length a 
allowing to simplify prefactors depending in a spurious man- 
ner on the spatial dimension d and have reminded the dimen- 
sionless stiffness parameter Coo = {b/l)^. 



'"* Definitions based on such "Einstein relations" 
more robust numerically. 



1171 are generally 



We note that a good approximation for a "freely rotating 
chain" with local stiffness potential is given by ID 

l + (cos(0)) 

1 — (cos(0)) 

which has been shown to be useful even for lattice models 
with discrete bond angles ll52l . In summary, a local stiffness 
energy does not change the scaling of the chain and sub- 
chain size with, respectively, A^ or s, but merely increases 
its amplitude. Note that a weak local chain stiffness with 
ip w 1 would be generated by disallowing chains to return 
after, say, two steps to the same lattice site (so-called "non- 
reversal random walks" [25 1) and similar local constraints. 
Since these are the only local stiffness contributions which 
may occur in the presented numerical work these small rigid- 
ity effects can be safely disregarded belowQ 

2.2.2 Polydispersity 

Let us relax the monodispersity constraint made above and 
assume a general (normalized) probability distribution pN 
for chains of length A^. We suppose that the mean chain 
length (A'^) is arbitrarily large and that all moments (A^'') = 
/(5° dN A"^ pn of the distribution exist. For subchain proper- 
ties everything remains as before. Let us call R^ the typical 
mean-squared chain size of chains of length A'^. Since one 
averages experimentally over chains of different length this 
total average depends now on the moment p of the distribu- 
tion which is probed. 



R 



\b^N, 



P' 



(19) 



(NP) {NP] 

with Np being the p-averaged chain length. Since for exper- 
imentally reasonable distributions (A^'^) ^ (A'^)'^ ^ Nj^, we 
have normally R^ ^ b^ (N) with numerical coefficients due 
to the moment taken. For the properties considered by us, it 
is the so-called "z-average" for p = 2 which matters most. 
The z-average is for instance probed by sedimentation in an 
analytical ultracentrifuge or in neutron scattering measure- 
ments of the intramolecular form factor F{q) as further dis- 
cussed in Sec. 12.2.41 1151. Note that all standard definitions 
and formulas 001 are recovered for a monodisperse length 
distribution pN = 5{N' -N). 

For the important case of Flory-distributed polymers with 

p^^pi exp(-juAf) with {N) = jx~^ (20) 

we have {NP) = p\ {N)P and thus R^ = {p+l)b^ (N). Such 
a Flory distribution is expected for systems of EP with an 
annealed length distribution where a constant finite scission 
energy E >0 has to be paid for the scission of each bond 
as described in Sec. 13. 51 This can be seen by minimizing the 
Flory-Huggins free energy functional ll9 HT7lll07]| 

F[pN]=V'£pN(Tlog(pN)+liN + E + 5iit^) (21) 

N 

'^ The only case where t he sma ll local chain rigidity qualitatively 
matters is presented in Sec. 15.4. 2l for the bond-bond correlation func- 
tion Pi (r) for large distances r 2> r, between bonds (87). 



with respect to the density pN ^ PPn/ (N) of chains of length 
A'^. The first term on the right is the usual translational en- 
tropy. The second term entails a Lagrange multiplier which 
fixes the total monomer density 



IA^Pa 



(22) 



N=l 



All contributions to the chemical potential of the chain /i^ 
which are linear in A'^ can be absorbed within this Lagrange 
multiplier. The scission energy E characterizes the imposed 
enthalpic free energy cost for breaking a chain bond. The last 
term 5/Jn encodes the remaining non-linear contribution to 
the chemical potential jUn which has to be paid for creating 
two new chain ends P^ For non-interacting random walks on 
the lattice this contribution is just a (model depending) con- 
stant entropic factor which renormalizes the imposed scis- 
sion energy to an A'-independent effective scission free en- 
ergy E + 5/Jn II109I . The minimization of Eq. (1211) under 
the density constraint, Eq. (|22] |. yields Eq. (|20] | with a mean 
chain length 

M-> = (N) = {{N") /plf" « ^pcxp{E/T) (23) 

as the reader will readily verify. 

2.2.3 Segmental size-distribution G{r,s) 

Due to the translational invariance in space and along the 
chain contours most perturbation calculations outlined be- 
low Il83ll84l i87 88 1 are more readily performed in recipro- 
cal space. The Fourier transform of a function /(r) is de- 
noted f{q) = ,^[f{r)] = ^ drj{r}e-''i'- and we write /(f) = 

■^[/('*)] = /o°/('*)^~" fo'" the Laplace transform of a func- 
tion f{s) with f being the Laplace variable conjugated to 
the arc-length s. We have introduced in Sec. 11.31 the prob- 
ability distribution G{r,s) of the end-to-end distance r of a 
subchain of arc-length i between the monomers n and m = 
n + s of a chain. The Fourier transform of this two-point in- 
tramolecular correlation function is thus in general G{q, s) = 
(exp(— i^ • r)) with the average being taken over all possible 
subchain vectors L = Lin'~Ln- For (infinite) Gaussian chains 
this becomes p^ 

G{q,s) = Go{q,s) = cxp{~(aqfs). (24) 

The index has been introduced for the general case where 
G{q,s) may differ from the Gaussian propagator Go{q,s) 
used in our perturbation calculations. Moments of the distri- 
bution G{r, s) are readily obtained from derivatives of G(^, s) 
taken at ^ = as recalled in Appendix IB . 1 1 II 1 1 61 . It follows 
for instance for the 2p-th moment of the distribution that 



r2/A 



= (-l)M^Go(,,.)L.o = ^^^^^''«^"(25) 



6P 



for Gaussian chains. We have thus 



Kpis) ^ 1 - 



pl6P {r}P) 
{2p+\)\ {bh)P 



= 0. 



(26) 



'^ As shown in Sec. 14.61 and Sec. IB. 61 the free energy co ntribution 
5/Xn may depend in general on the chain length N I6I7I106I . 



A very closely related characterization of G{r,s) is given by 
the standard non-Gaussianity parameter 



apis) = I 



pl6P 



^2p\ 



{2p+l)l{r 



2\P 



(27) 



which for Gaussian chains is identical to Kp{s). As further 
discussed in Sec. 15.51 OCp{s) has computationally the advan- 
tage that the effective bond length b must not be known 
a priori. Obviously, the general distribution G{r,s) is fully 

determined by either the dimensionless moments Kp{s) or 
, 



^p{s) ma. 

In Fourier-Laplace space the Gaussian propagator reads 



Go{q,t) 



1 



{aqy 



(28) 



If one needs to average over all bond pairs at a given dis- 
tance r irrespective of their curvilinear distance s, and one 
has thus to sum over all possible .s as in Sec. IB. 51 this cor- 
responds (for arbitrarily large chains) to setting f = for the 
corresponding Laplace variable. The summed up Gaussian 
propagator for infinite chains is thus Go{q,t = 0) = {aq)~^. 
For Flory-distributed Gaussian chains of finite mean chain 
length {N) = l/fJ. we have more generally a summed up 
Gaussian propagator 



G{q) = Goiq,t = 0) 



1 



(fl^)2+jU' 



(29) 



Inverse Fourier transformation yields ind = 3 the density 



G{r) = GQ{r,t = 0) 



1 



4na^r 



-^r/a 



(30) 



around a tagged reference monomer of all the monomers be- 
longing to the same chain. Obviously, for fi ^ °° one recov- 
ers the well-known density 



1 



„d-l/v 



ford > 1/v 



(31) 




1/R 



N 



thermal blob , , 

l/^(e) 1/a 



Fig. 2 Schematic representation of total structure factor S{q) and intra- 
chain coherent form factor F{q) for monodisperse chains (thin line) in 
double-logarithmic coordinates. For a large Lagrange multiplier £ ^^ oo 
we have S{q) ~ g(e = oo) ^ 1 for g ^ I/ct (dashed bold line). For fi- 
nite e the incompressibility is only felt for q<^l/^ i^)^^}^ ^^^ screen- 
ing length i^(e) setting the size of the "thermal blob" BI4I . 



2.2.4 Intramolecular coherent form factor F [q) 

The two-point correlation function G{q,s) may be probed 
experimentally through the intramolecular form factor F{q) 
which for monodisperse chains is defined by |14| 



1 ^ 
^^^^=N L (exp(-i^-(z:„ -£„))> 



1 



H,m— 1 



= -{C' + S') 



(32) 
(33) 



with C -- 



I^^^i cos(^ • r„) and S = Yin^i sm(q ■ r,,). The sec- 



n.m=i 



for infinite objects of inverse fractal dimension V |3 10|1_| 
This power-law dependence of the local density means that 
our random-walk polymer is a fractal set in the sense of 
Mandelbrot [19], i.e. the average mass (number of mono- 
mers) within a distance r of an arbitrary point of the set yv 

varies as a power of r. The r-dependence of the local den- R^fN) = V ((r — r ) 

sity reflects a type of spatial order that is not connected to 
translation or rotation symmetries but to the "dilation trans- 
formation" to a magnified system [3]. A structureless, uni- 
form material looks the same when magnified, provided the 
magnification is too weak to see the molecular constituents 
("lower cutoff')- The scaling exponent d—l/v characterizes 
the dilation symmetry in the same way that linear momen- 
tum characterizes translational symmetry and angular mo- 
mentum characterizes rotational symmetry 1101 . 



ond representation being an operation linear in A' has ob- 
vious computational advantages for large chain lengths. In 
principle, the characteristic polymer size may be obtained 
experimentally from F{q) in the Guinier regime for small q 
151. Expanding Eq. ^ yields F{q) =N{1- {Rg{N)qf/d) 
with 7?g (A') being the gyration radius defined as 



(34) 



For chains following Gaussian statistics (as our lattice chains 
with switched off interactions) where R^ = ((r„ —Lmf') ~ 
b^s, it follows by integration of Eq. (ODl that 7?^ (A^) ^R^/6^ 

da^N/3. For large wavevectors where the internal fractal 
chain structure is probed, the form factor decays as 



'^ The reader may verify the indicated scaling by direct integration 
of Eq. (|2} with respect to s for monodisperse chains with A' ^- oo. 



F{q) = 



[aqf 



' N^ for qR„ > 1 



(35) 



as sketched by the thin line in Fig. I2l I More generally, the 
form factor of monodisperse Gaussian chains is given by 
F{q) —NfD{x) withx = {qRg{N))^ and the Debye function 



/d W = -J (exp(-x) - 1 +x) . 



(36) 



For convenience of calculation, the Debye function is often 
replaced by the Pade approximation 



Fiq): 



N 



\ + {qR^)V2 {aqY + 2/N 



(37) 



where the last step holds for c/ = 3. 

For a general mass distribution pN the form factor is de- 
fined as asun 

P(^) = 7TJxT.P^ L (exp(-i£ •(£„-£,„))). (38) 

\"/ N=0 n.m=l 



In practice, Eq. 



corresponds to the computation of the 



averaged sum over contributions (CP' + 8^^ for each chain 
which one divides finally by the total number n^oa of labeled 
monomers 1.83 J . In the small-^ regime one obtains again a 
Guinier relation 



F{q): 



in 

(N) 



{Rg.z qf 



for /?g z ^ ^ 1 



(39) 



where /?g z stands now for the z-averaged gyration radius f6j 
where the moment p = 2 is taken over the standard radius 
of gyration R^iN), Eq. ([34]i. Note that for Flory-distributed 
Gaussian chains we have Rgz = da^/n and the form factor 
becomes [83 1 

Fiq) = 



[aqY+li 



(40) 



Eq. ( 140b reduces to Eq. ( I35l l for large wavevectors, as one 
expects, since in this limit F(^) must become independent of 
the length distribution /?n. Note that Eq. ( 140b has the same 
form as the Pade approximation for monodisperse chains, 
Eq. (O. 

2.2.5 Total monomer structure factor S(q) 

For non-interacting chains the local monomer density can 
freely fluctuate being restricted only by the chain connec- 
tivity. The isothermal compressibility Kj [1\ is thus given 
by the osmotic contribution due to the density of chains, i.e. 
TKj = \/{p/ {N)) diverges linearly with the typical chain 

'* The power-law scaling is obtained directly by Fourier transforma- 
tion of Eq. on which yields F{q) ~ l/q^l^ . This scaling only holds 
if the fractal object is "open", i.e. d > 1/v, and the scattering intensity 
is not dominated by the Porod scattering at the (possibly fractal) "sur- 
face" of a compact object 1 5 1 1 8 1 . The latter Porod scattering becomes 
relev ant, e.g., for self-avoiding polymers in strictly d = 2 dimensions 
I119U 20I and for non-olympic and un-kno tted rings in i/ = 3 or thin 
films of finite width (d = 2+) I 121I122I118 I. In both cases the poly- 
mers adopt compact configurations, 1/v = d, with a fractal surface of 
well-defined fractal surface exponent ds < d. The form factor thus de- 
cays as F{q)/N si I / {Rf^qf'' -''^ fTTSl . 



length (A') . More generally, one may characterize the fluc- 
tuations of the total density p{q) = L"=T^^P(~'^ '-«) i'^ 
Fourier space by means of the total structure factor 141 

1 '%ion 

S{q) = z L (exp(-i£- (£,,-£„,))) 



'mon „,,„=! 



I 



Y^cos{q-r„) + £sin(£ 



(41) 



with q being a wavevector commensurate with the simula- 
tion bojo and the thermal average (...) being performed 
over all configurations of the ensemble and all possible wave- 
vectors of length q = \q\. Since we have switched off all 
monomer interactions, monomers on different chains are un- 
correlated and, hence, 

^ = TT^ (42) 

S{q) F{q) 

as indicated for monodisperse chains by the thin line in Fig.|2l 
Due to the chain connectivity the fluctuations of chains (mea- 
sured by jct) and the fluctuations of monomers (measured by 
S{q)) may differ for polydisperse non-interacting systems: 
While S{q) ^ (N^) / (N) in the small-^ limit as seen from 
Eq. (|39] |, we have T Kjp = (N) for the compressibilitvPl 

2.3 Incompressibility constraint 

2.3.1 Dimensionless compressibility g 

Obviously, dense polymer solutions and melts are essentially 
incompressible, i.e. jct — > 0, and the above assumption that 
the total density can freely fluctuate is not very realistic. It 
is useful to introduce here a central dimensionless thermo- 
dynamic property characterizing the degree of density fluc- 
tuations on large scales, the so-called "dimensionless com- 
pressibility" 1831184118611 

g= — = lim {Tkt:p)= lim f lim5(^) ) . (43) 

At standard experimental polymer melt conditions g remains 
of course finite, say g w O.I, but typically well below unity 
as indicated by the bold dashed line in Fig.|2][8|. Note that 
we have defined g in the limit of asymptotically long chains 
to take off the trivial compressibility contribution due to the 
translational invariance of the chains mentioned above. For 
later reference we have also introduced in Eq. ( 143b the effec- 
tive "bulk modulus" v [ 1 15 1. As indicated, the dimensionless 
compressibility can be determined directly in experiments 
or in a computer simulation from the low-^ limit of the total 
monomer structure factor S(c{). This point is further elabo- 
rated in Sec. 14.41 



" If we use a cubic simulation box of linear dimension L the smallest 
possible wavevector is 2;i:/L. 

^^ The dimensionless compressibility g defined in Eq. l|43) for 
asymptotically long chains does not depend on this ideal chain contri- 
bution. Care is needed, however, if g is determined numerically from 
Siq) for polydisperse systems of finite (A'^) / (A') by e xtrapolation in 
analogy to the monodisperse case discussed in Sec. 14.41 



2.3.2 Lagrange multiplier £ 

Physically, the incompressibility of dense polymer systems 
arises of course due to the short-range repulsion of the mono- 
mers, i.e. it depends on non-universal physical and chemi- 
cal properties. From the theoretical and computational point 
of view it is, however, inessential how the incompressibil- 
ity at low wavevectors is imposedo This constraint could 
be achieved, at least in principle, by "simple sampling" [261 
of only those configurations respecting the chosen g. In this 
sense it is thus the "throwing away" of configurations from 
the extended configuration ensemble containing all possible 
linear chain paths on the lattice which creates the repulsive 
forces between chains, subchains and monomers Alterna- 
tively, one may design intricate local and global MC moves 
forcing the system through configuration space along a hy- 
perplane of constant g 1.27 ..291 . 

A more general and computationally more natural route 
is to use an extended ensemble ^7E\ and to impose the in- 
compressibility constraint through an external field with a 
Lagrange multiplier £ conjugated to the local monomer den- 
sity fluctuations. As may be seen in more detail in Sec. 13.41 
for a specific lattice model, this implies in practice that one 
has to pay an energy of order £ for the overlap of two mono- 
mers. While in the low-£ limit with g{e) ^ {N) the chains 
do not interact, i.e. S{q) =F{q) for all q, the incompressibil- 
ity constraint S{q) Ki g{£=°°) ^\ holds for all wavevectors 
q up to the monomeric scale {q w \/a) in the opposite limit 
£ — > oo. This is shown by the bold dashed line in Fig.|2] 



2.3.3 Thermal blobs of size ^ 

The situation is slightly more complicated for intermediate 
overlap penalties £ with 1 <C g[z) <^ (N) indicated by the 
thin dashed line. Since ^(£) is now a well-defined character- 
istic chain length in curvilinear space along the chain con- 
tour, it corresponds to a characteristic scale in real space, the 
"screening length" E, (e) of the density fluctuations defined 
as El 



1 2 a^ 
—a e 

2 ^ 2vp 



k2 ^ 2 " 



(44) 



where we use the effective bulk modulus v — \/{gp) follow- 
ing Eq. (143 br^ Generalizing Eq. (142b to systems with finite 

^' If the goal is to map a computational model onto a real polymer 
melt aiming to understand macroscopic properties, the starting point 
should be, in our opinion, to match the mechanical and thermodynamic 
properties in the low-g limit, e.g. the dimensionless compressibility g, 
rather than to fiddle with S{q) on the monomer scale. 

^^ The scale-free correlations described in this paper are thus akin to 
the effective "anti-Casimir forces" which arise in dense polymer melts 
due to the throwing away of configurations containing closed loops 
from an ex tended configuration ensemble with both linear chains and 
rings fTsTllSI . 

^^ From the scaling point of view a curvilinear length i » 1 translates 
quite generally to a spatial distance /- 2> a and a wavevector q<S^\/a 
according to r ~ !/<? ~ s^ with v being the inverse fractal dimension. 



compressibility the total structure factor is predicted to fol- 
low the so-called (static) "Random Phase Approximation" 

(RPA) aa], 



1 



S{q) 



1 

- + 
g 



1 



F{q) 



(45) 



(46) 



using Eq. (l35T l in the second step. The Ornstein-Zernike cor- 
relation equation 111] Eq. ( 146b justifies the above definition 
of the "screening length" ^ , i.e. density fluctuations decay 
in c/ = 3 exponentially as |41 

(p(r)p(0))-p2 = -^exp(-r/^). (47) 

We remind that ^ sets the size of the "thermal blob" IS cor- 
responding to a free energy due to the effective monomer 
interaction of order k^T .If one considers short subchains of 
arc-length s ^ g ox small distances r ^ ^ , the (sub)chains 
behave as if they were barely interacting, i.e. S{q) Ki F{q). If 
on the other side one focuses on the physical properties be- 
yond the thermal blob scale {s^ g, r^ ^,q <^ 1 /^ ) where 
the structure factor becomes constant, S{q) w g{£), one may 
renormalize all spatial distances by ^ and all curvilinear dis- 
tances by g and in these terms the system should behave as 
an incompressible packing of thermal blobs |3J. 

We emphasize that the perturbation results ||4| Eq. ( 1451) 
and Eq. ( 146b are supposed to apply only as long as the com- 
pressibility g is not too small and the screening length ^ 
remains a respectable length, in any case much larger than 
the lattice constant a^j Formally, this is expressed by the 
criterion [4,83,83 



p(^gl/2) 



pb^i 



«1, 



(48) 



with the Ginzburg parameter G^ being the small parameter 
of the standard perturbation theory. Note that Eq. ( 148b sets a 
lower bound to the correlation length ^ ^ ^/(pZ?^ )'''''" ^). 

Please also note that Eq. (|44]l, Eq. (|46]l and Eq. ( 1481) 
are consistent with relations given by Edwards |4|. (For in- 
stance, Eq. ( |48] | corresponds to Eq. (5.46) of Ref. |4|.) The 
only difference is that following 1 1 8 , 84 1 we have replaced 
the second virial coefficient V2 of the monomers by the effec- 
tive bulk modulus v and the bond length of the unperturbed 
chain by the effective bond length b. 

2.3.4 Correlation hole effects 

Let us for the clarity of the presentation return to monodis- 
perse chains in the large-£ limit, i.e. let us assume that g<^l 
and that thus the total monomer density p does not fluctu- 
ate. On the other hand, composition fluctuations of labeled 
chains or subchains may certainly occur, however, subject 

^'^ We remind that from the thermodynamic point of view the funda- 
mental property characterizing the solution is the dimensionless com- 
pressibility g and not the correlation length t, . 
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Repulsion 



Fig. 3 Effective interactions due to the incompressibility constraint: 
(a) Since the total density p can barely fluctuate, chains of length A' are 
known to repel each other due to an entropic penalty u^ » Pn/P ^^' 
by the self-density p^ ^ ^/^n- (b) Self-similar pattern of nested seg- 
mental coiTelation holes of decreasing strength !(* w s/pRf ^ Cj/y^in 
d = 3 dimensions aligned along the backbone of a reference chain. The 
large dashed circle represents the classical correlation hole of the total 
chain (s rs N). This is the input of some re cent approaches to model 
polymer chains as soft spheres I123I124I9 9I. We argue instead that 
the incompressibility constraint on all length scales s matters — and 
not just ioT s ^ N — leading to a short distance repulsion m* of their 
"subchain correlation holes" which increases with decreasing s. 



to the total density constraint. Composition fluctuations are 
therefore coupled and chains and subchains must feel an en- 
tropic penaky when the distance r between their CM be- 
comes comparable to their typical size 194 1821 . 

As sketched in Fig. Oa), let us first remind the well- 
known "correlation hole" effect for two test chains of length 
N in the bath |4,94|. The scaling of their effective interac- 
tion under the incompressibility constraint is obtained from 
the potential of mean force u{r,N) = — ln{p{r,N) / p{oo^N)) 
II117II with p{r,N) being the pair correlation function of the 
chains, i.e. the probability distribution to find the CM of the 
second chain at a distance r assuming the CM of the first 
chain at the origin (r = 0). Since the correlation hole is shal- 
low for large A', expansion of the logarithm leads to 



uir,N) w 1 



p{r,N) 
p{°o,N) 



1 



p-pN(r) PN(r) 



(49) 



with pt^{r) being the density distribution of the reference 
chain around its CM. This distribution scales as 



pNM«Pi^/(r//^^ 



(50) 



with p^ w N/R^ being the chain self-density, Eq. O, and 
/(x) a universal function which becomes constant for x <C 1 
and decays rapidly for jc ^ 1 13 The interaction penalty for 
two chains at r/7?N ^ 1 is thus given by 



which decreases as u 



:7V/p<-7V 



A-d/l 



(51) 



N -^ \/y/N md = ?> 19411821. Hence, 
although the incompressibility constraint couples the chains, 
the correlations become rapidly negligible with increasing 
chain length A'^ in agreement with our discussion in Sec. II. 11 
Interestingly, the above scaling argument does not only 
hold for chains {s = N — \) but also for the potential of 

^' For (to leading order) Gaussian chains fix) must also be Gaussian 
t3l as indicated by the bold line in Fig.|4] 
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Fig. 4 Potential of mean force u{r.s) for 3D polymer melts obtained 
using the classic al BF M algorithm without monomer overlap (e = oo) 
described in Sec. l3.3l for monodisperse chains of length N = 4096 and 
volume fraction (^ = 0.5 of occupied lattice sites 1 8 1]. Data for different 
s <N are successfully scaled tracing u{r,s)/ul as a function of r/R^ 
with i(* = s/pRJ. The potential is positive and roughly Gaussian (bold 
line) for small r/R^ where the (sub)chains repel each other It becomes 
weakly attractive for distances r/i?s ~ 0.7 (vertical arrow). 



mean force u{r,s) = —ln{p{r,s)/p{oo^s)) obtained in a sim- 
ilar manner from the pair correlation function p{r,s) of the 
center-of masses of subchains of arc-length s < N — I 
I84|. Since 



= [N/sYl'^-'^ = ^/nJs > 1 for i < TV, 



(52) 



the correlation hole effect strongly increases with decreasing 
s, albeit it remains always perturbative in c/ = 3. Note that 
the subchain correlation hole potential m* does not depend 
explicitly on the bulk compression modulus v. It is dimen- 
sionless and independent of the definition of the monomer 
unit, i.e. it does not change if A monomers are regrouped 
to form an effective monomer (p — !■ p/A, i -^ s/k) while 
keeping R^ fixed. 

That the effective correlation hole potential for chains 
and subchains is more than a heuristic scaling argument can 
be seen from Fig. |4l We present here the rescaled poten- 
tial of mean force a(r, s) for chains of length A' = 4096 ob- 
tained using the BFM algorithm described in Sec. 13.31 The 
correlation hole potential for the total chains is indicated by 
the squares. The predicted scaling is confirmed by the per- 
fect data collapse of u{r,s)/ul plotted as function of r/R^. 
Please note that, strictly speaking, u{r,s) for subchains ac- 
counts also for the attractive interaction between subchains 
on the same chain and therefore differs slightly from the ef- 
fective interaction potential of two independent chains of 
length N = s -\- \. This leads to the (very weak) attractive 
contribution to u{r,s) for subchains (barely) visible in Fig.|4] 
However, this additional effect does not affect the scaling on 
short distances, r <^ R^, which matters here. 
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2.3.5 Connectivity and swelling 

To connect two test chains of length A'^ to form a chain of 
length 2N an effective free energy m^ has to be paid and 
this repulsion will push the two halves apart from each other 
II94I . We consider next a subchain of length 5 in the middle of 
a very long chain. All interactions between the test subchain 
and the rest of the chain are first switched off but we keep 
all other interactions, especially within the subchain and be- 
tween the subchain monomers and monomers of surround- 
ing chains. The typical size 7?s of the test subchain remains 
essentially unchanged from the size of an independent chain 
of the same strand length. If we now switch on the interac- 
tions between the tagged subchain and monomers on adja- 
cent subchains of same length s, this corresponds to an effec- 
tive interaction of order a* as before. (The effect of switch- 
ing on the interaction to all other monomers of the chain 
is inessential at scaling level, since these other monomers 
are more distant.) Since this repels the respective subchains 
from each other, the corresponding subchain is swollen com- 
pared to a Gaussian chain of non-interacting subchains. It is 
this effect we want to characterize. 

2.3.6 Perturbation approach in three dimensions 

Let us return to systems of finite dimensionless compress- 
ibility g but let us focus on subchains of length 5 which are 
larger than the number of monomers g contained in the ther- 
mal blob, i.e. we focus on scales where the incompressibil- 
ity constraint matters. Interestingly, when taken ?As = g the 
subchain correlation hole potential becomes 

(53) 



size, Rl/b^s w 1 



must approach the asymptotic limit 



with Gz being the standard Ginzburg parameter already de- 
fined in Eq. ( 1481) . Hence, it follows for the subchain correla- 
tion hole potential that 

ul w G^ig/sfl'^-^ < Gz for t/ > 2 and 5 > g. (54) 

Although for real polymer melts as for computational sys- 
tems large values of Gz w 1 may sometimes be found, m* ^ 
\/\/s decreases rapidly with .s in three dimensions, as il- 
lustrated in Fig.[3fb), and standard perturbation calculations 
can be successfully performed. 

As sketched in Sec. 12.41 these calculations consider di- 
mensionless quantities '^[u*\ which are defined such that 
they vanish {"^[u* = 0] = 0) if the perturbation potential m* 
is switched off and are then shown to scale, to leading order, 
linearly with u* . For instance, for the quantity Kp{s), defined 
in Eq. ( 126b . characterizing the deviation of the subchain size 
from Flory's hypothesis one thus expects the scaling 



KpK] 



pRi' 



(55) 



The +-sign indicated marks the fact that the prefactor has to 
be positive to be consistent with the expected swelling of the 
chains. Consequently, the rescaled mean-squared subchain 



for large 5 from below. For 3D melts Eq. ( 155b implies that 
Kp{s) should vanish rapidly as l/(p^''v^)Ej Taking apart 
the prefactors — which require a full calculation — this 
corresponds exactly to Eq. dlOl) with a swelling coefficient 
Cs w \/pb^ . Note also that the predicted deviations are in- 
versely proportional to b^, i.e. the more flexible the chains, 
the more pronounced the effect. Similar relations "^[a*] ^ m* 
may also be formulated for other quantities and will be tested 
numerically in Sec.|5] 



2.4 Perturbation calculation 

2.4.1 General approach 

We remind that the first-order perturbation calculation of an 
observable ^ under a dimensionless perturbation potential 
C/ <C 1 (defined in units of T) generally reads 14] 



{^{^-U))o 

(l-C/)n 



K)o + ([/)o(^)o-(C/^)o (56) 



where averages performed over an unperturbed reference sys- 
tem of Gaussian chains are denoted (. . .)g. Obviously, (^) = 
(^)q for C/ = 0. For the reduced quantity ^ = si — {s^)q 
Eq. (156b simplifies to 



-{um) 



(57) 



For {s^)q 7^ one may introduce the dimensionless observ- 
able '^ = —SSj {s^)q which, of course, also obeys Eq. (|57] |. 
If {'^U)q w (f^)o' ^s is the case if self-energies and ultra- 
violet divergencies can be dumped into the reference of U , 
this is consistent with Eq. ( 155b . The observable s/ stands, 
for instance, for the 2/7-th moment s/ = rjm of the vector 
hm — Lm~Ln between two monomers n and m — n^s on 
the tagged chain as shown in Fig.[T][80 84| or for the scalar 
product si = /„ 'luJl^ of two bond vectors [87,89j. If, as in 
the latter case, we have {s/)q = for linear chains by con- 
struction, one only has to compute {s/) = i^SS) w — {U.'^)q. 
In this subsection we use b^ = yldao for the bond length 
of the unperturbed Gaussian reference chain which may a 
priori be different from the effective bond length b = y/lda 
defined and measured according to Eq. dlVb . The reference 
bond length bo is a parameter which may be suitably cho- 
sen or adjusted in a Hartree-Fock iteration scheme f23 1 (as 
shown in Sec. l5.3.3] l for the specific problem and observable 
considered. If one is, for instance, interested in predicting 
the effective bond length b for a weakly interacting system, 
a good trial value for bo should be the bond length / of non- 
interacting chains on the lattice [4 1 . To be consistent the per- 
turbation potential U must in this case vanish if the interac- 
tions are switched off (v = 0). If on the other side the aim is 
to characterize the deviations from Flory's ideality hypoth- 
esis in incompressible melts (^ ^ 1) one may naturally set 

^^ This is different in thin films where u^ » G^ decays only logarith- 
mically (94). 
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bo = b, i.e. one uses as reference the Gaussian chain which 
fits the chain on large scales (Sec. 12.4 .3T l. Obviously, in this 
case the perturbation potential U must vanish in this low- 
wavevector limit such that C^) w — {'^U)q — > 0. The last 
choice for bo turns out to be the best one in all cases where 
one does not need to predict the effective bond length b but 
where it can be determined independently by fitting some 
large scale intrachain property such as the typical chain end- 
to-end distance 7?n- 

The perturbation energy for a tagged test chain of length 
Ni is given by the potential 

rNi rk 

U= dk dlv{ni). (58) 

Jo Jo 

As further discussed in Sec. 12.4.21 the effective interaction 
v(r) between two monomers of the test chain arises due to 
the presence of the bath of surrounding chains which screens 
the direct excluded volume interaction v5(r). The calcula- 
tions are most readily performed in Fourier-Laplace space 
using definitions given in Sec. 12.2.31 See for instance the 
calculation presented in Appendix IB. 2l for the non-Gaussian 
contribution dG{r,s) = G{r,s) — GQ{r,s) to the subchain size 
distribution or Appendix IB .4 1 for the bond-bond correlation 
function Pi (i). 

2.4.2 Effective interaction potential 

We have still to specify the effective monomer interaction 
v{q) = J^[v(r)] in reciprocal space with q being the wavevec- 
tor conjugated to the distance between two monomers n and 
m of the tagged chain. Note that in general the test chain 
length A't and the mean chain length (A'^) of the bath may 
differ. Within linear response the effective pair interaction 
reads 111183 

-i- = i- + Fo(^). (59) 

v{q)p vp 

The first term stands here for the bare excluded volume in- 
teraction V between monomers. As we have already stressed 
above, Eq. ( 143b . thermodynamic consistency requires that v 
is set by the excess contribution to the isothermal compress- 
ibility of the solution: v = X/gp 1 18.83.86 J25j FI Fo{q) 
stands for the ideal chain intramolecular form factor of the 
bath of chains surrounding the reference chain. According to 
Eq. ( 138b the effective interaction v(^) depends thus in gen- 
eral on the length distribution /7n of the bath. We remind 
that for Flory-distributed melts the form factor is given by 
Eq. ( l40l) . Replacing /i by 2/N this corresponds to the Fade 
approximation, Eq. dJTJ) . of the awkward Debye function for 
monodisperse chains. 

Let us first assume that fo(^) ':^ S{q) Ki g, i.e. we assume 
^ <C 1/(7 in incompressible solutions {g <^1) and q <^ l/E, 
for systems with a well-defined thermal blob (g ^ 1). The 
effective interaction is then given by 

v{q)p ^TrTT^OTq<^l/^ and q<^l/a, (60) 

^^ The bulk modulus v only dominates v(^) for all q in extremely 
compressible systems where g » Fq (0) = {N^ ) / (A') . 



i.e. the effective interaction is given alone by the inverse 
structure factor of the bath and does not depend explicitly on 
the compressibility g of the solution^ According to Eq. ( 1391 ) 
the effective potential becomes in the low-wavevector limit 



vo = ^{q -^ 0) 






for^< l/Rg 



(61) 



i.e. Vo = 2/i/p for Flory-distributed melts and vq = 1/pA^ 
for monodisperse melts. Long test chains are ruled by vq 
which acts as a weak repulsive pseudo-potential with asso- 



ciated Fixman parameter z ^ vo^/M Ii4j||f_| It follows that 



z > 1 for TVt > {{Ny I {N)f w {NY 



(62) 



and the chains thus must swell and obey excluded volume 
statistics |3 7|13 We note that for a Flory-distributed bath 
Eq. (160b becomes 

v(?)p»x((flo^)^+M)for^<l/^ and^< Xja (63) 

which can also be used within the Fade approximation (/i = 
2lN) for the calculation of monodisperse systems f80''841. 
More importantly, the effective potential becomes for inter- 
mediate wavevectors 



v(^)p ; 



{a-mf 



forl//?g,z<^< 1/^ 



(64) 



and this irrespective of the length distribution p^ of the bath. 
Eq. ( 164b lies at the heart of the announced power-law swell- 
ing of (sub)chains, Eq. ^ or Eq. (|55]l I80l[83l [84ll87]l . 

For later reference in Sec. lB.4l we note that for a Flory- 
distributed bath of finite compressibility the pair potential 
reads 



v(^) =v- 



(fl0^)2+^ 



for ^ < 1 /a and jU^ < 1 (65) 



(ao^)2 + (ao/^] 

where |^ = a^gjl = a^j'^Mp following Eq. ( |44] |. Allowing 
to characterize wavevectors below and above 1 /^ , Eq. ( l65T l 
reduces to Eq. ( 161b for very low wavevectors and to Eq. ( 1641) 
in the intermediate wavevector range. In the limit of asymp- 
totically long chains (\i — > 0) Eq. (l65T l becomes lU 

q^ 
Hq) = V 2 !:.-2 f""" 1 /^s,^ <-q<-i/a (66) 



which corresponds in real space to ID 

exp(-r/^) 



v(r) = V I 5(r 



471^^ 



(67) 



^^ Since we shall set bg = b -dt the end of calculation and since the 
effective bond length b depends on g, the effective potential v{q) de- 
pends implicitly on g. 

^' Characterizing the excluded volume interaction free energy of a 
chain with itself the Fixman parameter of a chain of length A' of ex- 
cluded volume V may be defined more generally as 



z{N)- 



:v(iV/4)X = 



(v/fo'')A'2-''/2. 



Note that for a monodisperse melt of chain length N = g and di- 
mensionless compressibility g = 1/vp the Fixman parameter and the 
Ginzburg parameter become identical, ^(A' = g) = Gz- 
^^ For this reason an upper bound is indicated, e.g., in Eq. ( I106t . 
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i.e. the effective potential consists of a strongly repulsive 
part v5(r) of very short range (r w a), and an attractive part 
of range £, stemming from the compression of the reference 
chain by the bath chainsl3 Using Eq. ( 166b it follows that 



/ dr v(r) = v(^) 



= 



(68) 



which is commonly taken as a proof that "there is no ex- 
cluded volume interaction among the segments whose mean 
separation is larger than ^" [4|. Unfortunately, Eq. ( 168b does 
not imply mathematically that all other moments, say the in- 
tegral over v(r)r^'', should also rigorously vanish. It is thus 
incorrect to state that all possible correlation functions must 
be short-ranged. However, it remains relevant that v(^) ^ q^ 
vanishes with decreasing wavevector and the same applies 
for the total perturbation U{q) to the Gaussian reference^ 

2.4.3 Free energy for high compressibilities 

For later use in Sec. |4] we reformulate here a perturbation 
calculation result obtained long ago by Edwards lU using 
Eq. ( [66l) which allows to predict thermodynamic properties 
of melts with sufficiently large compressibilities g. Integrat- 
ing twice with respect to the density p the osmotic pres- 
sure given by Eq. (5.45) or Eq. (5 .11.5) of fH one obtains for 
monodisperse melts the free energy per monomer 

/3/(/3) = /3e,eif + ^log(paVA^) + ^v2p 



1 1 

127r|V 



with^^ = 



12v2p 



(69) 



and P = l/T being the inverse temperature, bo = I the bond 
length of the Gaussian reference chains and viiP) the sec- 
ond virial coefficient of a solution of unconnected mono- 
mers. The first term jSeseif is due to the (essentially constant) 
intrachain self-energy which shall be discussed in Sec. 14.21 It 
is due to the reference energy chosen in our numerical model 
Hamiltonian and it is normally not accessible experimen- 
tally. A similar intrachain energy contribution to the free en- 
ergy arises also from Eq. (5.43) of |4 1 if an upper cutoff qmax 
is introduced for the wavevectors q to avoid the ultra-violet 
divergence. Such an upper cutoff is justified by the discrete- 
ness of the monomers of real polymers. This leads necessar- 
ily to a non-universal free energy contribution which can be 
seen as an integration constant with respect to the integration 
of a measurable property such as the osmotic pressure or the 
compressibility. The second term in Eq. ( l69l) represents the 
translational invariance of monodisperse chains of length A^ 
(van't Hoff 's law). Due to this contribution the compressibil- 
ities depend in general on A'^ as will be discussed in Sec. l4.4F^ 
The (bare) excluded volume interaction between the mono- 
mers is accounted for by the third term. The underlined term 

^' See Fig. II.l of de Gennes' book [Sj. 

■'^ Hence, {"W) ^- in the large-scale limit for an observable "W which 
suggests bo = b for the Gaussian chain reference bond length. 

^' For general polydisperse melts of given partial densities p^v = 
PPn/ {N) the ideal gas contribution becomes Y,nPn^°Z{Pn)/P- 



in the second line represents the leading correction to the 
previous term due to the fact that the monomers are con- 
nected by bonds summing over the density fluctuations to 
quadratic order. As one expects |3 1, the corresponding cor- 
relations of the density fluctuations reduce the free energy 
by about one k-^T per thermal blob of volume ^^. Interest- 
ingly, according to Edwards H the chain connectivity, i.e. 
the presence of attractive forces between bonded monomers, 
does not change the excluded volume V2 — as one would 
expect naively — but rather gives rise to an additional term 
scaling differently with densitvFn Various thermodynamic 
properties are readily obtained from the quoted free energy 
and will be compared with our numerical results in Sec.|4] 
The underlined density fluctuation contribution to the free 
energy will be demonstrated numerically from the scaling 
of the specific heat cy (Sec. 14.3b . 

As the reader might have noticed we have written the 
free energy in Eq. ( 169b following Edwards assuming v = V2 
for the bare monomer interaction and bo = I for the bond 
length of the Gaussian reference chain. As already alluded to 
above (Sec. l2.4.Tb . one would nowadays rather set v = l/gp 
and bo = b using the imposed or measured dimensionless 
compressibility g and the measured effective bond length b. 
However, the choice of Edwards has a clear advantage: g 
and b may not be known with sufficient precision while the 
second virial coefficient V2 and the bond length / can always 
be calculated from the given model Hamiltonian. Accord- 
ing to Eq. (5.46) of Ref. |4| the stated free energy is sup- 
posed to hold in the limit where Gz ^ 1 with Gz ^ l/y^ 
being the Ginzburg parameter. As we shall see in Sec. |4] 
this restricts the validity of the related predictions to rather 
weak values of the (reduced) Lagrange multiplier e/T ap- 
plied to control the compressibility. In the range of validity 
of Eq. (|69] | it turns out that V2 w l/gp and I w b, i.e. the 
difference between both parameter choices correspond to ir- 
relevant higher order corrections. 

2.4.4 Subchain size distribution and its moments 

We turn now to the perturbation calculation predictions of 
intrachain conformational properties which will be compar- 
ed with our numerical data in Sec. |5] We focus first on the 
scale-free wavevector regime for arbitrarily long chains de- 
scribed by the effective interaction potential Eq. ( 1661) , i.e. 
effects related to the chain length distribution pN are irrele- 
vant. For the observable £/ = r^,„ with 1 <^n<m = n + s-^ 
Nt we indicate in Fig. [Sfa) the different interaction graphs 
one may compute in real space II80II 






-45 



n bt, 



I+=I- = 



VTAn^blp 



VTAn^blp 



blVs 



^* A free energy contribution » 1/^nP ^ 1/A'^" may be added to 
Eq. ( I69t if one insists on taking as reference for the connectivity con- 
tribution to the free energy the limit g — > A', i.e. ^ — > Rf^, where the 
chain connectivity becomes irrelevant. 
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Fig. 5 Sketch of relevant interaction diagrams for the perturbation cal- 
culation of a long test chain of length A't. The dashed lines indicate 
the effective monomer interaction v(r) between two monomers k and / 
on the chain whose distance is weighted using the Gaussian propaga- 
tor Go(n 1^ — l\)- (a) The calculation of the subchain size R^ requires 
the computation of four graphs, the dominant contribution stemming 
from the interactions of monomers within the subchain. The numerical 
factors indicate the relative weights contributing to the I/-/!? predicted 
by Eq. i fTot 1 80 84|. (b) The direct calculation of the bond-bond cor- 
relation function P\ (s) is simplified by choosing two bonds /j and (2 
outside the .s-segment. For symmetry reasons only the graph /q for the 
effective interactions between the monomers in the two dangling tails 
of lengths S\ and 52 gives a non- vanishing contribution. 



/o = 



3 



-bl^s. 



(70) 



\/7An^blp 

The diagram I{ corresponds to the standard graph computed 
by Edwards for the total chain {s = N^— 1) |4|. Consistent 
with Edwards its leading Gaussian contribution describes 
how the effective bond length is increased from bo to b under 
the influence of a small excluded volume interaction inside 
the subchain between n and m. Note that all other contri- 
butions proportional to ^/s correspond to the leading non- 
Gaussian corrections predicted in Eq. jlOl) . They only de- 
pend on bo and p but, more importantly, not on v in agree- 
ment with the scaling discussed in Sec. 12.3.61 The relative 
weights of these four contributions are indicated in Fig. [3 a) 
in units of —bQ^/s/{^/24n^bQp). The dominant correction 
stems from the interaction /[ within the subchain. The di- 
agrams /+ and /_ are obviously identical in the scale-free 
limitlfj Summing over all contributions this yields 

c^fb_ 

\fs \bo 
w here in the second line we have used the definition Cs '-- 
\j7ApK^ jpb^ already mentioned in Sec. ll.Sl and have set 



= b'-s 1 - 



(71) 



bi\\ 



12 v^ 

nb% 



= bi\\ 
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% 



Gzo 



(72) 



^^ Using Eq. J4) it follows that the interactions described by the 
strongest graph /; align the bonds /„ and /„, while the others tend to re- 
duce the effect 1 80 1. As shown in Fig.jSjb), it is better to plac e the b ond 
pair outsi de the subchain if one computes P\ [s) directly (Sec. 12.4.61 and 
Sec. lB.4t . For symmetry reasons only the interaction graph /q between 
the dangling ends matters for the alignment of the bond pair. Both pic- 
tures are consistent and lead to the same result. 



with Gzo = ^/yplb\p. 

For higher moments of the distribution G{r,s) it is con- 
venient to calculate first the perturbation deviations of the 
Fourier-Laplace transformation dG{q,t) = ^[^[dG{r,s)]] 
with dG{r,s) = G{r,s) — Go{r,s) and to obtain the moments 
from the coefficients of the expansion of this "generating 
function" in terms of the squared wavevector q^. As ex- 
plained in detail in the Appendix IB. 21 this leads to a devi- 
ation 



dG{r,s) = 



Inbls 



3/2 



exp 



3^\ _£^ /_fo_ 



2bh) VI v^ 



with n = r/bo\/s and the universal function 



^(") = V32 -;^ 



9« — -«' 



/(n)(73) 



(74) 



which allows to specify all moments of G(r,.s) 18411 . 

2.4.5 Adjusting the bond length of the reference chain 

The above perturbation result Eq. ( 1721) is of relevance to de- 
scribe the effect of a weak excluded volume v on a refer- 
ence system of ideal polymer melts with bond length bo = I 
where all interactions have been switched off (v = 0). It is 
expected to give a good estimation for the effective bond 
length b only for a small Ginzburg parameter Gz » Gzo <C 1 . 
For the dense incompressible melts we want to describe the 
latter condition does not hold and one cannot hope to find 
a good quantitative agreement with Eq. ( 172b . Note also that 
large wavevectors contribute strongly to the leading Gaus- 
sian term. The effective bond length b is, hence, strongly 
influenced by local and non-universal effects and is very dif- 
ficult to predict in general (Sec. 15.31) . 

Our more modest goal is to predict the coefficient of the 
l/v^-perturbation and to express it in terms of a suitable 
variational reference Hamiltonian characterized by a conve- 
niently chosen bo and the measured effective bond length b 
(instead of Eq. (|72l)). Following Refs. |43 ,84| we argue that 
for dense melts bo should be renormalized to b to take into 
account higher order graphsj^ Restating thus Eq. ( 1731) with 
bo = b the subchain size distribution may be rewritten 



SG{r,s) c, 
Go{r,s) s/s 

and for the 2/7-th moment of distribution this yields 



r^P' 



(2p+l)! 
6Pp\ 



ib's) 



?>{2Pp\pf c, 

2{2p+\)\^s 



(75) 



(76) 



^^ The general scaling argument discussed in Sec. 12.3.61 states that 
we have only one relevant length scale in this problem, the typical 
subchain size R^ ~ b^ itself. The incompressibility constraint cannot 
generate an additional scale. It is this size R^ which sets the strength of 
the effective interaction which then in turn feeds back to the deviations 
of i?s from Gaussianity. Having a bond length io in addition to the ef- 
fective bond length b associated with R^ would imply a second length 
scale foov^- 
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which reduces for p = 1 to Eq. jlQI) as stated in the Intro- 
duction. As a consequence the non-Gaussianity parameters 
Kp{s) and CXp{s) defined in Sec. l2.2.3] become 



and 

ap{s) = 



2{2p + l)\^s 



f3(2Pplpf_ \c^ 



(77) 



(78) 



'"'^^ \2{2p+lV. '^ V~^- 

Eq. (178b can be obtained from Eq. ( 177b by expanding the 
second moment (p = 1) in the denominator of the definition 
Eq. (|27ll. 

2.4.6 Bond-bond correlation function 

The bond-bond correlation function Pi (i) is a central prop- 
erty since it allows to probe directly the colored forces act- 
ing on the reference chain due to the incompressibility con- 
straint, Eq. ( 1176b . Using Eq. (|4]l Pi (s) may be obtained by 
differentiating the second moment R^ of the subchain size 
distribution G{r,s) with respect to the arc-length 5. For ar- 
bitrarily large chains and s ^ g this yields Pi (s) = cp/s^''^ 
with cp = CooCs/8 as announced in the Introduction. 

It is also possible to obtain Pi (s) directly by averag- 
ing the observable £/ = l^ -l^Jfi. Since for linear chains 
(^)o = by construction, the task is to compute Eq. ( I57brl 
Changing slightly the notations as indicated in Fig.[5jb), we 
consider two bonds Zj and I2 outside the .^-segment. The 
lengths of the two tails of the chains are denoted ^i and ^2. 
One of the tails, say S2 = Nt — s — S\ — \, may be fixed by 
the total length of the test chain. Placing the head of the first 
bond /) at the origin we consider the correlation function 
C\{r) = (ii(0) -hiir)) /i^ between two bonds separated by 
the distance r = \r\. Interestingly, it can be seen by sym- 
metry considerations that Ci (r) does only depend on the ef- 
fective interaction of the monomers in the first tail with the 
monomers in the second tail and not on the monomers in the 
intermediate strand of length s. Hence, we need only to cal- 
culate one interaction graph as opposed to the four graphs 
required by the calculation of Pi {s) throught R^. It is for this 
reason we have chosen the indicated positions of heads and 
tails of the bond vectors. With B{t) denoting the normalized 
distribution of the bond vector / of the polymer model con- 
sidered, the interaction graph in real space may be written 

Ci(r) = ~ j AU<M_M-x)B{l2) ^{lUi) 

X L £ /dridr2Go(ri+li,5i)Go(r2 -12,^2) 

si=0s2=0'' 

xv(r2+l-li) (79) 

where r^ points from the head of the bond li to the monomer 
si in the first tail and £2 from the tail of the bond I2 to 

^^ We remind that for closed cycles the ring closure implies long- 
range angular coiTclations even for Gaussian chains, hence for rings 
(sz/) » (^)q 7^ to leading order (87]. 



the monomer ^2 in the second dangling chain end. As the 
.s-segment is not implied in the perturbation of Ci(r), the 
constraint which consists in putting the two points on the 
same chain and putting a .s-strand between them introduces, 
to lowest order, the Gaussian propagator Go{r,s). Using Par- 
seval's theorem the bond-bond correlation function reads 



Pi{s) = I drCi{r)Go{r,s) = I -g-jCi{q)G{q,s) 



(80) 



with Ci (q) being the Fourier transform of Ci (r)|3 To sim- 
plify the notations we set from the start bo — b, i.e. we take 
the effective bond length as the bond length of the Gaussian 
reference chain. Although this is not strictly necessary, the 
calculation in reciprocal space may be strongly simplified 
by assuming the bond vector distribution to be a Gaussian 
nil) = Go{Ls=l) with b^ = 2da^ = f = {f). This implies 
that the chain is perfectly flexible, i.e. c„ = ib/lY = IrH 
Under these premises a bond vector / may be represented in 
reciprocal space as 



■^[IB{D] = idqB{q) w ia^ 2^ = - ^^ £ 



- d - 



(81) 



with S(^) = ^[B(/)] = exp(— (a^)^) being the Fourier trans- 
formed bond vector distribution Cj Let us denote the wave- 
vectors conjugated to the bonds Zj and I2 by q and q , re- 
spectively. The Fourier transform of the observable £/{li , I2) 
thus reads 



■11 '12 



^9J = -^£r^2- 



(82) 



Using Eq. ( I82l) for the observable and Eq. ( I66l) for the inter- 
action potential for Flory-distributed systems of given com- 
pressibility g one may integrate Eq. ( 180b in reciprocal space 
as shown in Appendix IB .41 In the limit of very long chains 
in — > 00) one obtains ind = 3 dimensions I86JI 



Pii^) = 



Cp 

p3/2 



Ay/2ne^''&ric{y/2ii)\ 



(83) 



as a function of the reduced arc-length u = s/g with erfc(ji:) 
being the complementary error function |44|. As one ex- 
pects, Eq. ( 1831) reduces to Eq. (|9]l for large u ^ I, i.e. irre- 
spective of the compressibility g the bond-bond correlation 

^* It can be shown that for infinite chains Ci (g) = A\{q) / [bq)-'- and, 
hence, 

C\ {r) = ^^^ exp(-r/.^) ind = 3, 

i.e. Ci {r) — > f or 1^ ^- at fixed distance r. Note that Ci (r) is a mere 
technical intermediate quantity which should not be conf used w ith 
the bond -bond con-elation function P\{r) discussed in Sec. 15.4.21 and 
Sec. 13] 

^' The complete formula for systems of general rigidity can be recov- 
ered by multiplying the final perturbation calculation result with c^ as 
may be seen by scaling considerations 1 87 89 1 or by simply comparing 
the result with the bond-bond correlations obtained using R^. 

■*" For a general bond vector distribution one may expand B[q) at low 
momentum as indicated in Sec. IB. II 
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function behaves as in the incompressible limit. In the op- 
posite limit where the structure within the thermal blob is 
probed Eq. (183b corresponds to the weaker decay 



P,{S): 



cp 



g'l^ ^' 



(84) 



This regime is consistent with the classical expansion result 
of the chain size in terms of the Fixman parameter z{s) w 
w^/b^ 1861 r^ We therefore refer to this limit as the "Fix- 
man regime". 

2.4. 7 Finite chain size ejfects 

To describe properly finite chain size corrections, Eq. ( 164b 
must be replaced by the general formula Eq. ( 159b . For mono- 
disperse chains {N — {N) — Ni) the form factor FQ{q) is given 
by Debye's function Eq. ( 136b . This approximation allows in 
principle to compute, e.g., the mean-squared total chain end- 
to-end distance, £/ = [r^^ — LiY- Ons verifies readily (see 
f4l, Eq. (5. III. 9)) that the effect of the perturbation may be 
expressed as 



(^)o(t/)o-(-«^t/)o = 



dq 



[In) 

/ As s^{N -s) ex-pi- {a^qY 
Jo 



3 v(^) 4{aoqfal 



(85) 



We take now first the integral over s. In the remaining in- 
tegral over q small q wavevectors contribute to the V^- 
swelling while large q renormalize the effective bond length 
of the dominant Gaussian behavior linear in A'^ (as discussed 
above). Since we wish to determine the non-Gaussian cor- 
rections, we focus on small wavevectors q <^ 1 /^ , i.e. the 
effective interaction potential is given by Eq. ( 1601) . We thus 
continue the calculation using v(A',x)p = 1/(A'^p/d(ji:)) with 
/d(x) being Debye's function and x = {Rg{N)q)^ = {aoq)^N. 
This allows us to express the swelling as 



I 



iilN 






(86) 



We have set here bo = b in agreement with the renormal- 
ization of the reference bond length discussed above. The 
numerical integral I{xu) = Jq" dx. . . over x is slowly conver- 
gent at infinity. As a consequence the estimate I{°°) = 1.59 
may be too large for moderate chain lengths. In practice, 
convergence is not achieved for values Xu{N) w {b/^)'^N 
corresponding to the screening length ^ . 

We remark finally that for various properties numeri- 
cal integration can be avoided replacing the Debye function 
by the Pade approximation, Eq. ( 1631) . This has been done 
for instance for the calculation of finite chain size effects 
for the bond-bond correlation function Pi {s,N) discussed in 
Sec. l5ATFI 

*' Omitting all prefactors we remind |4| that to leading order R^ rn 
b^s{l +z{s) ...). Using Eq. ^ it follows that Pi {s) ~ l/g/s. 

^^ It is interesting to compare the numerical value /(oo) w 1.59 ob- 
tained for the r.h.s. of Eq. ( [86l with the coefficients one would ob- 




Flg. 6 The BFM is an efficient lattice MC algorithm for coarse- 
grained polymer chains where monomers are represented by cubes on 
a simple cubic lattice (of lattice constant cr) connected by a set of al- 
lowed bond vectors, (a) The classical BFM assumes that lattice sites 
are at most occupied once. The panel shows the recently proposed 
variant with finite excluded volume penalty |86|. An energy e has to 
be paid if two cubes totally overlap. A corresponding fraction is as- 
sociated with the partial monomer overlap, as sketched for two cube 
corners, (b) Using local MC jump attempts to the next (filled circles) 
and next-nearest (open circles) neighbors we investigate in Sec. |6] the 
influence of the incompressibility constraint on the dynamics of over- 
damped polymer melts without topological interactions, (c) Using hard 
and structureless walls systems of reduced effectiv e dime nsion d <3 
may be investigated as outlined in Sec. l6.2l and Sec. 17.2.11 



3 Bond-fluctuation model 

3.1 Introduction 

The theoretical predictions sketched above should hold in 
any dense homopolymer solution assuming that the chains 
are asymptotically long, i.e. at least N/g ^ 10^ and even 
better N/g ^10^. The computational challenge is to equili- 
brate and to sample such configurations using an as simple 
as possible coarse-grained model for polymer melts I21'i29]. 
In this study we use the BFM, an efficient lattice MC algo- 
rithm proposed as an alternative to single-site self-avoiding 
walk models by Carmesin and Kremer in 1988 [28|. As il- 
lustrated in Fig. |6] the key idea of the model is to increase 
the size of the monomers which now occupy whole unit 
cells on a simple cubic lattice connected by a specified set 
of allowed bond vectors. While the multitude of possible 
bond lengths and angles allows a better representation of 
the continuous-space behavior of real polymer solutions and 
melts, the model remains sufficiently simple retaining thus 
the computational efficiency of lattice models without being 
plagued by ergodicity problems 1129 P . The BFM algorithm 
has been used for a huge range of problems addressing the 
generic behavior of long polymer chains of very different 
molecular architectures and geometries: statics and dynam- 

tain by computing Eq. ([85j either with the effective potential \{q) 
for infinite chains given by Eq. l l66t or with the Pade approximation, 
Eq. (|63). Within these approximations of the full linear response for- 
mula, Eq. l |59t . the coefficients can be obtained directly without nu- 
merical integration yielding overall similar values. In the first case we 
obtain 15/8 Ri 1.87 and in the second 11/8 » 1.37 1 84 1. While the first 
value is clearly not compatible with the measured end-to-end distances, 
the second yields a reasonable fit, especially for small A' < 1000. 
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ics of linear l50l|5T]|52l|53]|54l|55l|4g7T][80ll82l[83l[84l[85]l^ 
ll26yi27lll28lfT29lll30lll31l and cyclic IIJ2, 121 133 1 homopoly 
mer melts, polymer blends |[134"135','136''137,1381, gels and 
networks 1 1 39 1 , glass transition [ 1 40 , 1 4 1 1 42 1 43 , 1 44 1 , poly- 
mers and copolymers at surfaces H 14511 14611 . brushes |147, 
fT48]fT49l . thin films n30lll51lll52lfT53 T95^. equilibrium poly- 
mers 01O5II1O6II1O7II154I and other problems related to monomer 
and chain self-assembly fl55"1561. For recent reviews on 
the BFM algorithm see Refs. 129,301 . 

Throughout this paper all lengths and densities are given 
in units of the lattice constant a, time scales are given in 
units of the Monte Carlo Step (MCS) and Boltzmann's con- 
stant ^B is set to unity. Taking apart some paragraphs in 
Sec.|4]we assume a temperature T = I. If not specified oth- 
erwise the chains are monodisperse of length A'^. 

We define first the classical BFM variant without mono- 
mer overlap (e = °°) and explain then how dense configura- 
tions may be obtained using a mix of local and global MC 
moves (Sec. 13.31 ). The generalization of the BFM Hamilto- 
nian to finite monomer overlap penalties e is presented in 
Sec. 13.41 Finally, we turn in Sec. l3.5l to polydisperse equilib- 
rium polymer systems with annealed size distribution. 
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Fig. 7 Diffusion time Tn = Rf^/6Df^ vs. chain length A' for different 
versions of the BFM without monomer overlap. All data are for our 
standard volume fraction ^ = 0.5. The BFM version with topology 
conserving local "L06 moves" is represented by stars. All other data 
sets use topology violating local "L26 moves" 1 84 1. The L26 dynamics 
(diamonds) is essentially Rouse-like which allows the determination of 
Dn although the monomers have not yet moved over i?N for the largest 
chains considered. Additional "slithering snake" (SS) moves increase 
the efficiency of the algorithm by approximately an order of magnitude 
(squares, bold line). A power-law exponent 1.6 ±0.1 (dashed line) is 
found if "double bridging" (DB) moves are included. 



3.2 Classical BFM without monomer overlap 

The classical implementations of the BFM idea do not per- 
mit monomer overlap, i.e. each monomer occupies exclu- 
sively a unit cell of 2^ lattice sites on a (i-dimensional sim- 
ple cubic lattice Il28ll29ll30l . The fraction (/) of occupied lat- 
tice sites is thus = T^p with p being the (i-dimensional 
monomer number density. A widely used choice for the al- 
lowed bond vectors for the 3D variant (li = 3) of the BFM in- 
troduced by the Mainz condensed matter theory group around 
K. Binder 15011511 15205311126111271112811129111301 is given by 





where P stands for all the possible permutations and sign 
combinations of a lattice vector. This corresponds to 108 
different bond vectors Z of 5 possible bond lengths (2, v5, 
\/6, 3, \/T0) and 100 angles between consecutive bonds. 
The smallest 13 angles do not appear for the classical BFM 
because excluded volume forbids the sharp backfolding of 
bonds. If only local hopping moves to the 2d = 6 nearest 
neighbor sites are performed — called "L06 moves" |84| — 
this set of vectors ensures automatically that polymer chains 
cannot cross. (The corresponding "L04 moves" for the 2D 
variant of the BFM are represented in Fig. |6lb) by filled 
circles.) Topological constraints, e.g. in ring polymers [1321 
112111133] or polymer gels II139L hence are conservedr^ 

^-^ Following 1 50 51] we keep lists of the monomer positions in ab- 
solute space, their corresponding lattice positions and of the indices 
1 < (■ < 108 of the bond vectors connecting the monomers of the chains. 
Since the bond vector index can be encoded as a byte, this allows a 
rather compact storage of the configurations. Predefined tables allow 



Consequently, several authors report reptation-type dy- 
namics for chain lengths N^ N^f^ 10^ at a standard "melt" 
volume fraction (/) = 8p = 0.5 ||50ll5T]l56ll57l Fl As may be 
seen from the stars indicated in Fig.|7] the relaxation time 7n 
obtained using L06 moves |84| becomes similar to the rep- 
tation theory prediction Tn ^ A^^ (dash-dotted line). We have 
used here — as elsewhere if not stated otherwise — periodic 
simulation boxes of linear dimension L = 256 containing 
Wmon = pL^ = 2^" w 10^ monomers. This large system size 
allows to eliminate finite-size effects even for the longest 
chain lengths studied. The relaxation time 7n = R^/6Dj^ has 
been estimated here (for historical reasons) using the self- 
diffusion coefficient Z)n obtained either from the monomer 
MSD h{t) or the chain CM MSD hj^{t). Other operational 
definitions of 7^ exist 05111561 which lead to similar, verti- 
cally slightly shifted results. The last data point for A^ = 1024 
has to be taken with care as usual in computational as in ex- 
perimental studies 01570 . Being obtained by extrapolation 
using the expected shape of h]^{t)/t in log-linear coordi- 



the rapid verification of the excluded volume condition on the periodic 
lattice. Following Miiller 1 134 30] we use a Wigner-Seitz representa- 
tion of the cubic lattice where a cube is not represented by 8 entries on 
the lattice but just by one variable in the cube center. This variable can 
be a boolean if we are only interested in homopolymers or an integer 
if we deal with a mixture of different monomer types. Although the 
Wigner-Seitz representation of the BFM algorithm is about a factor 3 
slower than the original implementation, it has the advantage that the 
code becomes more compact and can be more readily adapted to the 
various polymer architectures or interaction potentials of interest. 
'*'' The BFM version by Shaffer 1 54 55] assumes a different bond vec- 
tor set which leads to a higher chain stiffness c„ and, hence, to smaller 
values of the entanglement length A 'e- N ote that if one applies a local 
stiffness potential, as discussed Sec. 12.21 one finds quite ge nerally that 
topological effects become more pronounced I52I121I133I . 
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Fig. 8 Two-dimensional sketch of the global MC moves used: (a) 
Monomers move collectively along the chain if a slithering snake move 
is performed. Effectively, this amounts to removing a monomer (the 
striped one on the left), connecting it to the other end of the chain and 
leaving the middle monomers unchanged. Therefore, density fluctu- 
ations and constraint release only occurs at the chain end. At higher 
volume fractions it is thus necessary to add local hopping moves, as 
shown at the top of the panel to maintain the computational efficiency 
[49 1 . (b) Connectivity altering double bridging (DB) moves are very 
useful at high densities allowing us to extend the accessible molecular 
mass up to A' = 8192 |84|. Since density fluctuations do not couple to 
DB moves, local moves again must be added. 



nates, it corresponds to a lower bound for T^. Note that the 
relaxation time appears thus to increase even more strongly 
with N than the standard reptation theory predicts [4J . We do 
not pursue this issue here (which has also been observed in 
MD simulations) the important point being merely that local 
topology conserving moves are too inefficient to equilibrate 
and sample large -A^ polymer melts. 

Note that the classical BFM without monomer overlap 
and using L06 moves is strictly speaking not ergodic, since 
some configurations may be easily constructed which are not 
accessible starting from an initial configuration of stretched 
linear chains. Although topology conservation is irrelevant 
for the present work (taking apart the preliminary results 
presented in Sec. 16.51) we keep the set of allowed bond vec- 
tors, Eq. ( 1871) . for consistency with previous work. 



MSD hn{t) even for our largest chains with TV = 8192. As 
shown in Fig.|7] we find 7n ^ 530N^ for L26 dynamics. This 
is, of course, still prohibitive for the longest chains we aim 
to characterize. 

L26 moves yield configurations not accessible with L06 
moves. Concerning the static properties both system classes 
are practically equivalent. This has been confirmed by count- 
ing the number of monomers which become "blocked" (in 
absolute space or with respect to an initial group of neigh- 
bor monomers) once one returns to the L06 scheme. Typi- 
cally, we find about 10 blocked monomers for a system of 
2^*^ monomers II84L While the few blocked monomers are 
irrelevant for static properties they obviously matter if dy- 
namical properties are probed r^ 

3.3.2 Slithering snake moves 

In addition to these local moves one slithering snake move 
per chain is attempted on average per MCS corresponding 
to the displacement of A'^ monomers along the chain back- 
bone, as sketched in Fig. Ufa). Note that in our units two 
displacement attempts per MCS are performed on average 
per monomer, one for a local move and one for a snake 
move™ Interestingly, a significantly larger slithering snake 
attempt frequency would not be useful since the relaxation 
time of snakes without or only few local moves increases 
exponentially with mass as shown in 14811491 due to the cor- 
related motion of snakes expected in analogy to the acti- 
vated reptation limit for real polymer melts mentioned in 
Sec. 11.41 II46II47II . In order to obtain an efficient free snake 
diffusion [with an A'-independent curvilinear diffusion coef- 
ficient D^{N) - N° and Tn w N^/D^{N) - N^ 14911331 it is 
important to relax density fluctuations rapidly by local dy- 
namical pathways. As shown in Fig. |7] (squares), we find 
a much reduced relaxation time 7n w 40A^^ which is, how- 
ever, still inconveniently large for our longest chains. Note 
that most of the CPU time is used by the local moves and the 
computational load per MCS remains A'^-independent. 



3.3 Local and global topology violating MC moves 

To equilibrate and sample BFM melts such as the ones pre- 
sented in Fig. 121 we have replaced the realistic but very slow 
L06 moves by a mix of local topology violating so-called 
"L26 moves" §^ with global "slithering snake" 132113311491 
and "double bridging" 1361137117311291 MC moves as shown 
in Fig. El 

3.3.1 Local L26 moves 

Already the use of local moves to the 3'^ — 1 = 26 next and 
next-nearest lattice sites surrounding the current monomer 
position [Fig. Ob)] dramatically speeds up the relaxation dy- 
namics, especially for A^ > 512, as can be seen from the di- 
amonds indicated in Fig.JTl Since the dynamics is to leading 
order of Rouse-type, as further discussed in Sec.|6l the dif- 
fusion coefficient Dn can readily be estimated from the CM 



3.3.3 Double bridging moves 

Double bridging (DB) moves are found to be very useful at 
high densities and help us to extend the accessible molecu- 
lar masses close to A^ w 10^. As for slithering snake moves 
we use all 108 bond vectors to switch chain segments be- 
tween two different chains. Only chain segments of equal 
length are swapped to conserve monodispersity. Topological 
constraints are again deliberately violated. Since more than 
one swap partner is possible for a selected first monomer, 
delicate detailed balance questions arise II29II . This is par- 
ticularly important for short chains. To avoid the compu- 
tation of weights [29 1 one simple solution to this problem 

'^^ The same problems arise if slithering snake or double bridging 
moves are used to equilibrate the configurations. 

*^ It is computationally more efficient for large A' to take off a 
monomer at one chain end and to paste it at the other leaving all other 
monomers unaltered. Before static or dynamical measurements are per- 
formed the original order of the monomers must then be restored. 
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Table 1 Various properties for monodisperse BFM melts at volume 
fraction (j) = 0.5: the dimensionless compressibility g, the root-mean- 
square bond length /, the effective bond length h, the acceptance rate 
A, the local monomer mobility W and the self-diffusion coefficient Dn. 
The dynamical data refer to local L26 moves to the nearest and next- 
nearest lattice sites which yields essentially Rouse-like dynamics. The 
values of the dimensionless compressibility g and the effective bond 
length h for asymptotically long chains have been obtained using ex- 
trapolation schemes discussed in Sec.|4]and Sec.[5] respectively. 



is to refuse all moves with more than one swap partnerri 
The configurations are screened with a frequency /db for 
possible DB moves where we scan in random order over 
the monomers. The frequency should not be too large to 
avoid (more or less) immediate back swaps and monomers 
should move at least over a couple of lattice constants be- 
tween subsequent DB moves. In the example presented in 
Fig-El(spheres) a frequency /db =0.1 is used. Empirically it 
is found that 7n w l3N^^^^^■^ using the diffusion coefficient 
Dn obtained from the monomer MSD h{t). For N = 8192 
this corresponds to 3 x 10^ MCS. This allowed us even for 
the largest chains to observe monomer diffusion over several 
7?N within the 10^ MCS which were feasible in 2007 on our 
XEON-PC processor cluster |[84ll . The influence of /db on 
the performance has yet not been explored systematically, 
but preliminary results suggest smaller DB frequencies for 
future studies. The power-law exponent w 3/2 remains ro- 
bust as is also confirmed by MD simulations of a bead-spring 
model coupled to DB moves ll84l . This finding begs for a 
systematic theoretical investigation. 



3.3.4 Summary of static properties 

Some static properties obtained at our reference volume frac- 
tion (j) = 0.5 assuming no monomer overlap (£ = oo) are indi- 
cated in Table [U Averages are performed over all chains and 
typically 1000 configurations. Taking apart the systems for 
A' = 8192, chains are always much smaller than the linear 
box size L = 256. In the large A^-limit we obtain an average 

''^ Since we have a finite number of possible BFM bonds an even 
simpler option is to select randomly one bond vector / for a given 
first monomer at a lattice position r and to check whether a suitable 
monomer of another chain exists at r -I- /. Obviously, using only one 
instead of 108 bonds reduces the number of DB moves performed, but 
since correlated moves forth and back in topology space are pointless, 
this is not of disadvantage. 



bond length (|/|) w 2.604, a root -mean-squared bond length 
I w 2.636 and an effective bond length b « 3.24 as will be 
shown below in Sec. 15.31 This corresponds to a ratio c„ = 
{b/lY ~ 1 -52 and, hence, to a persistence length /p = /(c<x> + 
l)/2 w 3.32. The swelling coefficient Cs = ^/TAJk^ /b'^p de- 
fined in Sec. ll.5l is thus Cs w 0.41. Especially, we find from 
the zero wavevector limit of the total structure factor S{q) 
a dimensionless compressibility g w 0.246 [Eq. (143 H which 
compares well with real experimental melts. From the mea- 
sured bulk compression modulus v = \/g{p)p w 66 and the 
effective bond length b one estimates a Ginzburg parameter 
Gz = ^J^jb^p w 0.96. Following Ref. [18] the interaction 
parameter v is supposed here to be given by the full inverse 
compressibility and not just by the second virial coefficient 
V2 = 27 of the BFM monomers [J 27 J . 



3.4 BFM with finite excluded volume penalty 

3.4.1 Definition of Hamiltonian 

Fig. |6] shows how finite energy penalties may be introduced 
in the BFM algorithm [861. The overlap of two cube corners 
on one lattice site (A^ov = 1) corresponds to an energy cost 
of e/8, the full overlap of two monomers (A^ov = 8) to an 
energy £. More generally, with A'ov being the total number 
of interacting cube corners the total interaction energy of a 
configuration is 



E = ItVo, 



(88) 



With the energies of the final (Ef) and the initial configura- 
tions (Ei) we accept a proposed MC move according to the 
Metropolis criterion with a probability min(l,exp[— (Ef — 
E\)/T]) II26II27L If the overlap penalty is the only energy 
scale as in the studies presented here, one may, of course, 
either vary the overlap parameter e or the temperature T. 
For the presentation of thermodynamic properties in Sec.|4]it 
will be more naturally to use T as the control parameter and 
to fix arbitrarily £ = 1 . The inverse temperature fi = 1/T and 
the dimensionless overlap strength x = e/T are thus numeri- 
cally equal. (Both notations are kept for dimensional reasons 
and for future generalization to models with more than one 
energy scale.) In other parts of this review, especially Sec.|5] 
and Sec.|6] it will be more natural to set temperature to unity, 
T = 1, using the overlap strength e as the control parameter. 

3.4.2 Second virial coefficient 

To illustrate this finite excluded volume interaction we in- 
dicate the second virial of an imperfect gas of unconnected 
monomers, V2 = J d5_{l — e"^®' ^), which is shown below 
to be useful for roughly characterizing the effective strength 
of the potential. 5 stands for a possible lattice vector be- 
tween the centers of two interacting cubes. It is easy to see 
that there are 8 vectors corresponding to A^ov = 1 as shown 
in Fig. ISa), 12 to A^ov = 2 (overlap of two cube corners). 
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6 to A'ov = 4 (overlap of two faces), and 1 to A'ov = 8 (full 
overlap). Settings = e/T this leads to a second virial 

V2(jc) = 8x(l-exp(-jc/8)) 
+ 12x (l-exp(-jc/4)) 
+ 6x (l-exp(-jc/2)) 
+ lx(l-exp(-jc)) (89) 

given in units of the lattice cube volume a^ . We note that the 
second virial becomes constant, V2 = 27, for large x ^ 1 as 
expected |127|. This second virial coefficient is about half 
the effective bulk modulus v w 66 indicated above. In the 
opposite limit we have 



27 

V2(x) w 8x -x^ for X <C 1. 

16 

We shall see in Sec. l4.4l that V2 w v in this limit. 



(90) 



3.4.3 Implementation 

Since a lattice site may be occupied now by more than one 
monomer, it is not possible to use a compact boolean oc- 
cupation lattice as for the classical BFM. Instead we have 
mapped Eq. ( [88]) onto a Potts spin model i26ll 



(91) 



L d 



with constant monomer number n^^a = Hr^il) = L^P- We 
use the Wigner-Seitz representation of the BFM on the cubic 
lattice 034 30 1 where an integer spin variable S{r) counts 
the number of BFM monomers (5 = 0, 1 , 2, . . .) with cubes 
centered at a Wigner-Seitz lattice position r. Since we have 
now to compute the interaction between cube centers instead 
of cube corners, the coupling constant J characterizing the 
interaction between two spins depends only on the relative 
distance 5: 



J{5) = £ < 



1/8 if5^ = P(l,l,l) for cube corners, 
1 /4 if 5 = P( 1 , 1 , 0) for cube edges, 
1 /2 if 5 = P( 1 , 0, 0) for cube faces, 
1 if 5 = P(0, 0, 0) for full overlap, 
otherwise. 



(92) 



Since the interaction is still short-ranged and the values of 
J are readily tabulated, this remains an efficient rendering 
of the monomer interactions. Note that the first term on the 
r.h.s. of Eq. ( |9T1) contains a constant self-interaction contri- 
bution of the njnon monomers with themselves for 5^ = 0, 
which is subtracted by the second term™ 



'^^ Attractive interactions similar to the ones used in I134I135I136I 
11371 may be easily added to the Potts spin formulation of the soft BFM. 
The simulation of polymer blends requires additional Potts spin latti ces 
as the two lattices used to obtain the chemical potential in Sec. 14.51 
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Table 2 BFM solutions with overlap penalty e = 10 for different vol- 
ume fractions (^ = 8p. The indicated dynamical properties — accep- 
tance rate A, local mobility W and self-diffusion constant Dn — have 
been obtained using local L26-moves for chains of length N = 1024. 



3.4.4 Equilibration and system properties 

As start configurations we have used the equilibrated BFM 
configurations without monomer overlap (e = °°) described 
in Sec. 13.21 1|84|1 . As one may expect, the configurational 
properties are found essentially unchanged for x ^ 5 (Ta- 
ble[T|. Local L26-moves need to be added to the snake moves 
for X > 1. Otherwise the slithering snake motion will be- 
come ineffective [49|. Simple slithering snakes without lo- 
cal moves are sufficient, however, for smaller penalties. We 
have crosschecked our results in this regime for A^ = 2048 
and A^ = 8192 using boxes of linear size L = 512 by starting 
our simulations with Gaussian chains at x = and increas- 
ing then the penalty. Table [T]present some system properties 
obtained for our reference volume fraction (/) = 0.5 such as 
dimensionless compressibility g or the effective bond length 
b. Averages are performed over all chains and at least 100 
configurations. The chain lengths A'^ = 64, A^ = 1024 and 
N = 2048 have been studied with particular care. Density 
effects have been studied more briefly. As will be discussed 
in Sec. 14.21 we have sampled weak overlap penalties (x ^ 1) 
for A^ = 8192 to investigate the intrachain contributions to 
the mean energy. We have also probed various densities for 
A^ = 1024 and £ = 10 as summarized in Table |2l This was 
done to check for density effects on the deviations to Rouse 
dynamics as discussed in Sec. 16.41 



3.5 BFM with annealed mass distribution 

3.5.1 Motivation and context 

As discussed above (Fig. [T]), the equilibration and sampling 
of strictly monodisperse polymer melts is a delicate issue. 
An elegant way to test the computed conformational proper- 
ties is given by associating a finite "scission energy" £ > to 
each bond which has to be paid if a bond is broken II17II105I 
106 107 108 , 109 154|. Since we are only interested here in 
linear polymer melts, the formation of closed cycles and the 
branching of the chains is not permittedri As sketched in 
Fig. |9] we relax thus the constraint that bonds can never 
break, i.e. that the connectivity matrix (defining which monomers 
are connected by BFM bonds) is quenched, and allow the 

'^^ In systems of experimental relevance cl osed cycles ar e suppressed 
by the non-negligible rigidity of the chains I158I159I160I . 
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Fig. 9 Self-assembled linear EP: (a) Bead-spring representation of 
a worm-like micelle assuming a finite scission energy E and a fi- 
nite barrier height B 1 17|. The scission energy determines the static 
propeities and fixes the ratio of the scission and recombination rates, 
kg ~ exp{—{E+B)/k^T) and k,- ~ exp{—B/k^T). Both energy scales 
are sketched vs. a generic reaction coordinate q |117|. The formation 
of closed rings and the branching of chains are not allowed in the pre- 
sented studies, (b) Two-dimensional projection of EP modeled using 
the classical BFM without monomer overlap |107|. Chains consists 
of symmetrically connected lists of bonds. The pointers of end-bonds 
point to NIL. The breaking of a saturated bond ibond requires to set the 
pointers of the two connected bonds iboiul and jhond = pointer(ibond) 
to NIL. Setting pointer(-2)=4 and pointer(4)=-2 connects the two end- 
monomers imon=2 and jmon=4. 



polymerization of the chains and their respective monomers 
to take place under condition of chemical equilibrium. Such 
systems of self-assembled EP are not only useful for com- 
putational purposes but are also of high experimental rele- 
vancelfj An important example sketched in Fig. |9ja) is that 
of some surfactant molecules forming long giant worm-like 
micelles which break and recombine constantly at random 
points along the sequence [17 20|. Similar systems of EP are 
formed by liquid sulfur [il02iil03 l, selenium Iil04i and some 
protein filaments j9l- Although with respect to their static 
properties EP behave very much like quenched polymers, 
the constant reorganization of the chain connectivity offers 
an additional relaxation pathway reducing strongly the re- 
laxation times 01711104 inMHT09l . Obviously, EP are intrin- 
sically polydisperse with an annealed length distribution /7n 
minimizing the free energy of the system (Sec. 12.21) . Since 
both the bonding energy per chain —E{N — 1) and all other 
free energy contributions to the chain chemical potential fj.^ 
are extensive with repect to the chain length A' — at least ac- 
cording to Flory's ideality hypothesis — one expects a Flory 
distribution decreasing exponentially with chain length A'. 



^^ In the surfactant literature fTT] EP are often referred to as "liv- 
ing polymers" (LP) although this is potentially confusing since they 
are distinct from systems that polymerize stepwise, in the presence of 
a fixed number of initiators, for which this term has previously been 
reserved |,16ril. Since LP are held together by strong covalent carbon- 
to-carbon bonds, they do not break in the middle of the polymer chain. 



3.5.2 Spatial monomer moves 

The EP systems presented in this study have been obtained 
with the classical BFM algorithm without monomer over- 
lap (e = oo, r = 1) at the standard melt volume fraction 
= 8p = 0.5. Using again the Wigner-Seitz representation 
we mark (only) the index of each monomer on the peri- 
odic lattice of linear size L = 256. Thus the indices of the 
monomers in the neighborhood of a reference monomer are 
readily obtained which is helpful for the recombination of 
bonds described below. Only local L06 or L26 monomer 
moves have been used since the breaking and recombina- 
tion of the chains reduce the relaxation times dramatically 
compared to monodisperse systems 1 109 1. Additional global 
MC moves as described in Sec. l3.3l mav be added, however, 
in future studies. 

3.5.3 Connectivity pointer list 

Self-assembled EP are only transient objects and it is thus 



inefficient to base the data structure on the chains 11051 . 
rather it should be based on the (saturated or unsaturated) 
bonds of each monomer II107I . As sketched in Fig.JQjb), this 
allows via a linear pointer list between the bonds to avoid 
all sorting procedures. Using the assumption that no branch- 
ing of chains is allowed, the two (possible) bonds of each 
monomer imon are called ibond = imon and ibond = —imon. 
No specific meaning (or direction) is attached to the sign: 
this is merely a convenience for finding the monomer from 
the bond list: imon = \ibond\. Pointers are taken to couple 
independently of sign and the bonds are coupled by means 
of a pointer list in a completely transitive fashion. Only two 
simple operations are thus required for breaking or recom- 
bining bonds. Unsaturated bonds at chain ends point to NIL. 
Only these bonds may recombine. A minor caveat attached 
to this data structure arises if the ends of a given chain are not 
allowed to bind together as in the presented studies. Since 
there is no direct chain information in the data structure we 
have to check this constraint before every recombination by 
working up the pointer list which only adds four lines to the 
source codepH 

3.5.4 Connectivity altering moves 

As sketched in Fig.|9ja), EP systems are not only character- 
ized by the monomer density p and the finite scission en- 
ergy E which determine the static properties but also by a 
barrier height B >0 which only influences the scission and 
recombination rates. This barrier is taken into account by 
setting an attempt frequency (Ob = exp(— B/T) for choosing 
randomly one bond ibond out of the In^non bonds of the sys- 
tem. This frequency is a convenient tool for testing the dy- 
namics of EP at different lifetimes of the chains 1 108, 109] 
162] although for the static properties discussed the choice 

^ ' For higher E the simulation becomes actually faster per MCS since 
the number of recombinations goes down like the squared density of 
unsaturated end monomers (p/ (A'))^ ~ exp{—E) 11071 . 
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Table 3 Various properties of EP obtained by means of the 3D BFM 
algorithm without monomer overlap {£ = oo) at volume fraction (j> = 
8p = 0.5: imposed scission energy E, the mean chain length {N}, the 
ratio F(0) = ^A'^) / (A^) comparing the first and the second moment of 
the number distribution, the root-mean-squared bond length /, the z- 
averaged end-to-end distance R^.z and radius of gyration i?g,z obtained 
using Eq. l |19t with p = 2. For all scission energies we have used pe- 
riodic simulation boxes of linear size L = 256 containing «,„„„ = 2^" 
monomers. 




4 5 6 7 
x=N/<N> 

Fig. 10 Normalized chain length distribution /?n of linear EP for dif- 
ferent scission energies E obtained using the classical BFM algo- 
rithm without monomer overlap: (a) The main panel demonstrates 
the collapse of the rescaled distribution y = pN (A') as a function of 
X = N/ (N). The exponential decay (solid line) implied by Flory's ide- 
ality hypothesis is (to first order) consistent with our data, (b) First four 
moments of the distribution vs. E. (c) Replot of the data of panel (a) 
in log-linear coordinates focusing on short chains. The data points are 
systematically below the exponential decay (solid line) for x ^ 1. 



of B is irrelevant. The bond ibond corresponds to a monomer 
imon = \ibond\ at a position r. Depending on whether the 
bond ibond is saturated or unsaturated we try to break it or 
to connect it to a suitable nearby unsaturated monomer. A 
delicate detailed balance problem arises [ 109] if «„ > 1 un- 
saturated monomers are available for recombination (the no- 
closed-cycle condition having been verified). If one chooses 
now one of these monomers at random, a weight l/n„ has 
to be taken into account for the reverse breaking process. 
Choosing in addition to the reference bond ibond a trial bond 
vector / allows to avoid these weights. If ibond is unsatu- 
rated one searches for an unsaturated monomer only at the 
position r + / (taking into account the periodicity of the lat- 
tice). Since monomer overlap is forbidden there is at most 
one unsaturated monomer at this position. If this is case and 
if no closed cycle is formed, the recombination is accepted 
since the energy change is —E < 0. To satisfy detailed bal- 
ance a saturated bond can therefore only be broken if its 
bond vector is identical to the trial bond L Applying the 
Metropolis algorithm [26] a scission is performed whenever 
the value of a random number between and 1 is smaller 
than exp{—E/T). The fact that we only probe one lattice 
vector I for possible recombinations and not all possible 108 
obviously strongly reduces the number of recombination and 
scission events. However, since we are interested in uncor- 
related changes of the connectivity list, a broken monomer 
must anyway move over a certain distance, say 10a, before 
a new recombination attempt is made. Otherwise there is a 
strong chance that the newly created chain end monomer re- 
combines with its previous partner II108II1091I162II . 



3.5.5 Some computational results 

As summarized in Table [3] we have sampled EP systems 
with scission energies up to £ = 15, the largest energy cor- 
responding to a mean chain length (N) w 6183. As for the 
athermal classical monodisperse BFM systems we obtain a 
dimensionless compressibility g = 0.24, an effective bond 
length b = 3.244 and a swelling coefficient Cs = 0.41. As 
discussed in Sec. 12.2.21 one expects EP melts to be Flory 
distributed [Eq. (I20b 1 if the chemical potential of the chains 
is extensive with respect to their mass, fi^ ^ N. The main 
panel of Fig.[TO]presents the normalized length distribution 
/?N for different scission energies E as indicated. A nice 
data collapse is apparently obtained if /^n (^) is plotted as 
a function of the reduced chain length x = N / {N) using the 
measured mean chain length (A^) . At first sight, there is no 
sign of deviation from the exponential decay indicated by 
the solid line. The mean chain length itself is given in panel 
(b) as a function of E together with some higher moments 
{NP) = Y.N^'^Pt^ of the distribution. As indicated by the 
dashed line, we find (A^'^) ^ sxp{pE/2) as expected from 
Eq. (123 b [17 107 1 . The data presented in the first two panels 
of Fig. [To] is thus fully consistent with older computational 
work [ 105 106 107 108 109 llOlllllI which has let us to 
believe (incorrectly) that Flory's ideality hypothesis must 
hold rigorously. Closer inspection of the histograms reveals, 
however, deviations for small jc ^ 1 . As can be seen from 
panel (c), the probability for short chains is reduced with re- 
spect to the Flory distribution (solid line). We shall further 
investigate this depletion in Sec. 14.61 

We note finally that all EP systems presented here have 
been sampled within 4 months while the sample of monodis- 
perse configurations for N = 8192 alone required about 3 
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Fig. 11 Mean overlap energy per monomer y = e/e vs. x = e/T 
for several chain lengths A' as indicated. The energy decreases 
monotonously with increasing x. The decay becomes Arrhenius-like 
forx » 10 (dash-dotted line). The dashed line indicates the energy pre- 
dicted from the second virial of soft BFM beads, Eq. <93t . The main 
figure demonstrates the weak chain length dependence on logarithmic 
scales, especially for strong excluded volume interactions (x > 1). In- 
set: Same data plotted with linear vertical axis emphasizing the higher 
mean energy for long polymers (N > 64) for x <?; 1 caused by a self- 
energy contribution Cseif/E Ri 0.18. The self-energies are indicated by 
the filled triangles. Eq. {95) is represented by the bold line. 



years on a similar XEON processor. EP are therefore very 
interesting from the computational point of view, allowing 
for an efficient test of theoretical predictions. 



4 Thermodynamic properties of BFM melts 

4.1 Introduction 

To characterize the soft BFM model introduced in Sec. 13.41 
we will first investigate thermodynamic properties such as 
the mean overlap energy per monomer e, the specific heat 
Cv, the dimensionless compressibility g or the excess chain 
chemical potential /Jn as functions of the reduced overlap 
strength x = e/T = ej3. For the small-x limit these proper- 
ties have been calculated long ago by Edwards |4| as sum- 
marized in Eq. (169 l l. Various thermodynamic properties ob- 
tained from the quoted free energy will be compared with 
our numerical results |86|. To demonstrate that deviations 
from Flory's ideality hypothesis are also present in thermo- 
dynamic properties we will investigate in detail in Sec. 14.61 
the chemical potential in systems of annealed EP using the 
classical BFM algorithm without monomer overlap. 



4.2 Mean overlap energy 

From the numerical point of view the simplest thermody- 
namic property to be investigated here is the mean interac- 
tion energy per monomer, e = (E) /n^on, due to the Hamil- 
tonian, Eq. dSSl l. Fig. [TT] presents the dimensionless energy 
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Fig. 12 Reduced mean energy e/e (spheres) and self-energy eseif/£ 
(triangles) as functions of the number density pforA' = 8192,L = 512 
and X = 0.001. As shown by the dashed line, e(p) is a superposition 
of the mean field energy 4p and the (essentially) constant self-energy 
^seif/£ ~ 0.18A'''x°p''. Inset: e/e — 4p as a function of chain length 
\/\/N — 1 for our reference density p = 0.5/8 and for a single chain 
(p = 0). The linear slope (bold line) is expected from the return prob- 
ability of Gaussian chains. 



y — e/e for BFM melts ((/> = 0.5) for different chain length 
A'^. Decreasing the overlap strength x starting from config- 
urations obtained using the classical BFM (Sec. 13.31) , the 
interaction energy increases first exponentially for large x 
and levels off for x ^ \ where the monomers freely over- 
lap. The data for unconnected beads (A'^ = 1) represented by 
the filled spheres and polymer chains (A' ^ 1) are broadly 
speaking similar, especially for large overlap penalties, x > 
1 . Interestingly, the mean energy of polymer melts increases 
more strongly for x ^ 1 as can be seen better from the log- 
linear data representation chosen in the inset of Fig.[TT] Also 
shown in the inset is the mean intrachain self-energy per 
monomer Cseif (filled triangles) obtained for the largest chain 
length available for a given x. In fact about half of the energy 
of polymer melts for all x is due to the self-interactions of 
the chains 1.86.1 . For x^l the self-energy becomes eseif/£ ~ 
0.18 which is exactly the observed energy difference be- 
tween polymer and bead systems. 

Before addressing this point let us consider the energy of 
soft BFM beads (A^ = 1) for which the second virial coeffi- 
cient V2{x) has been given in Eq. ( |89] l. Since e = da {^f{^)) 
the mean energy becomes to leading order 11631 
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corresponding to the first term in the third line of Eq. 
Eq. ( |93l l is represented by the dashed line in Fig. [TT] It cor- 
responds to an Arrhenius behavior with .V w p exp(— ji;/8)/2 
for jc > 1 (dash-dotted line) and to y — > i8p = 4p for x < 1 . 
This simple formula predicts well the bead data over the en- 
tire range of x and also yields a remarkable fit for polymer 
chains with larger overlap penalties. 
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The energy difference between polymer chains and beads 
for X <C 1 is accounted for by the first free energy contribu- 
tion indicated in Eq. ( |69l l. This contribution is further inves- 
tigated in Fig. [T2] presenting data for such small x that the 
entropy dominates all conformational properties. The self- 
energy of a chain is thus given by the probability p{s,d_) that 
a random walk of 5 BFM bonds returns to a relative position 
5 with respect to a reference monomer at r. Hence, 
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where the first sum runs over all positions with non- vanish- 
ing coupling constant J(5) as defined in Eq. ( l92l ). The proba- 
bility p{s,d_) and the weights J{8_)p{s,8_) can be tabulated in 
principle for small s. Since the return probability decreases 
strongly with s, these model-specific small-5 values dom- 
inate the integral, Eq. ( 194b . As can be seen from the in- 
set of Fig. [12] for single chains (p = 0), eseif ~ 0.18e for 
large A'. The weak A'^-dependence visible in the panel stems 
from the upper integration boundary over the Gaussian re- 
turn probability which leads to a chain length correction lin- 
ear in f = l/\/N — 1 (bold line). Also shown in the panel 
are energies for our reference density p = 0.5/8. They are 
shifted vertically by the mean field energy 4p assuming that 
density fluctuations of different chains do not couple. The 
main panel presents the mean energy e and the mean self- 
energy fiseif as functions of the density p for chains of length 
A^ = 8192. The self-energy (triangles) stays essentially p- 
independent. The total interaction energy sums over the self- 
energy and mean-field energy contributions as shown by the 
dashed line. The self-energy contribution can only be ne- 
glected for volume fractions larger than unity. 

Summarizing Eqs. ( |69l l and ( 1901 ) the energy should scale 
to leading order in x as 



};w0.18+4p- 



24^/2 ,Jxp 



Pp 



for X < 1 



(95) 



where the two x-independent contributions have already been 
discussed above. The underlined term stems from the den- 
sity fluctuation contribution in Eq. ( |69l l. Eq. ( 1951 ) is indicated 
by the bold line in the inset of Fig.[TT] It yields a reasonable 
description for small x. Since the energy is dominated by 
the two constant contributions to Eq. ( 195b for x < 0.001 and 
since higher expansion terms become relevant for x > 0.1, 
the predicted y^-decay corresponds unfortunately only to 
the small-x regime. To show that it is indeed the density fluc- 
tuation term which dominates the temperature dependence 
for X ^ 1 we will consider now the specific heat cy, i.e. the 
second derivative of the free energy with respect to j3 . 



4.3 Energy fluctuations 

The fluctuations of the interaction energy are addressed in 
Fig. [13] displaying the enthalpic contribution to the specific 




Fig. 13 Specific heat per bead cy vs. x. The dashed hne indicates the 
energy fluctuations predicted from the second virial, Eq. j96t . which 
fits the data of soft BFM beads {N = \) over six decades. While 
the chain length does not matter for strong excluded volume interac- 
tions, the energy fluctuations are found to increase strongly with N for 
-t <g; 1. For short chains we observe cy ~ pN^/^x^ as can be seen for 
A' = 16 (thin solid line). The chain length effect drops out for large 
A' where cy » p'^^x^'^A^ (bold line) as suggested by Eq. <97t . Inset: 
cy /{p^l^x^l^) as a function of the reduced chain length u = N/g{x). 



heat per monomer, cy = ^P^f^aiPfiP)) 11631 . Using the 
second virial of soft BFM beads, Eq. ( [89] ). one obtains 
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as represented by the dashed line. In the large-jc limit this 
yields an exponential decay cy ~ px^ cxp{—x /S) / 16 (dash- 
dotted line) while for x ^ 1 a power-law limiting behavior is 
obtained: cy w j^px^ ^ x^- Eq. ( 196b predicts the energy fluc- 
tuations of BFM beads for essentially all x, slightly underes- 
timating the maximum of cy at x w 10. Chain length effects 
are small for large x where Eq. ( |96b can be used to fit the spe- 
cific heats of polymer melts. Strong A^-effects are, however, 
visible for x ^ 1 where cy increases monotonously with A^. 
This can better be seen from the inset where the specific heat 
is plotted as a function of the reduced chain length u = N/g 
with g being the dimensionless compressibility determined 
in Sec. l4.4m For large chains with m ^ 1 this increase levels 
off at an A'-independent envelope 
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due to the density fluctuation contribution in Eq. ( I69l) . In 
contrast to Eq. ( 195b for the mean energy the density fluctua- 
tion term does now correspond to the leading contribution to 
the numerically measured property. This increases the range 
where the density fluctuation contribution can be demon- 
strated to over three decades in x. Eq. ( 197b is indicated by 
the bold lines in the main panel and the inset of Fig. 1131 

^^ Since e and cy correspond to different derivatives of the free en- 
ergy / with respect to j3, there is obviously no inconsistency in the 
finding that cy reveals larger chain length effects than e. 
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Fig. 14 Dimensionless compressibility gt^{x) = limg^ooS{q,N) vs. x 
for different A' using the same symbols as in Fig.[T3] Main panel: Raw 
data as obtained from the low-wavevector limit of the structure factor. 
Chain length effects become irrelevant for x > 0. 1 if N > 64 and for 
X > 0.001 if A' > 2048. The data are compared to the simple second 
virial approximation l/v2(x)p (dashed line) which reduces to l/(8xp) 
for X ^ 1. As one expects, the compressibility levels off for large x 
and becomes identical to the value g w 0.25, known for the classical 
BFM 1 84 1 (dash-dotted line). Inset: As suggested by Eq. (|98) the ex- 
cess part of the inverse compressibility l/gf^{x) — \/N becomes chain 
length independent, i.e. the data points for all A' collapse. The bold line 
indicates g(x) = limA/_j„„gN(x) from Table[T] 



4.4 Compressibility 

The key control property characterizing the decree of cou- 
pling of the polymer chains due to the imposed penalty e 
is the dimensionless compressibility g{x) = limN-^oogti{x) 
of asymptotically long chains. As suggested by Eq. ( 1431) , 
we compute first the dimensionless compressibility g-^{x) = 
limcj-^0 S{q,N) from the low-^ limit of the total monomer 
structure factor for different overlap penalties x and chain 
lengths A' (see below for details). These raw data are pre- 
sented in Fig. [14] as a function of x. As one expects, g-^{x) 
decreases monotonously with x. Note that the structure fac- 
tor S{q,t) measures the complete compressibility, not just 
its excess contribution. As can be seen, e.g., from Eq. ( 169b 
or from the virial expansion of polymer solutions [3J, the 
compressibility can be written in general as 
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for all X with g^^{x^N) being the excess contribution to the 
c omp ressibility which may, at least in principle, depend on 
A^0As can be seen from the inset of Fig. [141 <2ZZ rescaled 
compressibilities collapse, however, on one A'^-independent 
master curve if one plots l/gex{x,N) = l/g-f^{x) — l/N as 
a function of x, even the compressibilities obtained for un- 
connected beads (A'^ = 1). Within numerical accuracy the A'^- 
dependence observed for g-^ix) can therefore be attributed 
to the trivial osmotic contribution and the excess compress- 
ibility gex "^ I^ is thus identical to the compressibility g{x) 

'^ Small corrections may ari se a s they do arise for the chemical po- 
tential as shown below in Sec.l4.6l 




Fig. 15 Total structure factor S{q) as a function of wavevector q for 
N = 2048 for different overlap penalties x = e/T as indicated. For 
comparison, we have also included the single chain form factor F(q) 
for X = 0.00 1 . The low-wavevector limit of the structure factor is used 
to determine the dimensionless compressibility gNlJ-)- Only for x < 3 
does the structure facto r dec ay monotonously with q as suggested by 
the RPA formula, Eq. ijTOO). G{q) becomes essentially constant for 
smaller temperatures except for wavevectors corresponding to the first 
sharp diffraction peak. The box size L = 256 allows only a direct and 
fair determination of gN (x) for x > 0. 1 . We have been forced to increase 
the box size to L = 512 for smaller x as may be seen for an example 
with X = 0.1 (dash-dotted line). As shown by the bold dashed line, the 
RPA formula is used to improve the estimation of gt^(x) for small x. 



of asymptotically long chains. The bold line indicated in the 
inset presents the best values of ^(x) summarized in Table[T] 
These values have been obtained from the excess compress- 
ibilities for the largest chain length available for .r > 0.001 . A 
precise numerical determination of gex(jc) becomes impossi- 
ble for even smaller overlap penalties. We thus have used for 
the smallest x-values sampled the theoretical prediction 
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for X <^ I due to the postulated free energy, Eq. ( 169] ). The 
prefactor V2(ji:)p representing the bare monomer interaction 
is indicated by the dashed line in the main panel of Fig. [T4] 
Hence, g{x) w l/(8ji:p) = 2/x for weak interactions. The un- 
derlined term is the leading correction due to the density 
fluctuation contribution to the free energy. It implies that the 
excess compressibilities for polymer melts and unconnected 
beads cannot be completely identical. However, the differ- 
ence is far too small to be measurable in the limit where 
Eq. ( |99l l applies. Although this result is unfortunate from 
the theoretical point of view, the data collapse observed in 
the inset suggests that it is acceptable to numerically esti- 
mate the long chain compressibility g{x) by computing the 
structure factors of rather short chains. 

We now turn to the total structure factor 5'(^,A'^) shown in 
Fig. [15] to explain how the compressibilities have been ob- 
tained. Only chains of length A^ = 2048 are presented for 
clarity. Since the wavevectors q used for computing S{q) 
must be commensurate with the cubic simulation box of 
linear dimension L, i.e. the smallest possible wavevector is 



26 





O 


x^O.001 




+ 


x^O.01 




D 


x=0.1 




X 


x=1.0 


nCL 




y=1/(1+Q^) 




.+■+++■++ 



10"'' 10"' 10° 10' lo'^ 10^ 10* 

I I 

10"' 10° 10' 10' 

Q=q^ with^^=bVl2 

Fig. 16 Rescaled total structure factor S{q)/g{x) as a function of the 
reduced wavevector Q = qt, for chain lengtii A' = 8192 and several 
X < 1 as indicated. The screening length t, of the thermal blob is ob- 
tained according to Eq. l l44t using g (x) a nd b{x) from Table [T] The 
bold line compares the data with Eq. nOlt . If replotted as indicated in 
the inset the data collapse on the bisection line. Deviations from the 
RPA formula become visible for larger x as shown for x = 1 (crosses). 



2n/L, it thus is important to have a sufficiently large box 
for a reasonable determination of gn{x)- Note that around 
and above qml monomer structure and lattice effects be- 
come important. Being interested in universal physical be- 
havior we focus on wavevectors g <C 1. For comparison, 
we have also included the single chain form factor F{q) 
for X = 0.001 (bold line). Note that the qualitative shape of 
F{q) — decaying monotonously with q from its maximum 
value F{q = 0,N) = N — depends very little on the overlap 
penalty x (not shown). We remind that the "random phase 
approximation" (RPA) formula EH 
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relates the total structure factor to the measured form fac- 
tor. Eq. ( 1 100b is of course consistent with Eq. ( |98] | in the 
^ ^ limit. It allows to directly fit for the excess compress- 
ibility gex{x,N) fa g{x) using the measured structure factor 
S{q,N) and form factor f(^, A'), at least in the x-range where 
the RPA approximation applies. As may be seen from the 
figure, S{q,N) indeed decreases systematically with x, i.e. 
with decreasing g{x) . For x < 3 it also decays monotonously 
with q, again in agreement with Eq. dlOOb . Interestingly, this 
becomes qualitatively different for larger excluded volume 
interactions {x > 3) where the total structure factor is es- 
sentially constant (in double-logarithmic coordinates), very 
weakly increasing monotonously with q. The RPA formula 
apparently does not apply in this limit in agreement with 
Eq. ( 148b . Fortunately, this is of no concern for our main 
purpose — to compute g{x) — since in precisely this limit 
the compressibility is readily obtained from a broad plateau 
(even for much smaller boxes) which in addition becomes 
chain length independent, as we have already seen from the 
inset of Fiq. (fT4l) . Using boxes with L — 256 it is possi- 
ble to directly measure the plateau values for x < 0.3. For 
smaller jc we have simulated boxes with L = 512 containing 
nmon w 8.4 • 10^ monomers and corresponding to a smallest 



wavevector q w 0.01. This box size becomes again insuffi- 
cient for the smallest reduced overlap penalties x we have 
simulated, as shown in Fig. [l5]for x = 0.001 (dashed line). 



It is for these values where the RPA formula, Eq. dlOOb . al- 
lowing to fit the deviation from the (barely visible) plateau, 
has been particulary useful. 

As already stated in Sec. 12.31 for intermediate wavevec- 
tors (where q corresponds to distances much smaller than the 
radius of gyration and much larger than the monomer size) 
the general RPA Eq. ( 1 1001) may be rewritten as 
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S{q) g{x) 

where we have used that the form factor becomes F{q) k, 
2/(aqY [4|. This assumes that corrections to Gaussian chain 
statistics may be ignored |82,83| and that finite chain size 
effects are negligible. From the numerical point of view the 
approximated RPA Eq. jlOlb has the disadvantage that the 
effective bond length b(x) = \/6a[x) needs to be determined 
first. As shown in Fig. [I6l it has the advantage that it al- 
lows for an additional test of the values of g(x) and b{x) 
indicated in Table [T] The main panel presents the rescaled 
structure factor S{q)/g{x) for chains of length A^ = 8192 as 
a function of Q = ^^ with ^ being obtained from g{x) using 
Eq. ( 144b . All data collapse on the master curve 1/(1 + Q^) 
indicated by the bold line provided that the wavevector q re- 
mains sufficiently small and no local physics is probed. That 
the used compressibilities are accurate is emphasized further 
in the inset where g{x)/S{q) — 1 is plotted as a function of 
Q^ using only sufficiently small wavevectors q. According 
to Eq. dlOlb all data should collapse on the bisection line in 
double-logarithmic coordinates if the correct compressibili- 
ties are used. This is indeed the case. Please note the weak 
deviations visible for x = 1 which is due to the breakdown 
of the RPA formula for large x mentioned above. 

4.5 Chemical potential: Gaussian contribution 

According to Flory's ideality hypothesis the chemical po- 
tential /Jn of polymer melts is expected to be extensive with 
respect to their mass |3|. Fig. 1171 presents the reduced ex- 
cess contribution to the chemical potential, y = lit^/TN, ob- 
tained using thermodynamic integration (as explained be- 
low) for three chain lengths N = I, 64, and 2048 as func- 
tions of x = e/T. As one expects, y{x) increases first linearly 
with X and then levels off. Chain length effects are again 
small on the logarithmic scale chosen in the plot. For large 
X the chemical potential becomes slightly larger for beads 
(y w 2.64) than for long chains where y w 2.1 as shown by 
the dash-dotted line. That the chemical potential of poly- 
mer chains is reduced compared to unconnected beads is 
expected due to the (effectively) attractive bond potential. 
For X ^ 1 this reduction should be described by the density 
fluctuation contribution to the free energy [Eq. ( |69b 1 which 
corresponds to an excess chemical potential 

dp 
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Fig. 17 Reduced chemical potential y = fij^/TN measuring the re- 
versible work for bringing a test chain of length N into a bath of 
chains of equal length A' at volume fraction tj) =0.5. Increasing linearly 
(dashed line) for small x it levels off for large x 2> 1 (dash-dotted line). 
The dashed line shows the second virial approximation y » V2(x)p for 
unconnected beads , fittin g successfully the data below x » 1 . The bold 
line indicates Eq. (T02j. Inset: The chemical potential has been ob- 
tained by thermodynamic integration over the excluded volume inter- 
action of an inserted ghost chain L129J . 
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for X <C 1 with V2 (x) being the second virial of unconnected 
beads. The dashed line in Fig. [T7] presents the leading con- 
tribution V2(a:)p for unconnected beads, the bold line in ad- 
dition the underlined connectivity contribution in Eq. ( 1102b . 
Surprisingly, it turns out that the simple second virial ap- 
proximation provides a better fit over the entire x-range than 
the full prediction. That the density fluctuation contribution 
overestimates the reduction of the chemical potential for x > 
1 is in agreement with Eq. ( 148b . For x^\ where Eq. ( 1102b 
applies in principle the relative correction, scaling as \Jx/ p, 
becomes unfortunately too small to allow a fair test of the 
theory using the measured chemical potential. 

We now explain how the chemical potential presented in 
Fig. [17] has been obtained. Obviously, the simple insertion 
method due to Widom |24| becomes rapidly inefficient with 
increasing x. Generalizing the method suggested in 11291 



13011 we have performed the thermodynamic integration 
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over discrete values of the affinity X = exp(— esgj3/8) char- 
acterizing the excluded volume interaction of a ghost (g) 
chain that is inserted into an equilibrated system (s). (A'sg) 
refers to the mean number of lattice sites where system and 
ghost monomer cube corners overlap at a given interaction 
A. Generalizing the Potts spin mapping, Eq. ( 19 lb , of the 
excluded volume interactions for homopolymers presented 
above, we use now two spin lattices, ^s (r) describing (as be- 
fore) the interaction of the system monomers and So (r ) the 
ghost chain. The spin lattices are kept at the same tempera- 
ture T and are both characterized by the same penalty e = 1 



which has to be paid for a complete overlap of two system 
monomers or two ghost monomers. The interaction of both 
spins is described by 

AE,g = £5s(r)£/,g(5)5g(r+ 5) (104) 

with coupling constants /sg(^) ^ £sg defined as in Eq. ( 1921 ) 
taken apart the energy parameter e which is replaced by the 
tunable interaction energy Egg. Starting with decoupled sys- 
tem and ghost configurations at £sg = 0, i.e. X = \, we in- 
crease the interaction parameter up to Egg = e, i.e. A(e) = 
exp(— ej3/8), always keeping the coupled system at equilib- 
rium. Monitoring the distribution of the number A'sg of over- 
laps between system and ghost cube corners we use multi- 
histogram methods as described in II129II130I to improve the 
precision of the integral. The mean overlap number (A'^sg) is 
shown in the inset of Fig.[l7]as a function of A for A^ = 2048 
and two inverse temperatures x = 3 and x = 100. Starting 
from A = 1 the overlap number decreases monotonously 
with increasing coupling between system and ghost mono- 
mers. Interestingly, a power-law behavior 

{N,,)/N^X'/' (105) 

is found empirically for large x ^ 10 (dashed line). Fitting 
this power law and integrating then analytically Eq. ( 1103b 



provides a useful crosscheck of the numerical integration 
using the multihistogram analysis. This is a technically im- 
portant finding, since the multihistogram analysis requires 
overlapping distributions of N^^ and hence much more equi- 
librated intermediate values A as indicated for x = 100. 



4.6 Chemical potential: Non-extensive corrections 



4.6.1 Theoretical predictions 



We have seen in Sec. 14.51 that the chemical potential jx^ 
of a test chain of length A' plugged into a melt of chemi- 
cally identical monodisperse polymers of the same length 
increases essentially linearly with A^. Focusing now on clas- 
sical BFM systems where £ = oo and setting temperature to 
unity (r = 1), we show how small non-linear deviations may 
be captured numerically |88|. More generally, the challenge 
is to characterize the chemical potential ii„ of a test chain 
of length n immersed into a bath of A^-chains of an arbitrary 
(normalized) length distribution /?n ■ We remind that accord- 
ing to Flory's ideality hypothesis one expects |[3l, 



IJ.n = 11 n for g <^ n <^ (N) , 



(106) 



with jU > being the effective chemical potential /:)er mono- 
mer. The upper boundary (A') indicated in Eq. ( 1106b is due 
to the well-known swelling of extremely large test chains 
where the bath acts as a good solution |3,7|. As we have 
pointed out in Sec. 11.51 Flory's hypothesis has been chal- 
lenged by the discovery of long-range correlations imposed 
by the incompressibility constraint. These correlations lead 
to the systematic swelling of chain segments as further dis- 
cussed in Sec. \5\ II80II86II . We question here the validity of 



28 




Fig. 18 Non-extensive deviation of the cJiemical potential Sun/c^ as a 
function of tlie test ciiain length n as predicted by perturbation calcula- 
tion 1 88 1. The reference for the chemical potential /!„ is set here by the 
chemical potential of Gaussian chains of same effective bond length 
b. Flory-distributed (FD) melts are indicated by solid lines, monodis- 
perse (MD) melts by dash-dotted lines (Fade approximation). The bold 
lines refer to test chains in melts of constant (mean) chain length with 
(N) = 100. The deviation changes sign at n Ri [N). The thinner lines 
refer to test chains of same length as the typical melt chain, n = (N) , 
where the asymptotic Gaussian behavior (5^„ = 0) is approached sys- 
tematically from below. 



Flory's hypothesis for an important thermodynamic prop- 
erty, the chemical potential /!„ of a test chain inserted into 
a three-dimensional melt. One expects that the correlation 
hole potential (Fig. O leads to a deviation 



5jUn = Mn - M« ~ u*{n) ^ +l/\/n for n < (A^) 



(107) 



that is non-extensive in chain length and this irrespective of 
the distribution p^ of the bath. Assuming a quenched Flory- 
size distribution, Eq. ( 1201) . it can be demonstrated as shown 
in Appendix IB .61 that 



5Mn « -^ (1 - 2iin) for g < n < (TV)^ 



(108) 



where we have set c^ = 3cs/8 using the swelling coefficient 

Cs = y^TA/n^/pb^ defined in the Sec.[T3] Eq. ( flOSl l is rep- 
resented by the bold solid line in Fig. [18] (A corresponding 
prediction for a monodisperse bath [88] is indicated by the 
dash-dotted line.) As anticipated by Eq. ( 11071) . the first term 
in Eq. (1108b dominates for short test chains. The second term 
dominates for large test chains with n > l/(2/j) becoming 
non-perturbative for n^ 1/m^- Both contributions to 5jin 
decrease with increasing «o Interestingly, while 5/in de- 
creases at fixed (A'^), it increases as SfJ-n/c^ = —\l\fn (thin 
solid line) for a test chain with n = (A'^). The chemical po- 
tential of typical chains of the bath approaches thus the ideal 
chain limit from below. 

Flory-distributed polymer melts are obtained naturally in 
systems of self-assembled linear EP where branching and 

^'' This corresponds to an effective enhancement factor of the parti- 
tion function quite similar to the 5;Un = — (y — 1 ) log («) in the standard 



excluded volume statistics with y f 
walk susceptibility exponent (S). 



1.16 > 1 being the self-avoiding 



the formation of closed rings are forbidden \\1\ . Since the 
suggested correction, Eq. ( 1108b . to the ideal chain chemical 
potential is weak, the system must remain to leading order 
Flory distributed and Eq. (1108b should thus hold|3 This cor- 
rection implies for the annealed length distribution of linear 
EP that (to leading order) 
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where both the lower (g <^ N) and the upper limit (A'^ <C 
(A') ) of validity are irrelevant in the large- (A') limit. Note 
that Eq. dllOb is properly normalized, i.e. the prefactor /i of 
the distribution remains exact if 5/iN is given by Eq. ( 1108b . 
At given jU the first moment increases slightly 



and, more generally, one expects for the pth moment 
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p\ p\ 

with r{x) being the Gamma function 1441 . The non-expo- 
nentiality parameter Lp = 1 — (A''') /p\ {N)P thus scales as 

Lp =WpC^^ with 

w„ = {r{p+\/2) + ^pp\-2r{p + ?>/2))/p\ (113) 

being a /j-dependent geometrical factor0 

4.6.2 Computational results for EP melts 

Eq. (IllOb and Eq. ( 11131) allow us to demonstrate numeri- 



cally the prediction, Eq. ( 1108b . from the observed non-ex- 
ponentiality of the length distribution of EP melts obtained 
as described in Sec. 13.51 As demonstrated in the main panel 
of Fig. [To] these EP melts are indeed essentially Flory dis- 
tributed, i.e. the length distribution p-^ decays to leading 
order exponentially with the reduced chain length N / (N) . 
However, deviations for N/ (N) <C 1 are visible in Fig. [TOl c) 
in qualitative agreement with the predicted positive devi- 



ation of the chemical potential, Eq. ( 1107b . The curvature 
of — log(/?N) is further analyzed in Fig. [19] Motivated by 
Eq. (1109b . we present in panel (a) the functional 



V[p^] = -log{pj^)-ixN+ login) 



(114) 



^' The chemical potential of a test chain does depend on the length 
distribution pN of the bath. However, for an infinite macroscopically 
homogeneous systems it is independent on whether this distribution is 
annealed or quenched, i.e. if it is allowed to fluctuate or not. This fol- 
lows from the well-known behavior of fluctuations of extensive param- 
eters in macroscopic systems: the relative fluctuations vanish as l/\/V 
as the total volume V — > oo. The latter limit is taken first in our calcu- 
lations, i.e. we consider an infinite number of (annealed or quenched) 
chains. The large-A' limit is then taken afterwards. 

^^ The polydispersity index /, i.e. the ratio of weight average and 
number average, becomes / = {N^) / (N) = 2(1 — W2Cf, ,fji) < 2. In 
this sense, the distribution becom es narrow er just as for dilute good 
solvent EP where / = 1 -I- l/y < 2 I107I164I . 
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Fig. 19 Characterization of pN using the functionals (a) V[pn] ~ Sjit^ 
and (b) Vt^fpisj] = 2V[pn] — V[p2n] which should both vanish for per- 
fectly Flory-distributed ideal chains as indicated by the horizontal lines 
in the panels. Data points for different scission energies E collapse if 
VfpNl/c^j-yT' and W[pN]/cfi\/]u(2— l/\/2) are plotted vs. the reduced 
chain length x = fiN. For small x both functionals decay as l/\/x as 
shown by the d ash-d otted lin es. Th e bold lines correspond to the full 
predictions Eq. (TT6) and Eq. HlSt for V[pn] and W[pn]i respectively. 



where the second term takes off the ideal contribution to the 
chemical potential. The last term is due to the normalization 
of pn and eliminates a trivial vertical shift depending on the 
scission energy E. Consistently with Eq. (Illlb . the chemical 
potential per monomer fi has been obtained from the mea- 
sured mean chain length (A'^) using 



M = WMi + c 
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Note that /i and 1/ (N) become numerically indistinguish- 
able for ii > 7. If the Gaussian contribution to the chemical 
potential is properly subtracted one expects to obtain directly 
the non-Gaussian deviation to the chemical potential, 5/iN ~ 
y [/)n]- Due to Eq. (1108b the functional should thus scale as 

V[p^]/c^^^{l~2x)/V^ (116) 

withx = fiN as indicated by the bold line in the panel. This is 
well born out by the data collapse obtained up to x w 5. Ob- 
viously, the statistics deteriorates for x ^ I for all energies 
due to the exponential cutoff of /?n. Unfortunately, the statis- 
tics of the length histograms decreases strongly with E and 
becomes too low for a meaningful comparison for i? > 9. It 
is for this numerical reason that we use Eq. (1115b rather than 
the large -ii limit ;U = 1 / (A'^) since this allows us to add the 
two histograms for E = 5 and E = 6 for which high preci- 
sion data is available. Otherwise these energies would devi- 
ate from Eq. (1116b for large x due to an insufficient subtrac- 
tion of the leading Gaussian contribution. 

Since the subtraction of the large linear Gaussian contri- 
bution is in any case a delicate issue we present in panel (b) 
of Fig.[T9]a second functional. 
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Fig. 20 Non-exponentiality parameter Lp for different moments p as 
a function of mean chain le ngth ( TV) . Lp is finite decreasing with chain 
length as suggested by Eq. <113t . The vertical axis has been rescaled 
successfully using the p-dependent weights Wp indicated in the inset. 



where by construction this contribution is eliminated follow- 
ing a suggestion made recently by Semenov and Johner [94 1. 
The normalization factor /J appearing in Eq. ( 11171) elimi- 
nates again a weak vertical scission energy dependence of 
the data. Obviously, W [pf^] = for perfectly Flory-distributed 
chains. Following Eq. dlOSI) and Eq. ( 1109b one expects 

^[/^n] 1-0^06. 

with X = ixN. Eq. (1118b is indicated by the bold line which 
compares again rather well with the presented data. 

The functionals presented in Fig. [19] require high-accur- 
acy histograms. That /?n is only approximately Flory dis- 
tributed can be more readily seen using the "non-exponential- 
ity parameter" Lp = I - (NP) /pi {N)'\ Obviousl y, L,, = 
for rigorously Gaussian chains. As stated in Eq. (1113b . we 
expect the non-exponentiality parameter to decay as Lp = 
WpCfiy/Ji ^ \/ \J{N), i.e. as the correlation hole potential of 
the typical melt chain. The main panel of Fig. [20] presents 
Lp/wpC^ as a function of (A') w l//i. The predicted power- 
law decay is clearly demonstrated by the data. Note that the 
scaling of the vertical axis with the p-dependent geometrical 
factors Wp allows to bring all moments on the same master 
curve. As can be seen from the inset of Fig. [20] this scaling is 
significant since Wp varies over nearly a decade. Deviations 
from the predicted scaling are visible, not surprisingly, for 
small {N) < 10. Since the coefficient c^ is known the data 
collapse on the theoretical prediction (bold line) has been 
produced without any free adjustable parameter. 



4.7 Summary 

In this section we have discussed various thermodynamic 
properties of BFM melts at our reference volume fraction 
^ = 8p = 0.5. For large overlap penalties x = ejT ~^ 1 all 
thermodynamic properties remain similar to the known val- 
ues for athermal BFM systems (x = oo) ||I29III30ll84l . From 
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the computational point of view this is of some interest since 
it suggests the equilibration and sampling of configurations 
at a reasonable cost at overlap penalties between x = I and 
X = 10 where the slithering snake algorithm is still efficient. 
Decreasing the monomer interaction further we have found, 
as one expects, good agreement with the prediction by Ed- 
wards, Eq. ( |69l l. for small overlap penalties (x -^ 1). In this 
limit the thermodynamics is essentially determined by the 
second virial contribution, Eq. ( l89l) . Interestingly, being the 
second derivative of the free energy with respect to the in- 
verse temperature the specific heat c\ allows the verifica- 
tion of the predicted scaling of the free energy contribu- 
tion —kBT/^^p due to the chain connectivity (Fig. [13]). The 
central control parameter of this work, the dimensionless 
compressibility g = limAr_j.oo(lim^_j.o5'(^,A'^)), is discussed 
in Sec. 14.41 While the determination of g{x) is numerically 
trivial for large overlap penalties, more care is needed for 
X < 1 (Fig. [ni. Note that the static RPA, Eq. (051), breaks 
down in the opposite limit {x ^ 1) where the structure fac- 
tors of beads and polymer melts become ultimately identical. 

For a test chain of length n = N immersed into a bath 
of monodisperse chains of length A^ we have confirmed by 
thermodynamic integration in Sec. 14.5 I that the chemical po- 
tential /In is essentially extensive with respect to the chain 
length, as expected from Flory's ideality hypothesis. In order 
to demonstrate that the small deviations predicted by the per- 
turbation theory |88| exist, we have investigated in Sec. 14.61 
essentially Flory-distributed (Fig. 1101) self-assembled EP sys- 
tems. The detailed analysis of the measured length distribu- 
tion pn (Fig. [19]) has allowed us to show the existence of 
a correction dlJ.„, Eq. dlOSI) , scaling essentially as the cor- 
relation hole penalty 5n„ w u*{n) for n <^ (N). The ob- 
served deviations for EP systems II88I beg for an improved 
numerical verification by means of the thermodynamic in- 
tegration method (Sec. 14.51) for test chains immersed into 
monodisperse melts 



5 Intramolecular conformational properties 

5.1 Introduction 

We turn now to the description of intrachain conformational 
properties of BFM melts at volume fraction (/) = 0.5 and tem- 
perature T = I comparing our numerical results to the the- 
oretical predictions announced in Sec. 11.51 and Sec. |2] We 
discuss first the properties on monomeric level (Sec. 15.2b 
and show then how the effective bond length b for asymp- 
totically long chains can be obtained using the theoretical 
input developed in Sec. [2] We present then the bond-bond 

'^ The presented results for dense polymer solutions may also be of 
relevance to the chemical potential of dilute polymer chains at and 
around the ©-point which has received attention recently |74 165 1. The 
reason for this connection is that (taken apart different prefactors) the 
same effective interaction potential \{q) ~ q^ enters the perturbation 
calculation in the low wavevector limit. A non-extensive coixection 
5/Xn ~ +\/-/n in three dimensions is thus to be expected. 




Fig. 21 Various bond propeities vs. overlap penalty e for volume frac- 
tion ^ = 0.5 and temperature T = \. The data for the root-mean-square 
bond length /(e) and the effective bond length b{e) = limA/^^&N(e) 
for asymptotically long chains are given in Table [T] The da sh-do tted 
line indicates the effective bond length as predicted by Eq. l ll20t as- 
suming bo = I{e) for the bond length of the referenc e cha in. The bold 
line shows the fix points obtained by iteration of Eq. U20t . 



correlation function as a function of curvilinear distance s 
(Sec. 15.4.1b and as a function of the spatial distance r be- 
tween both bonds (Sec. l5.4!2] l. Higher moments of the dis- 
tribution G{r,s) and its deviations dG{r,s) from Flory's ide- 
ality hypothesis will be analyzed in Sec. 15.51 and Sec. 15.61 
respectively. We turn finally to the characterization of an ex- 
perimental relevant observable, the intramolecular form fac- 
tor F(^). 



5.2 Bond properties 

We begin by characterizing local-scale features of the al- 
gorithm described in Sec. [3] By definition of our version 
of the BFM algorithm the bond length is allowed to fluc- 
tuate between 2 and \/T0. One expects that switching on the 
overlap penalty e will suppress large bonds due to the in- 
creasing pressure. The mean bond length is characterized by 
the root-mean- square length /. (Other moments yield similar 
results.) The mean bond length rapidly becomes (A^ > 20) 
chain length independent |50|. As can be seen from Fig. 1211 
I shows a monotonous decay between e w 3 and e w 20. 
As other local properties, the bond length becomes constant 
in the small-e and large-e limits (dashed lines). The value 
l{e = 0) = 2.718 gives the lower bound for the effective 
bond length b{e) of asymptotically long chains (stars) ob- 
tained below. 

Defining the bond angle 6 between two subsequent bonds 
by the scalar product cos(0) = e„ ■ e„+i of the normalized 
bond vectors e, = /,/ 1 /, | , the local chain rigidity may be char- 
acterized by (0) and (cos(0)). Note that (0) and (cos(0)) 
can be regarded as chain length independent, just as the mean 
bond length. The local rigidity is negligible for £ ^ 1, i.e. 
(0) w 90° and (cos(0))) w due to the symmetry of the 
distribution p{d) with respect to 90°. The rigidity then in- 
creases around e w 1 and becomes constant again for large e 
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where (0) w 82.2° and (cos(0)) w 0.106 |86|. The increase 
of the local rigidity for larger excluded volume interactions 
is of course expected due to the suppression of immediate 
backfoldings corresponding to bond angles 6 > 143° | f52]| . 
The distribution p{d) therefore becomes lopsided towards 
smaller 6 (not shown). According to Eq. dlSb the effective 
bond length of a "freely rotating" (FR) chain is given by 
b{e) = Z(£)^?ii with CFR = (1 + (cos(0)))/(l - (cos(0))). 
This simple model, indicated by the crosses in Fig.|2T] yields 
a qualitatively reasonable trend (monotonous increase of the 
effective bond length at e w 1) but fails to fit the directly 
measured effective bond lengths quantitatively [,86 J . 



5.3 Mean-squared chain and subchain size 
5.3.1 Total chain size /?n 
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1 for different e. The chains only remain Gaus- 
sian on all scales and all A' for extremely small e. For e > 0. 1 one ob- 



One way to characterize the total chain size R^ is to measure .^^^^^ ^^(^^ j^ ^^^^y li^^^^ly -^ agreement with Eq. (TTD- This can be 



the second moment of the chain end-to-end distance R 



N 



((''A' -ri)^). We consider the ratio feN(e) =/?N(e)/V^- 1 
to compare the measured chain size with the ideal chain be- 
havior which is commonly taken as granted |3,4,39| and 
which is the basis of our perturbation calculation (Sec. 12.41) . 
The task is to extrapolate for the effective bond length /?(£) = 
limAr_j.oo ^N (s) of asymptotically long chains. The ratio /?n (s) 
for N = 64 and N = 2048 and the asymptotic limit b{£) — 
obtained by extrapolation as described below — are pre- 
sented in Fig. [211 Obviously, b]^{£) -^ l{£ = 0) for all N 
in the small-£ limit. Z?n(£) increases then in the intermediate 
e-window before it levels off at e w 10. The swelling due 
to the excluded volume interaction is the stronger the larger 
the chain length, i.e. bj<s{e) increases monotonously with A^. 
This swelling thus cannot be attributed to a local rigidity as 
described, e.g., by the freely-rotating chain model. 

The A'-effect can be seen better in Fig. |22] where we 
have plotted b^ for several penalties e as a function of f = 
1 l^/N — 1. The choice of the horizontal axis is motivated by 
Eq. jlOb suggesting the linear relation 



bl{e)Kb^{e) (l-c(e)cs(e)0 for7V/^>l 



(119) 



with Cs = ^lAjii^ I pb^ being the swelling coefficient de- 
fined in Sec. ll.Sl and c(e) an additional numerical prefactor 
of order unity. This prefactor has been introduced in agree- 
ment with Eq. (|86] |. The reason for this coefficient is that 
the corrections to Gaussian behavior differ slightly for inter- 
nal chain segments [as described by Eq. (llOb I and the total 
chain size which is characterized in Fig.|22l We remind that 
c -^ I{°°) ~ 1 .59 for N ^°°. However, since this value corre- 
sponds to the limit of a very slowly converging integral [i84i 
it is better to use Eq. ( |119l l as a two-parameter fit for b{e) and 
c(e) and to crosscheck then whether the fitted c is of order 
unity. As shown in the figure for three overlap penalties, this 
method can be used reasonably for overlap penalties as low 
as e w 0.1, albeit with decreasing £ it systematically under- 
estimates the "true" ^(e)-values indicated in Table[T] Please 
note that N/g w 400 for e = 0. 1 and A^ = 8192. Chains with 



used for a simple two-parameter fit for the effective bond length b{e) 
as indicated for £ = 0. 1, 1.0 and oo. 



N ^ 8192 would be required to use this method for even 
smaller e. In this limit it is better to use as a first step the 
value bj^{£) of the largest chain length simulated as a (rather 
reasonable) lower bound for b{e). 



5.3.2 Subchain size Rs 

The mean-squared subchain size R^ is presented in Figs. [23] 
and [24] where we focus on melts without monomer over- 
lap (e = oo) [841 . We show how the effective bond length 
b for asymptotically long chains may be extrapolated from 
the measured R^ for finite arc-length 5 and finite chain length 
N using the perturbation prediction Eq. (llOb . 

As already remarked by Auhl et al. Ii73l . Fig. [23] shows 
clearly that the chains are swollen, i.e. R'^/s increases sys- 
tematically and this up to very large arc-length s. In agree- 
ment with Eq. (llOb , the Gaussian behavior (horizontal line) 
is approached from below and the deviation decays as m* ^^ 
l/\/s. The indicated bold line corresponds to ^ = 3.244 and 
Cs w 0.41 which fits nicely the data over several decades in .s 
— provided that chain end effects can be neglected {s <^ N). 
Note that a systematic underestimation of the true effec- 
tive bond length would be obtained by taking the largest 
R^/s K, 3.23^ value available, say, for monodisperse chains 
of length A^ = 2048. Interestingly, R^/s does not approach 
the asymptotic limit monotonically 1173 11841 . Especially for 
short chains one finds a non-monotonic behavior for i — > 
A'^. This means that the total chain end-to-end distance 7?n 
shows even more pronounced deviations from the asymp- 
totic limit. We emphasize that the non-monotonicity of R^/s 
becomes weaker with increasing A^ and that, as one expects, 
the inner distances, as well as the total chain size, are char- 
acterized by the same effective bond length b for large .s or 
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Fig. 23 Mean- squared subchain size R^/s vs. s for different chain 
length A'. Log-linear coordinates are used to emphasize the power law 
swelling over several orders of magnitude of s. The data approach 
the asymptotic limit (horizontal line) from below, i.e. the chains are 
swollen. This behavior is well fitted by Eq. i fTot ior \ <^s ^N (bold 
line). Non-monotonous behavior is found for s — *■ A', especially for 
small N |73|. The dashed line indicates the measured total chain end- 
to-end distances b^ showing even more pronounced deviations. The 
dash-dotted line compares this data with Eq. (86). 



A^P^ The non-monotonic behavior may be qualitatively un- 
derstood by the reduced self-interactions at the chain ends 
which lessens the swelling on these scales. Our perturba- 
tion prediction Eq. ( 186b is indicated in Fig.|23]by the dash- 
dotted line. As already remarked, this prediction contains 
a slowly converging numerical integral of order unity. The 
dash-dotted line corresponds to the value c = /(oo) w 1.59 
predicted for asymptotically long chains, the dashed line to 
the slightly smaller value c = 1 .45 fitted in Fig. |22] In sum- 
mary, it is clear that one should use the subchain size R^ 
rather than the total chain size 7?n to obtain in a computa- 
tional study a reliable fit of the effective bond length b. 

The representation chosen in Fig.|23]is not the most con- 
venient one for an accurate determination of b and c^. How 
precise coefficients may be obtained according to Eq. dlOl) is 
addressed in the Fig.|24]for chains of length A^ = 2048. In a 
first step R'^/s should be plotted as a function of l/y^ just 
as in Fig. |22] This allows for a first rough estimation of b. 
Since data for large .s are less visible in this representation, 
we recommend for the fine-tuning of b to switch then to log- 
arithmic coordinates with a vertical axis Ki (s) = I — R^/b^s 
for different trial values of b. The correct value of b is found 
by adjusting the axes such that Ki {s) extrapolates linearly as 
a function of Cs/^/s to zero for large .sO We assume for the 
fine-tuning that higher order perturbation corrections may be 
neglected, i.e. we take Eq. jlOb literally. The plot shows that 
this method is very sensitive, yielding a best value b = 3 .244 
for e = oo that agrees with the theory over more than one 

'* We note en passant that this is not the case for the compact chain 
conformations adopted by chains in strictly 2D melts where the effec- 
tive bond length associated to the chain ends is smaller than the one for 
subchains even for A? ^- oo 1 120|. 

^^ For the bead-spring model presented in Ref. 1 84] the empirical 
swelling coefficient Cj is shown to differ for yet unknown reasons by a 
few percents from the predicted value \J 2\ j ■n} j pb^ . 
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Fig. 24 Replot of the mean-squared subchain size as K\ (s) = 1 — 

Rl/b^s vs. c^/y/s with Cs = \JlApi? j pb^ , as suggested by Eq. i fTot . 
for different trial effective bond lengths b as indicated. Only chains of 
length N = 2048 and e = oo are considered for clarity. This procedure is 
very sensitive to the value chosen and allows a precise determination. 



order of magnitude without curvature. Using this method 
we have determined the effective bond lengths b{£) for a 
broad range of overlap penalties e as indicated in Table [T] 
and shown by the stars in Fig. l2T]|.86J . 



5.3.3 Predicting the effective bond length 

Up to now, we have used the theoretical results ioft the ef- 
fective bond length b(z), rather than to predict it from the 
known thermodynamic properties and local model features 
such as the bond length Z(e). As reminded in Eq. d72] l the 
increase of the effective bond length for weakly interacting 
and asymptotically long polymer melts has been in fact cal- 
culated long ago by Edwards |J4l. Following a suggestion 
made by Muthukumar and Edwards pS | this perturbation 
result may be rewritten as a recursion relation 
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Gi with Gi ; 



1 



^/i^P 



(120) 



with bi being the bond length after the ;-th iteration step, 
g the measured dimensionless compressibility and Gi{bj,g) 
the relevant Ginzburg parameter quantifying the strength of 
the interaction acting on a chain segment of length s = g. 
Since Gz » Gi becomes small for large compressibilities 
g{s), one expects good agreement with our data for small e. 
If we set bo = /(e) this yields after one iteration {i = 1) the 
dash-dotted line indicated in Fig. [21] i.e. a reasonable predic- 
tion is only achieved up to £ w 0.01. The predictive power 
of Eq. ( 11201) can be considerably improved over nearly two 
decades up to e w 1 if one applies the formula iteratively 
using the effective bond length b, obtained at step / as input 
for the Ginzburg parameter for computing /?,+i. This recur- 
sion converges rapidly as shown by the bold line indicated 
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in Fig.EUobtained a fter 2 iterations Note that Gz < 0.34 
for e < 1 where Eq. ( 1120b fits our data nicely. The fix-point 
solution of Eq. (1120b does not capture correctly the leveling 
off of ^(e) setting in above e w 1. Since the Ginzburg pa- 
rameter becomes there of order one, this is to be expected. 
In summary, we have shown that the iteration of Eq. ( 1120b 
allows a good prediction for b{£) for Gz(e) <C 1. If reliable 
values for compressibilities g{e) ^ 1 are available, this is 
the method of choice if one cannot afford to simulate very 
long chains. 



5.4 Bond-bond correlation function 

5.4.1 Bond-bond correlation function Pi {s) 

As we have seen in Fig. |24l to demonstrate the deviations 
from Flory's ideality hypothesis starting from the subchain 
size 7?s requires to subtract a large Gaussian contribution 
b^s. Unfortunately, this requires as a first step the precise 
determination of the effective bond length b{e) for asymp- 
totically long chains which might not be available. Indeed 
we have used in the preceding Sec. 15.3.21 the fact that the 
scaling of Ki (s) critically depends on this accurate value 
to improve the estimation of b{£). Hence, it would be nice 
to demonstrate directly the scaling implied by our key pre- 
diction without any tunable parameter. The trick to achieve 
this is similar to our demonstration of the density fluctu- 
ation contributions to the free energy, Eq. ( 1691 ), presented 
in Sec. 14.31 We consider the curvature of R^, i.e. its sec- 
ond derivative with respect to s. Using Eq. (|4| this second 
derivative is obtained directly from the bond-bond correla- 
tion function Pi [s) computed by averaging over all pairs of 
monomers («, m = n + s). We remind that P\ {s) is generally 
believed to decrease exponentially as in Eq. ( 115b . This text- 
book belief is based on the assumption that all long range in- 
teractions are negligible on distances larger than ^ . Hence, 
only correlations along the backbone of the chains should 
matter and it is then straightforward to work out that an ex- 
ponential cutoff is inevitable due to the multiplicative loss of 
any information transferred recursively along the chain IJl. 
The bond-bond correlation function P\{s) obtained for 
monodisperse chains at (/) = 0.5 and £ = oo is presented in 
Fig- ESI The power-law decay with exponent ft) = 3/2 pre- 
dicted by our key perturbation result Eq. ^ is perfectly 
confirmed by the larger chains (A'^ > 256). As can be seen 
for N — 16, exponentials are compatible with the data of 
short chains, however. This might explain why the power- 
law scaling has been overlooked in older numerical stud- 
ies, since good statistics for large chains (A' > 1000) has 
only become available recently. However, it is clearly shown 
that Pi (s) approaches systematically the scale-free asymp- 
tote with increasing A^. The departure from this limit is fully 

^^ Essentially the same result is obtained up to £ Ri 1 if one sets di- 
rectly bQ = b using the measured effective bond length as start value of 
the iteration. 




Fig. 25 The bond-bond correlation function Pi{s)/cp as a fun ction o f 
the curvilinear distance s for a broad range of chain lengths N f80'84). 
Provided that 1 ^ ^ ^ A', all data sets collapse on the power law slope 
with exponent O) = 3/2 (bold line) as predicted by Eq. ([9}. The dash- 
dotted curve P\ {s) » exp(— ^/1.5) shows that exponential behavior is 
only compatible with very small chain lengths. The dashed lines cor- 
respond to the theoretical prediction, Eq. | |12U . for short chains with 
N=16, 32, 64 and 128 (from left to right). 



Eq.(83) 




u=s/g 

Fig. 26 Pi {s) for different overlap penalties e f861. Inset: Pi (.s) in log- 
log coordinates. The data approaches a power law behavior, Pi (s) ~ 
l/.s'", with exponent ft) = 1/2 for small e (dashed line) and o) = 3/2 
for e > 1 (dash-dotted line). Main panel: Pi{s)/ [cp{/g^/~] vs. u = 
s/g as suggested by Eq. ([83}. For large u, where an incompressible 
packing of thermal blobs is probed, all data collapse onto the dash- 
dotted line as predicted by Eq. {9}, i.e. Pi (s) becomes independent of 
the compressibility g{e). 



accounted for by the theory if chain end effects are care- 
fully considered (dashed lines). Generalizing Eq. (|9]l using 
the Pade approximation, Eq. ( 1631) , perturbation theory yields 



Pi{s) = 



cp 1+3m + 5m 



„3/2 



1 



;i-«r 



(121) 



where we have set u = \Js/N |f84l . For m ^ 1 this is con- 
sistent with Eq. (|9]l. In the limit of large s ^ N, the corre- 
lation functions vanish rigorously as Pi{s) °^ (1 — m)^. Con- 
sidering that non-universal features cannot be neglected for 
short chain properties and that the theory does not allow for 
any free fitting parameter, the agreement found in Fig.l25lis 
rather satisfactory. 
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The bond-bond correlation function Pi{s) for different 
overlap penalties e is presented in Fig.|26]for chains of length 
A'^ = 2048 1 86 1. As can be seen from the unsealed data shown 
in the inset, Pi (s) approaches a power law with exponent 
(0=1/2 (dashed line) in the limit of weak overlap penal- 
ties in agreement with Eq. ( I84I) . For £ > 1 our data is com- 
patible with an exponent ft) = 3/2 (dash-dotted line) as sug- 
gested by Eq. ^. Hence, we have demonstrated without any 
tunable parameter that Flory's ideality hypothesis is system- 
atically violated for all segment lengths s and all overlap 
penalties e. As suggested by Eq. ( 183b . the main panel of 
Fig.l26lpresents Pi {s)/{cp/g^i^) as a function of the reduced 
arc-length u = s/g using the dimensionless compressibilities 
g{e) and effective bond lengths ^(e) from Table[T] The data 
collapse is remarkable as long d& \ <^s ^N . The relation 
Eq. ( 183b is indicated by the bold line; it is in perfect agree- 
ment with the simulation data|3 The asymptotical behavior 
with ft) = 1/2 for M ^ 1 and ft) = 3/2 ioru^ \ is shown by 
the dashed and dash-dotted lines, respectively. As predicted 
by Eq. (|9]l, one recovers the power law Pi [s) = cp/s^'^ irre- 
spective of the blob size g. This suggests that the exponent 
ft) = 3/2 is not due to local physics on the monomer scale, 
since for i ^ ^ ^ 1 distances much larger than the monomer 
or even the thermal blob are probed. 

5.4.2 Distance dependence of angular correlations 
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Fig. 27 The angular correlations are characterized by the bond-bond 
correlation function Piir) = {l„-l_,„)/l^ averaging over all pairs of 
bonds of a chain of same distance r = \r\. The data shown has been 
obtained for EP of a broad range of scission energies E. If Pi (r) (A') is 
traced as a function of r/ r^ al l data points collapse. Confirming the 
general scaling idea Eq. (I122t the bold line indicates the predicted 
power-law asymptote Eq. ([123) for r/r, ^ 1, i.e. the bond-bond corre- 
lations are truely long-ranged. The plateau P\ (r) (A') = c* = 0. 14 seen 
for r/r^ 2> 1 is due the the enhanced weight of stretc hed c hain seg- 
ments (87). The complete perturbation prediction Eq. <125t given by 
the thin line interpolates perfectly between the power-law asymptote 
(bold line) and the plateau (dashed line). 



We have seen in the previous paragraph that Pi (s) decays for 
s/g ^ 1 as a power law with an exponent (0 = dv = 3/2. 
Since the decay of Pi(.s) resembles the return probability of 
a random walk Go{r -^ 0,s), it is tempting ||74| to attribute 
the observed effect to "local self-kicks" at a distance r Ki a 
involving the bonds Z„ and /,„ themselves (or their immediate 
neighbors). Accordingly, the bond-bond correlation function 
should reveal a 5 (r) -correlation if sampled as a function of 
the distance r = \r\ between bond pairs. This interpretation 
turns out to be incorrect, however, and we will show that the 
power law in i simply translates as 



i-"o 



r/b) 



-m/v 



(r/b)- 



(122) 



As demonstrated analytically in Sec. IB. 51 one expects indeed 
for incompressible solutions of infinite chains that 



Pi(r) » P„{r) = °° 3 for ^ < r < r. 



(123) 



I2npr^ 

as suggested by Eq. ( 11221) . i.e. the angular correlations are 
genuinely long-ranged. As discussed in Sec. IB .5 .41 the up- 
per cutoff r* arises due to the enhanced weight of stretched 
chain segments which align bond pairs for distances r ^ r, . 
Unfortunately, this cutoff increases rather slowly with chain 
length ill 

r*«^(7V)^/''</?N«^(A^)V (124) 

The simulation of computationally challenging chain lengths 
thus is required to demonstrate numerically the predicted 

^' Eq. l l83t has also been applied success fully to the Pi{s) obtained 
from single-chain-in-mean-field simulation 11661. 



power-law decay of Pi (r). Generalizing Eq. ( 1123b for Flory- 
distributed EP one obtains according to Eq. ( 12201) that 



Pi{r)=Po.ir)h{x)+c^lJ. 



(125) 



with h{x) = (1 + 2x)'^ exp{—2x) being a scaling function of 
X = y^r/2a and c* a phenomenological constant set (in 
practice) by the finite local chain rigidity. The angular cor- 
relation function Pi (r) for EP systems is thus predicted to 
level off for large reduced distances M 

We present in Fig. [27] numerical results obtained for sys- 
tems containing Flory-distributed EP. As indicated in the 
sketch. Pi (r) is obtained by averaging of all intrachain bond 
pairs of same distance r. To avoid trivial correlations the sec- 
ond bond /,„ is outside the subchain between the monomer 
n and m defining the vector r [87 89|. Model-depending 
physics not taken into account by theory obviously becomes 
relevant for short distances corresponding to subchains of 
a couple of monomers. For clarity, we have thus omitted 
data points with r < 5. For small distance r ^r^, the (un- 
sealed) bond-bond correlation function Pi (r) is found to de- 
cay strongly as expected from Eq. ( 1125b . Note that in this 
limit Pi (r) does neither depend on the mean chain length 
nor the chain length distribution, i.e. the same behavior is 
observed for monodisperse chains. See Ref. |87| for details. 
To scale away the (A'^) -dependence in the large-r limit, we 
trace in Fig. [27] the rescaled bond-bond correlation function 

*^ For monodisperse systems one finds instead P\ (r) ^ {r/N)~ for 
very large distances r with R^ <S^ r < Nl due to the increasing contri- 
bution of stretched bonds [87]. 



35 



lo" 


- 




^ 


j^^ too small s ! 


2 10-' 


; / 




f 


O P="" 
D P=2 
O P=3 
A p=4 
< P=5 


10" 


\J^ f= 


=0.5, 


N 


=2048, e=~, b=3.244 



10^ 



, 1/2 

X ~c /s 

s 

Fig. 28 Test of Eq. <77t where the rescaled moments -S'p(i) are plotted 
vs. X ~ Cs/\fs. Setting b = 3.244 for all moments, all data sets extrap- 
olate linearly to zero for small x. The saturation at large x is due to 
the finite extensibility of subchains. Since this effect becomes more 
marked for larger moments, the fit of b is best performed for p = I. 



Pi (r) (A'^) as a function of r/r^, with r^=b (N) ' using the 
mean chain lengths (A'^) indicated in Table [3j Note that the 
error bars (not shown) become clearly much larger than the 
symbol size for large bond energies E > 12. It is fair to state, 
however, that all data points collapse nicely on the one mas- 
ter curve indicated by the thin line predicted by Eq. ( 11251) . 
The asymptotic power law Eq. i 



11231) is indicated by the bold 
line. That Pi (r) for EP becomes constant for r/r* ^ 1 con- 
firms a non-trivial prediction of the theorvr^ 



5.5 Higher moments 

The preceding discussion focused on the second moment of 
the subchain size distribution G{r,s) and its derivatives with 
respect to s. We have also computed higher moments {r^'') 
with p <5 ||84]| . (For clarity we focus here on monodisperse 
chains.) We compare these moments in Fig. |28] with Eq. ( 177b 
tracing the cummulant Kp{s) defined in Eq. ( 126b as a func- 
tion of 

3{2Pplpf c, 

2(2p+l)! ^s 

All data sets collapse nicely on the prediction (bold line) for 
small x|3 It is important that the same effective bond length 
b is obtained from the analysis of all functions Kp(s) as illus- 
trated in Fig. |28] Otherwise we would regard equilibration 
and statistics as insufficient. 

The failure of Flory's hypothesis can also be demon- 
strated by means of the non-Gaussianity parameter Up (s) de- 
fined in Eq. ( 1271) which compares the 2/7-th moment with the 



^^ The clearly visible plateau can be used to determine the coefficient 
c*. The same coefficient can then be used for fitting the large-r behav- 
ior of monodisperse systems. Note that the best fit value c-^ = 0.14 
is close to P) (i- = 1 ) = (/„ • /j,^! ) //^ Ri 0. 10, the independently deter- 
mined bond-bond correlation between adjacent bond vectors. 

^ The curvature of the data at small s is due to the finite extensibility 
of the subchains which becomes more marked for higher moments. 



small s ! 




Fig. 29 Non-Gaussianity parameter ap{s) computed vs. Cs/y'i. Per- 
fect data collapse for all chain lengths A' indicated is obtained for each 
p. A linear relationship over nearly two orders of magnitude is found 
as theoretically expected, Eq. <78t . The data curvature for small s be- 
comes more pronounced for larger p. 



second moment (p = 1). In contrast to the related function 
Kp{s) this has the advantage that here two measured prop- 
erties are compared without any tunable parameter, such as 
b. Fig. |29] presents (Xp{s) vs. c^/\/s for three moments. For 
each p we find perfect data collapse for all A^ which con- 
firms the expected linear relationship Cip{s) w M*(i)|2j The 
lines indicate the theoretical prediction Eq. (178b . The pref- 
actors 6/5, 1 1 1/35 and 604/105 for p = 2, 3 and 4 respec- 
tively are nicely confirmed. They increase strongly with p, 
i.e. the non-Gaussianity becomes more pronounced with in- 
creasing p. Hence, b should be best fitted by the second mo- 
ment where the non-Gaussian behavior is the weakest. 

Figure [30] presents a similar correlation function which 
measures the non-Gaussian correlations of different spatial 
directions. It is defined by Kxy{s) = \ — (x^ y^) / (x^^ (y^) 
for the two spatial components x and y of the vector r as il- 
lustrated by the sketch given at the bottom of Fig. [30] Sym- 
metry allows to average over the three pairs of directions 
(x,y), (x,z) and (x,z). Following the general scaling argu- 
ment given in Sec. [2] we expect Kxy{s) w u*{s) w Cs/y/s 
which is confirmed by the perturbation result 



K„{s) 






(127) 



This is nicely confirmed by the linear relationship found 
(bold line) on which all data from both simulation models 
collapse perfectly. The different directions of subchains are 
therefore coupled. As explained at the end ofAppendix lB.il 
A'vv (s) and a2 {s) must be identical if the Fourier transformed 
subchain size distribution G{q,s) can be expanded in terms 
of q^ and this irrespective of the values the expansion coeffi- 
cients take. Fig. [30] confirms, hence, that our computational 
systems are perfectly isotropic and tests the validity of the 

^' If one plots a,,{s) as a function of the r.h.s. of Eq. l l78t all data 
points for all moments and even for too small s collapse on one master 
curve just as in Fig. l l28t . 
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Fig. 30 Plot of Kj,y{s) = 1 — (x^y^) / (x^) (y^) averaged over all pairs 
of monomers (n,m = n + s) and three different direction pairs as a func- 
tion of Cs/yJ. As indicated by the sketch at the bottom of the figure, 
Kjiy{s) measures the correlation of the components of the subchain vec- 
tor r. All data points collapse and show again a linear relationship 
Kxy ~ M*. Different directions are therefore coupled! No curvature is 
observed over two orders of magnitude confirming that higher order 
peiturbation corrections are negligible. Noise cannot be neglected for 
large s > 100 and finite subchain-size effects are visible for s k^ I. 
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Fig. 31 Deviation SG{r,s) = G{r,s) — Go{r,s) of the measured sub- 
chain size distribution from the Gaussian behavior Go(n.s) for several 
s ^N and A' = 2048. As suggested by Eq. l l75t . we have plotted y = 
{5G(r,s)/GQ{r,s))/(cs/\/s) as a function of « = r/h-/s. The data col- 
lapse confirms that the deviation scales linearly with w* rs Cs/v^. The 
bold line indicates the universal function f{n) predicted by Eq. <74t . 



general analytical expansion. The correlation function K^y 
is of particular interest since the zero-shear viscosity should 
be proportional to (a^,,) -- (x^y^) = {^) (j^) (1 - Kxy{s)) 
where we assume following Edwards |l4]| that only intrachain 
stresses contribute to the shear stress Oxx-. Hence, our results 
suggest that the classical calculations [4J — assuming incor- 
rectly Kxy = — should be revisited. 



5.6 Corrections to the subchain size distribution 

We turn finally to the subchain size distribution G{r,s) itself 
which is presented in Fig. [3T| From the theoretical point of 
view G(r, i) is the most fundamental property from which all 
others can be derived. The normalized histograms G{r,s) are 
computed by counting the number of subchain vectors be- 
tween r — dr/2 and r + dr/2 with dr being the width of the 
bin and one divides then by the spherical bin volume. Since 
the BFM model is a lattice model, this volume is not inr^dr 
but given by the number of lattice sites the subchain vector 
can actually point to for being allocated to the bin. Incorrect 
histograms are obtained for small r if this is not taken into 
account. Averages are taken over all (sub)chains, just as be- 
fore. Clearly, non-universal physics must show up for small 
vector length r and small curvilinear distance s and we con- 
centrate therefore on values r ^ a and 5 > 3 1 . When plot- 
ted in linear coordinates as in Fig. 11 of 0841 . G{r,s) com- 
pares roughly with the Gaussian prediction Go(r,5) given by 
Eq. (|2|, but presents a distinct depletion for small subchain 
vectors with n = r/by/s <C 1 and an enhanced regime for n w 
1 . To analyse the data it is better to consider instead of G{r,s) 
the relative deviation dG{r,s)/GQ{r,s) = G{r,s)/Go{r,s) — 1 
which should further be divided by the strength of the sub- 
chain correlation hole, Cs/\/s. As presented in Fig. 13 II this 
yields a direct test of the relation Eq. ( 175b derived in Ap- 
pendix IB .21 The figure demonstrates nicely the scaling of 
the data for all s. It shows further a good collapse of the 
data close to the universal function f{n) predicted by theory 
(bold line). Note that the depletion scales as l/« for small 
subchain vectors (dashed line). The agreement of simulation 
and theory is by all standards remarkable. Obviously, error 
bars increase strongly for n ^ 1 where Go{r,s) decreases 
strongly. The regime for very large n where the finite exten- 
sibility of subchain matters has been omitted for clarity. We 
emphasize that this scaling plot depends very strongly on 
the value b which is used to calculate the Gaussian reference 
distribution^ 



5.7 Intramolecular form factor 7^(17) 



The form factor F{q) defined in Sec. 12.2.41 is an important 
property since it allows to make a connection between the- 
ory and simulation on the one hand and experiments of real 
systems on the other hand |5|. Figure [32] presents the (un- 
sealed) form factors obtained for four different scission en- 
ergies E for our EP model at (/) = 0.5. The three differ- 
ent ^-regimes are indicated. Details of the length distribu- 
tion /7n matter in the Guinier regime which probes the to- 
tal coil size. Non-universal contributions to the form factor 
arise in the "Bragg regime" at large wavevectors. Obviously, 
the larger E the wider the intermediate Kratky regime where 

^^ As shown in Ref. 1841 . a similar plot can be achieved which does 
not require a precise value of b if the reduced deviation y is plotted as 
a function of r/Rs using the measured subchain size R^. 
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Fig. 32 Intramolecular form factor F{q) vs. wavevector q for EP of 
various scission energies E 1 83 1. The ideal chain form factor for Flory- 
distributed polymers, Eq. l l40t . is indicated by the solid line. In the 
Kratky regime between the total chain and monomer sizes the form 
factor expresses the fractal dimension of the Gaussian coil, Eq. <35K as 
shown by the dashed line. Experimentally, this is the most important 
regime since it is, e.g., not affected by the (a priori unknown) poly- 
dispersity. The computational data reveal an additional regime at large 
wavevectors corresponding to the monomer structure ("Bragg regime") 
which is not treated by the theory. 




Fig. 33 Kratky representation of F{q)q^ vs. q for monodisperse 
quenched polymers (crosses) and equilibrium polymers (open sym- 
bols). The non-monotonous behavior predicted by the theory is clearly 
demonstrated. The ideal chain form factor (thin line), overpredicts the 
dip of the form factor lAq^i 0.7 by about 20%. The bold line indicates 
the prediction for infinite chains, Eq. l lT4T i. The data is fitted assuming 
the effective bond length b ^ 3.244. 



chain length, polydispersity and local physics do not con- 
tribute much to the deviations of the form factor from ideal- 
ity. A very similar plot has been obtained for monodisperse 
polymers (not shown). Not surprisingly, it demonstrates that 
the form factors of both system classes become indistin- 
guishable for large wavevectors. Note that the F{q) obtained 
for EP of different E can be brought to collapse by tracing 
F{q)/F{0)asa function ofQ = qRg_^ with F (0) = (N^) / (N) 
and /?g z the measured z-averaged radius of gyration as indi- 
cated in Table [3] 1831 . A similar plot can again be obtained 
for monodisperse chains using F{0) =N and Q = qRg{N). 



Interestingly, a careful inspection of Fig.[32]reveals that 
Eq. ( 1401) overestimates systematically the data in the Kratky 
regime. This can be seen more clearly in the Kratky repre- 
sentation given in Fig.[33]in linear coordinatesli] We present 
here the systems with the longest masses currently available 
for both monodisperse {N = 4096 and Af = 8192) and EP 
systems {E = 14 and E = 15) II83L The non-monotonous 
behavior is in striking conflict with Flory's hypothesis. The 
difference between the ideal Gaussian behavior (thin line) 
and the data becomes up to 20%. This difference is qualita- 
tively expected from the perturbation theory result, Eq. ( 114b . 
derived in Appendix IB. 31 The predicted scaling can be eas- 
ily understood by means of a simple scaling relation follow- 
ing the discussion in Sec. 12.3.61 As the arc-length 5 is re- 
lated to the wavevector q and the ideal form factor Fo(?) by 
.v(^) ^ 1/I^P ^ Fo{q), it follows for the difference of mea- 
sured and ideal form factors that 1.82] 



F{q) 






1 * 

- XM* 
5 



1^ (128) 
P 



which agrees with Eq. ( 1141) ri The important feature of this 
l^p-correction is that it depends neither on the strength of 
the excluded volume interaction nor on the effective bond 
length b. Hence, it must be generally valid, even for semidi- 
lute solutions rn The prediction for infinite A'^, Eq. ( 114b . can- 
not capture the decrease of the form factor for small q lead- 
ing to the Guinier regime where F{q) is determined by the 
finite size of the simulated chains. A clearer evidence for the 
theory should thus be obtained by a different comparison be- 
tween theory and simulation, which accounts for the finite -A'^ 
effects. 

This is achieved in Fig. [34] which focuses on deviations 
5{\/F{q)) = \/F{q) — \/FQ{q) for monodisperse systems 
with FQ{q) being the Debye formula. As above in our dis- 
cussion of the subchain distribution G{r,s) and its moments 
one should avoid to use as Gaussian reference the ideal chain 
form factor F{){q) expressed in terms of the effective bond 
length b since the latter property might not be sufficiently 
accurate. A variation of a few percents breaks the scaling 
and leads to qualitatively different curves 183 1 as in Fig. |24] 
for the moment K\{s). Since such a precision is normally 
not available (neither in simulation nor in experiment) it is 
of interest to seek for a more robust representation of the 

^^ Such "Kratky plots" are commonly used to represent neutron scat- 
tering experiments and to test Flory's ideality hypothesis. In applica- 
tions, the existence of a Kratky plateau appears to be elusive |5|. Pos- 
sible causes for deviations are effects of chain stiffness and finite chain 
thickness. In special cases, the scattering signals from chain stiffness 
and thickness may comp ensate one another, leading fortuitously to an 
extended Kratky plateau I167II68I . Apparently, "Kratky plots have to 
be interpreted with care" [Sf 

^'* This is the same scaling dependence which leads to Pi{r) ~ 
1/p |r|'' ~ 1/p (bs'^I^Y K, Pi (s) for the bond -bond correlation function 
as a function of r or s as investigated in Sec. 15.4.11 

^' It applies then for g <g; 1 /(^ with £, being the semidilute blob length 
|3|. We have checked that the result from the renormalization group 
theory for semidilute solutions [7, 71 ] takes the same form as Eq. i]M 
with an amplitude 0.03124 that is within 0.03% of our 1/32 [83 1. 
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Fig. 34 Scaling attempt of the non-Gaussian deviations for monodis- 
perse polymers in terms of tiie measured radius of gyration Rg {N) . As 
suggested by Eq. l ll4t . tiie difference \/F(q) — l/Fo(?) of tiie mea- 
sured and the ideal chain Debye form factor has been rescaled by the 
factor Np/p^ and plotted as a function of 2 = '?^g- We obtain perfect 
data collapse for all chain lengths included. (Obviously, data points in 
the Bragg limit </ w 1 do not scale.) Note that the power law slope, 
'"(2) = 2^/32, predicted by Eq. SlAJ . can be seen over more then 
one order of magnitude. In the Guinier regime, the difference increases 
more rapidly, m(Q) ~ Q^ (dashed line), as one expects from a standard 
analytic expansion in Q^. 



form factor deviations which does not rely on b. The refer- 
ence chain size is thus set in Fig. |34] by the measured ra- 
dius of gyration Rg{N) which is used for rescaling the axis 
and, more importantly, to compute Fo{q). (A virtually indis- 
tinguishable plot has been obtained for EP.) The scaling of 
the vertical axis is suggested by Eq. ( 114b which predicts the 
difference of the inverse form factors to be proportional to 
A'"^^. Without additional parameters {Rg is known to high 
precision) we confirm the scaling of 



m{Q) 



N 



N 



F{q) Fi^){q))p^ 



(129) 



as a function of Q = qRg (N) with p^ = N/R^ being the self- 
density. Importantly, our simulations allow us to verify for 
2^5 the fundamentally novel Q^ behavior of the master 
curve predicted by Eq. (fT4] | and this over more than an or- 
der of magnitude! In this representation we do not find a 
change of sign for the form factor difference {dF(q) is al- 
ways negative) and all regimes can be given on the same 
plot in logarithmic coordinates. In the Guinier regime we 
find now in{Q) <== Q^ which is readily explained in terms of 
a standard expansion in Q^ since the first two terms in Q" 
and Q^ must vanish by construction because of the defini- 
tion of /?g(Ar), Eq. 



5.8 Summary 
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Fig. 35 The displacement auto-correlation function C{t) of dense 
overdamped colloids is known to reveal a negative algebraic long-time 
tail due to long-range dynamical interactions caused by the weak com- 
pressibility of the solution |11]: If aparticle (filled sphere) is displaced 
at time f = by a distance m(0) along the x-axis, this creates a density 
dipole (bold peaks) decaying in time by cooperative diffusion (dashed 
line). The gradient of the chemical potential associated to the decay- 
ing field generates a drag pulling the particle (open sphere) back to its 
original position. We argue 1 90 1 that a related mechanism causes scale- 
free correlations of the displacements of (sub)chains in dense polymer 
solutions without topological constraints. 



to the predicted deviations with respect to Flory's ideality 
hypothesis. From the computational side the most important 
point is discussed in Sec. 15. 3. 2] where we show how the ef- 
fective bond length b for asymptotically long chains should 
be extrapolated from subchains of finite arc-length with g ^ 
.s ^ A'^. That the deviations are indeed due to long-range 
interactions has been demonstrated from the scaling of the 
bond-bond correlation function Pi (r) ^ \/P in Sec. 15.4.2] 

The most important finding from the experimental side 
concerns the scale-free deviations demonstrated for the form 
factor F{q), Eq. ( fT4] l. We have shown that the Kratky plot 
does not exhibit the plateau expected for Gaussian chains in 
the intermediate wavevector range. These deviations should 
be measurable by neutron scattering experiments of flexible 
chains. Unfortunately, finite persistence length effects may 
mask the predicted behavior if the chains are not sufficiently 
long. An experimental verification — e.g. following the lines 
of the the promising recent study fl681 — would be of great 
fundamental interest; it could also delineate the conditions 
where the predicted conformational corrections to ideality 
are relevant in real polymer systems and must be considered 
in understanding their structure, phase behavior and equi- 
librium dynamics. It is to the latter point we turn now our 
attention. 



6 Scale-free dynamical correlations in polymer melts 

6.1 Introduction 

6.1.1 Prelude: Overdamped colloids 



Dense, essentially incompressible simple liquids with con- 
served momentum are known to exhibit long-range corre- 
In this section we have investigated various intrachain static lations of the particle displacement field 1 1 1 , 169 1. As first 
properties of polymer melts comparing our numerical data shown in the MD simulations by Alder and Wainwright II169I . 
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the coupling of displacement and momentum fields mani- 
fests itself by an algebraic decay of the VCF, 



c(0 = (v(0-v(0))^+i/^d(0'-^+iA 



ad 



(130) 



with vit) being the particle velocity at time t and ^i{t) ^t'^ 
the typical particle displacement with exponent a = 1 /2|_3 
Interestingly, even if the momentum conservation is dropped, 
as justified for overdamped dense colloidal suspensions 111], 
scale-free albeit much weaker correlations are to be expected 
due to the incompressibility constraint II11II170I . As illus- 
trated in Fig. |35] the motion of a tagged colloid is coupled 
to the collective density dipole fieldO 



8p{r,t)'. 



1 u{0)-r 



exp[-(r/^d(0)'] for f> 0,(131) 



created by the colloid's own displacement m(0) at f = 0. Af- 
ter averaging over the typical displacements of the test par- 
ticle and assuming a Cahn-Hilliard response proportional to 
the gradient V/j (r) of the chemical potential jU (r) of the den- 
sity field 111, this leads to a negative algebraic decay 



C(t) 



-1/^d 



\d+2 



l/t'^wiiho) = id + 2)a, (132) 



e.g., ft) = 5/2 in d = 3 dimensions. This phenomenological 
scaling picture agrees with more systematic mode-coupling 
calculations lilllll70l[r7T1l . It has been confirmed computa- 
tionally by means of Lattice-Boltzmann simulations linil . 
MD simulations 01721 and even M C sim ulations with local 
moves I90II as demonstrated in Sec.l6.2f 



6.1.2 Polymer melts without topological constraints 

The finding of long-range static correlations in dense poly- 
mer solutions discussed in Sec.|5]begs the question of whe- 
ther a similar interplay of chain connectivity and incom- 
pressibility leads to similar scale-free and ^-independent dy- 
namical correlations. To avoid additional physics we focus 
on model systems where hydrodynamic f92ll and topolog- 
ical constraints may be considered to be negligible, as in 
the pioneering work by Paul et al. [501(51], or are deliber- 
ately switched off II54II55| [841. As we have summarized in 
Sec. 11.51 deviations from Rouse-type dynamics have been 
reported for such systems in various numerical and experi- 
mental studies. The effective exponent j3 w 0.8 characteriz- 
ing the short-time CM diffusion, Eq. (fT3T i. is of course in- 
consistent with the key assumption of the Rouse model that 
the random forces acting on the chains are uncorrelated. We 
attempt here to clarify this problem using the BFM variant 
with topology non-conserving local L26-moves described in 
Sec. [3j Obviously, it would not be possible to equilibrate 

^^ The discussion is strongly simplified. Strictly speaking, it is not 
the particle diffusion which sets the dynamical l ength scale but the 
diffusive propagation of the transverse momentum 11691 . 

^' Prefactors are omitted for simplicity. Especially, we do not dis- 
tinguish between the self diffusion of the particle and the collective 
diffusion of the field characterized by rather different coefficients ILUJ. 

^^ Interestingly, the same power-law exponent (O = (d + 2)/ 2 is a lso 
seen in Brownian dynamics simulations of the Lorentz model 11731 . 



BFM configurations with chain lengths up to A' = 8192 us- 
ing local moves. Taking thus advantage of the configurations 
obtained using global MC moves we demonstrate that the 
VCF Cn(0 associated to the chain CM displacements does 
not vanish, but instead decays as 



Cn(0 « - 



7n 



Pn 



/(f/7N) 



(133) 



with f{x) being a universal scaling function. Note that the 
postulated Eq. ( 1133b does not depend explicitly on the com- 
pressibility of the solution. The squared characteristic "ve- 
locity" (/?n/7n)^ arises for dimensional reasons. The sec- 
ond pref actor p^/p stems from the interaction of chains and 
subchains imposed by the incompressibility constraint^ As 



one expects from Eq. ( 1132b . the scaling function decays for 
x> 1 asf{x)^ 1/x™ with an exponent ft) = ((i-f- 2) /2. More 
importantly, it will be shown that this long-time behavior is 
preceded for jc <C 1 by a much weaker algebraic decay with 



ft) = («i + 2)a = 5/4 for fi? = 3 



(134) 



due to the much slower relaxation (a = 1 /4) of the collective 
dipole field of subchains which was generated by the initial 
displacement of the tagged subchain at f = 0. The gradient 
of the chemical potential V/i (r) pulling the reference sub- 
chain back to its original position is of course not only due to 
the density fluctuation of the subchain density field but also 
to the tensional forces along the chains caused by the dis- 
placement. However, since subchain density fluctuations and 
tensions are coupled, their free energy contributions cause 
chemical potential gradients of the same order IcbTs'^N^ and 
are thus identical from the scaling point of view [;90II91| . 

6.1.3 Outline 

Reminding first how a VCF may be defined and determined 
in a Monte Carlo simulation [174], we confirm in Sec. 



that Eq. ( fT32] l holds for dense BFM beads. Section 163] fo- 
cuses then on polymer melts demonstrating the dynamical 
coupling of (sub)chains. Excluded volume and density ef- 
fects are discussed in Sec. 16.41 Preliminary simulation data 
for melts with topological constraints presented in Sec. 16.51 
indicate that the early-time behavior, Cn (0 '^ — A^ ^ i j - (^+2)/4 
for t <C 7n, is preserved before entanglement effects set in. 
A perturbation calculation prediction for this time window is 
outlined in Sec. 16.61 In demonstrating Eqs. ( 1133b and ( 11341) 
our study focuses on one mechanism explaining the strik- 
ing deviations from the Rouse behavior observed in the lit- 
erature. See Refs. [175,99,176 92 93] for related theoretical 
studies. We stress that the presented MC simulations ||90]l are 
necessarily incomplete since important additional dynami- 
cal correlations arise due to the incomplete screening of hy- 
drodynamic interactions mentioned briefly in Sec. 17.2.31 192] 

^^ That the effect is proportional to the correlation hole penalty, wj, w 
Pn/P^ may be guessed from the static coiTelations we have discussed in 



Sec. 12 The fact that the polymers behave as weakly interacting colloids 
for large times t/T^ 2> 1 confirms this scaling. 
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Fig. 36 Diffusion of BFM beads (N = I) without monomer overlap 
(e = oo) in d = 1 + , rf = 2+ and d = 3 dimensions obtained using L26- 
moves at the reference "melt density" (^ = 0.5) 1 174'|. Inset: Although 
the bead diffusion is essentially free, small deviations are visible if 
h{t) / (2dDt) is plotted in log-linear coordinates. Main panel: Collapse 
of displacement correlation functions C{t, St) for rf = 1+ and various 
time increments 5t (open symbols). The "final" function C(f) is ob- 
tained by adding the first decade of data for each St and logarithmic 
averaging. The cummulants for each effective dimension a gree n icely 
with the predicted power-law exponent a = {d + 2)/2, Eq. ( I132t . 



6.2 Diffusion of dense BFM beads 

Since in MC simulations there is no "monomer mass", no 
(conserved or non-conserved) "monomer momentum" and 
not even an instantaneous velocity, it might at first sight ap- 
pear surprising that a well-posed "velocity correlation func- 
tion" (VCF) can be defined and measured. To illustrate that 
this is indeed the case is the first purpose of this subsection. 
The second is to verify that the negative analytic decay of the 
VCF expected for overdamped colloids, Eq. ( 11321) . is also of 
relevance for dense BFM beads (A'^ = 1) diffusing through 
configuration space by means of local hopping moves on 
the lattice as shown in Fig. lUb). The systems presented in 
Fig. [36] correspond to three different effective dimensions 
d. The effectively ID systems {d = 1+) have been obtained 
by confining the beads to a thin capillary of square cross- 
section, the effectively 2D systems (d = 2+) by confining 
the beads to a thin slit as shown in Fig. |6jc). The distance 
H = 4 between parallel walls was chosen to allow the free 
crossing of the beads. (For the ID case this is important, 
of course, since ordering the beads along the capillary al- 
ters dramatically the dynamics.) The data has been obtained 
for beads without monomer overlap (e = oo) by means of 
L26-moves on a 3D cubic lattice with half of the available 
lattice sites being occupied ((/) = 0.5). We average over the 
2^" beads contained in each configuration. 

One standard measure characterizing the monomer dis- 
placements is the mean-square displacement (MSD) h{t) = 
((l(f) -£(0))^) displayed in the inset of Fig. [36] (with r{t) 
being the particle position at time t). As one expects, the 
monomer displacements become uncorrelated for f ^ 1, i.e. 
h{t) w 2dDt with D being the monomer self-diffusion con- 
stant: D = 0.0187 for d=l+,D = 0.0154 for d = 2+,D = 



0.0382 for d = 3. The typical displacement ^d(0 ^ /j'/^(0 
thus scales as ^d{t) ^ t" with a = 1/2. However, small de- 
viations are clearly visible for short times if h{t) / {2dDt) 
is plotted in log-linear coordinates. The deviations are par- 
ticular strong for d = l^ . The effective random forces act- 
ing on the beads are thus not completely white. Obviously, 
one might try to characterize these deviations by fitting vari- 
ous polynomials to the measured MSD. However, since one 
needs to subtract the rather large free diffusion contribution 
from the measured signal to characterize tiny deviations this 
is a numerically difficult if not impossible route. 

In analogy to the bond-bond correlation function Pi (s) ^ 
d^R^ discussed in Sec. 15.41 allowing to make manifest de- 
viations from the Gaussian chain assumption, it is numeri- 
cally much better to directly compute the second derivative 
of h{t) with respect to time to avoid this large, but trivial 
contribution. How this can be done is illustrated in the main 
panel of Fig. [36] We sample equidistant series of configu- 
rations at time intervals dt = 1, 10, 100, ... as indicated by 
the open symbols. Each time series contains lO'* configu- 
rations. Averaging over all possible pairs of configurations 
(fo,fo + we compute C(f, 5t) = {u{to + t) ■ u{to)), /5t^, i.e. 
a four-point correlation function of the monomer trajectories 
with u{t) — r{t + 5t) — r{t) being the monomer displace- 
ment vector at time f in a time interval dt. By construction 
C{t, dt) = if both displacement vectors are uncorrelated. 
According to Eq. (1173b one expects 

1 



C(f , 5t) w -dfh{t) dt" for t > dt. 



(135) 



As can be seen for BFM beads in J = l"*" dimensions, the 5t- 
dependence indeed drops out for t/dt > 1 and we thus avoid 
the second index dt writing C{t) for the displacement (or 
velocity) correlation function. Obviously, the statistics dete- 
riorates for very large t/dt where fewer configuration pairs 
contribute to the average (taking apart that the signal itself 
decays). It is for this reason that a hierarchy of time series 
of different dt is needed. Taking for each dt only the first 
decade of data {2 < t / dt < 20), these data sets are pasted to- 
gether and averaged logarithmically. The exponents (0 = 3/2 
(thin line), (0 = 2 (dashed line) and (0 = 5/2 (bold line) pre- 
dicted by Eq. ( 1132b ioi d = l,d = 2 and d = 3, respectively, 
compare well with our data over several orders of magni- 
tude in time, especially for d = 1'^. Note that if one is sat- 
isfied with less orders of magnitude it is sufficient to check 
the exponents using just a time window 5f = 1 as may be 
seen from the open spheres. The superposition of data from 
different time series is just a numerical trick which reduces 
the number of configurations to be stored and the number of 
configuration pairs to be computed for a given time t. 



6.3 Polymer melts without topological constraints 

6.3.1 Mean-square displacements 

Having shown that scale-free dynamical correlations exist 
for dense BFM beads as expected for overdamped colloids 
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Fig. 37 Various MSD for chains of lengtlis JV = 64andJV= 8192 t90l. 
The open symbols refer to the monomer MSD hit), the filled symbols 
to the MSD /jn(0 of the CM of chains, the crosses to the MSD h^(t) 
of the CM of subchains of arc-length s = 64 of total chains of length 
N = 8192. The dashed line indicates the monomer MSD expected for 
Rouse chains for ? <g; Tn, the thin solid lines the free diffusion limit. 
As emphasized by the bold lines corresponding to the exponent /3 = 
0.8 suggested by Eq. <13t . coixelations are visible for th e short-time 
behavior of the CM motion. The stars indicate Eq. <141t using c = 1 , 
i.e. the time window f -C ?* ~ 10^^ is described by an exponent /3 = 3/4. 
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Fig. 38 Acceptance rate A, effective local mobility W and self- 
diffusion coefficient ND^ vs. the monomer overlap penalty e. The data 
has been obtained using L26-moves at volume fraction <^ = 0.5 |90|. 
All dynamical properties decrease monotonously with the interaction 
penalty but become essentially constant above our main working point 
at e = 10. The dashe d line corresponds to the diffusion coefficient ac- 
cording to Eq. l ll37t and assuming the indicated mobilities W . 



(Fig. [35T l, we turn now our attention to 3D melts of long and 
flexible homopolymers. As described in Sec. 13.41 we focus 
on systems with finite overlap penalty £ = 10 and volume 
fraction = 0.5 which we sample using local topology non- 
conserving L26-moves. As can be seen in Fig.[32]for chains 
of length A'^ = 64 and A^ = 8192, these systems are essen- 
tially of Rouse-type. The monomer MSD h(t) is indicated 



by the open symbols|_3 As expected from Rouse dynamics 
we obtain the A^-independent short-time asymptotics ID 

h{t) = b^iWtyi^N^ for 10 < f < Tn^N'^/W (136) 

indicated by the dashed line (2a = 1/2). As can be seen for 
A^ = 64, the monomers diffuse again freely with a power-law 
slope 1 (thin lines) for times larger than the Rouse time T^. 
We remind that it was neither computationally feasible nor 
our goal to sample for our larger chains (A^ > 1000) over the 
huge times needed to make this free diffusion regime acces- 
sible. Following Paul et al. II51L the short-time power law, 
Eq. ( 11361) , can be used to determine the effective monomer 
mobility. We obtain W{e = 10) = 0.003 for our main work- 
ing point. Mobilities for other penalties are listed in Table [T] 
As may be seen from Fig.[38j W(e) decays with increasing 
excluded volume just as the acceptance rate A (e) (spheres) 
of the Monte Carlo attempts, but the decay is even more pro- 
nounced for W(e), i.e. increasingly more accepted monomer 
moves do not contribute to the effective motion [SI]. 

The full symbols displayed in Fig. [37] refer to the MSD 
h-nit) = ((z:n(0~Z:n(0))^) of the CM i^{t) of chains of 
length A'^. As one expects for Rouse chains, the amplitude of 
hfi {t ) decreases inversely with A'^ and the diffusion appears 
to be Fickian (/jn(0 ^ t) at least for times f > f * w 10"^ A^" as 
indicated by the vertical arrow. Obviously, /jn(0 and h{t) 
merge for times beyond the Rouse time {t ^ Tn). Fortu- 
nately, since /jn(0 becomes linear for t* <C 7n, it is possi- 
ble for all A^ to measure the self-diffusion coefficient Z)n by 
plotting Nh^ {t ) /6t vs. time t. For our main working point we 
obtain ND^{e = 10) w 0.009. Values for other e are again 
given in Table[T]and are represented in Fig. 1381 (diamonds). 
These values compare nicely with the Rouse model predic- 
tion 14] 



n 



NDj^ = -a"- W 



(137) 



shown by the dashed line in Fig. [38] Similarly, it is possi- 
ble (at least for our shorter chains) to measure the longest 
Rouse relaxation time T^ by an analysis of the Rouse modes 
and to compare it with the Rouse model prediction |4 1 Tn = 
4N^/n^W. We obtain again a nice agreement between di- 
rectly and indirectly computed relaxation times (not shown). 
Up to now we have insisted on the fact that our systems 
are to leading order of Rouse type and we have character- 
ized them accordingly. However, deviations from the Rouse 
picture are clearly revealed for short times, especially for 
/JnCOI— I In agreement with Eq. ( [T3T i the short-time CM mo- 
tion may be characterized by a power-law exponent j5 w 
0.8 (bold lines). Since in our BFM version topological con- 
straints are irrelevant, this shows that the deviations obtained 
for the classical BFM algorithm with topological constraints 
[51 1 cannot be attributed alone to precursor effects to repta- 
tional dynamics (which indeed exist as shown in Fig.| 



''* For the large N sampled here it is inessential whether this average 
is performed over all monomers or only ov er a few monomers in the 
center of chains at / w N/2 as in 15015 11561 . 

'^ That the monomer MSD h{t) also deviates for very short times is 
due to the trivial lower cut-off associated to the discretization. 
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Before we turn to the more precise numerical character- 
ization of these deviations by means of the associated dis- 
placement correlation function, let us ask whether the ob- 
served colored forces acting for short times on the CM of 
the entire chain are also relevant on the scale of subchains 
of arc-length i (1 ^ i < A^). To answer this question we 
compute the MSD hs{t) = ((/:s(0 - Lsi^))^) associated to 
the subchain CM r\.{t) as shown in Fig. (37] for subchains of 
length i = 64 in the middle of total chains of length A^ = 
8192 (crosses). Since for short times the subchain does not 
"know" that it is connected to the rest of the chain, one ex- 
pects it to behave as a total chain of the same length (s = N). 
This is indeed borne out by our data which are well de- 
scribed by 



o 8t=i 

D 8t=10 

O 8t=100 

A 8t=1000 

< 8t=10000 



sh^{t) w Nh^{t) for f < T^ w s^ /W 



(138) 



for all chain length A' and subchain length i studied. The 
subchains reveal thus for sufficiently short times the same 
colored forces as the total chain as can be clearly seen from 
the example given in Fig. (37] For larger times the subchain 
becomes "aware" that it is connected to the rest of the chain 
and gets enslaved by the monomer MSD, i.e. as expected 
from the Rouse model we observe 

h,{t) w h{t) w b^{Wt)^/^s^N° for T^ < f < Tn (139) 

and, obviously, hs{t) w ^^(f) w h{t) w 6D^t for even larger 
times t > Tn. 

6.3.2 Locality and relevant exponent a 

Two comments are in order here. First, it should be noticed 
that Eq. ( 1138b expresses the fact that the effective forces 
acting on the N /s subchains of length s in a chain of total 
length A'^ add up independently to the forces acting on the 
total chain. In this sense Eq. ( 1138b states that the deviations 
from the Rouse picture must be local. We will explicitly ver- 
ify this below (Fig. |4T]l[__| Second, if one chooses follow- 
ing Eq. ( 1135b an arbitrary time window 5t to characterize 
the displacement correlations, this corresponds to dynami- 
cal blobs containing i w ^JW8t adjacent monomers, which 
must move together due to the chain connectivity. Eq. ( |139b 
implies now that the dipole fielq__| associated with the CM of 
these ,s-subchains and created at /= by a tagged 5-subchain 
must decay according to a typical displacement ^d(0 ^ t" 
with a = 1/4. It is this exponent a which is mentioned in 
Eq. (1134b . Since this exponent is smaller than for colloids 
(a = 1/2), the subchain field must decay more slowly and 
one thus expects a much weaker decay of the associated 
VCE 

^^ The "locality" of the correlations described by Eq. {T38) or 
Eq. ([143) does not imply that the displacements of subchains around 
the reference subchain displaced at f = are 5(r)-coiTelated. 

^^ We remind that, strictly speaking, it is not the dipole field associ- 
ated to the density of the subchain center-of-masses, but to their chem- 
ical potential. The use of the imprecise notion "density dipole field" 
might be excused by the fact that both fields are supposed to decay 
similarly with time t, chain length A^ and subchain length s. 




10" 10 10 

t [MCS] 

Fig. 39 VCF CN(?,5f) for two chain lengths, A^ = 16 (top) and A' = 
8192 (bottom), obtained using L26-moves |90|. We only indicate for 
each St the data window which is used to construct (by logarithmic 
averaging) the final VCF Cn(?). The bold lines represent the exponent 
CO = 5/4 which is generally observed for times 10 ^ f ^ Tn. For N = 
16 we also indicate the exponent ft) = 5/2 (dashed line) expe cted for 
larger times where the polymer coils behave according to Eq. <132t . 



6.3.3 Center-of-mass velocity correlation function 

Following the numerical strategy presented in Sec. 16.21 we 
characterize now more precisely the correlations seen for the 
chain MSD lt^{t) and the subchain MSD h^{t) by computing 
directly their second derivative with respect to time f, i.e. 
the associated correlation functions Cn(0 ~ ^r^%(0/2 and 

c,{t)^dXit)l^- 

The VCF Cn{t, 5t) = {uit + to) ■ a(fo)),^ /5t^ for the dis- 
placement vector u{t) = £n(^ + ^0 ^ £n(0 of the chain CM 
tj^it) is shown in Fig.[39]for two chain lengths, A'^ = 16 (top) 
and A^ = 8192 (bottom). Averages are again performed over 
all configuration pairs (fo,f + fo) in the set of 10^ config- 
urations sampled for each dt. As in Sec. 16.21 we find that 
Cj^{t,5t) ^ 5t'^ for t ^ 8t. For clarity, only the data subset 
is indicated for each 5t which is used to construct the final 
VCF C-^{t) (as shown below in Fig.l40l). The bold lines rep- 
resent the predicted short-time exponent (0=5 /A which can 
be observed for A'^ = 8192 over nearly five orders of mag- 
nitude. For A^ = 16 we also indicate the exponent ft) = 5/2 
(dashed line) for t ^ Tj^ where the polymers should behave 



as colloids according to Eq. ( 1132b . Note that the magnitude 
of the signal decreases strongly with A^, which together with 
the fact that fewer chains per box are available, makes the 
determination of Cj^it) increasingly more delicate. 

The A'^-dependence of the VCF Cn(0 is further analyzed 
in Fig.iO] The rescaled VCF y = -Cj^{t){b^p)/{WD^) is 
traced in panel (a) as function of the reduced time x = Wt 
using the monomer mobility W = 0.003 and the diffusion 
coefficient Dn = 0.009/A'^ determined above (Table [Hi. This 
scaling makes the axes dimensionless and rescales the verti- 
cal axis by a factor A^. As shown by the successful data col- 
lapse for chain lengths ranging from A' = 16 up to A' = 8192 
on the slope indicated by the bold line, the VCF scales ex- 
actly as Cn(0 ^ l/N for f <C 7n. This confirms the already 
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Fig. 40 Scaling of VCF Cn(0 for different A^ (90|: (a) Rescaled VCF 
y = -Ct^{t){h^p)/{WDt^) vs. reduced time x = Wt. For times t < Tn 
tiie VCF tiius clearly scales as Cf^{t) ~ l/N, i.e. the correlations must 
be due lo ~ N local independent events. The exponent a = 5/4 (bold 
line) is observed over up to five orders of magnitude in time, (b) Col- 
lapse of J = CN (f)/C 'TM with Ctn ~ 1/A''/^ as a function oi x = t/Tt^ 
confirming Eq. n33t . The large time behavior (x 2> 1) is described by 
an exponent co = 5/2 (dashed line). The exponent (O = 5/4 (bold line) 
for jc <g 1 is demonstrated over eight orders of magnitude. 



stated "locality" of the correlations, Eq. ( 11381) , i.e. the forces 
acting on subchains add up independently to the forces act- 
ing on the entire chain. We have still to motivate the pre- 
cise scaling used for the axes. According to Eq. (1133b the 
VCF scales as a function of the reduced time t/T^. Substi- 
tuting the typical chain size 7?n ~ bN^'^ and relaxation time 
7n « Rn/Dn « N^/W this reduces for (0 = 5/4 to 

y = -Ct^{t){b^p)/{WDt^) = c {Wty^l^N^ for t < 7n(140) 

with c being a dimensionless constant. The bold slope indi- 
cated in the plot corresponds to a value c = 1 . Interestingly, 
since Cn(0 ~ 4^% (0/2, it follows from Eq. ( 11401) that 



%(0 =6D^t 1 



16c 
Wp 



{Wt) 



-1/4 



(141) 



As may be seen in Fig. [37] (stars) for TV = 8192, Eq. dMll l 
with c = 1 provides an excellent fit of the measured h^{t). 
We also note that the second term in Eq. (1141b dominates the 
short-time dynamics for t <.t* = W'^ {\6c /9b'^ p)^ N° w 
10^. This is indicated by the vertical arrow in Fig. [37] Hence, 
for t -^ t* the stars correspond to a power-law slope with 
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Fig. 41 Cs(f) for subchains obtained for chains of length N = 1024 
using an overlap penalty e = 10 and L26-moves [90]. Main panel: For 
t^Ta the subchains scale as the total chains in agreement with Eq. ([12} 
as shown by the bold line. For intermediate times Ts <g; f <g; Tn the sub- 
chain displacements (1 '^s^N) follow the monomer MSD and thus 
Cs(f) ~ s°{bWf{Wty/^'^ (dash-dotted line). Inset: The latter regime 
can be better seen by setting T^ = s~/W and Cj^ = —{bW)^s^^ and 
plotting y = Cs{t)ICj^ as a function of x = t/T^. The dashed slope cor- 
responds to y = 1 /Sjc^/^ for X 2> 1 . 



exponent j3 = 2 — O) = (6 — fif)/4 = 3/4. This is close to the 
phenomenological exponent /3 = 0.8 from the literature. The 
central advantage of computing the VCF Cn(0 lies in the 
fact that it allows us thus to make manifest that (negative 
algebraic) deviations from the Rouse behavior exist for all 
times and not just for t <^t*. 

Returning to our discussion of Fig. [40l a) we emphasize 
that the VCF of shorter chains decay more rapidly for large 
times following roughly the exponent (0 = 5/2 expected for 
effective colloids. We verify explicitly in Fig. l40t b) that the 
bending down of the data is consistent with the announced 
scaling in terms of a reduced time x = t/Tj^ and a verti- 
cal axis y = C^{t)/Cj^ using the amplitude Cj^ = CN(f = 
Tn) w -(/?N/7N)2pt^/p - -I/nV^ stated in Eq. ([T33] l. The 
successful data collapse confirms that the only relevant time 
scale in this problem is the relaxation time 7n for which the 
deviations for short (CO = 5/4) and long times (ft) = 5/2) 
match. It is worthwhile to emphasize that the general scal- 
ing Eq. ( 11331) together with the locality of the deviations, 
Cn(?) ^ l/N, immediately imply the exponent ft). This can 
be seen by counting the powers of the chain length A'^, 

-l = (l/2-2)2+(l-£//2)+2ft), (142) 



thus ft) = (li + 2)/4 = 5/4 in agreement with Eq. ( 11341) and 
the numerically observed time dependence. Assuming that 
Eq. (1133b holds, the exponents for A' and t thus contain the 
same information. 

The VCF Cs(f ) for subchains of arc-length .s < A^ is pre- 
sented in Fig. [41] The subchain VCF C^{t) is obtained as 
the total chain VCF Cn(0' the only difference being that the 
CM £s(f) of the subchain defines now the displacement vec- 
tor u{t) = !\i{t + dt) — rj^{t) . The data presented in the main 
panel is rescaled as in Fig. I40t a) with D^ = 0.009/.S setting 
now the relevant diffusion constant used for the vertical axis. 
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Larger wavevectors (1/7?n ^ ^ <C 1/^) are experimen- 
tally more readily accessible |5|. Since in this regime the 
dynamical form factor probes the MSD of subchains of arc- 
length i-- 1 la^ with K (t ) replacing /zn (f ) in Eg . d 1 46b I 177II . 
it is of relevance to verify whether it is possible in practice to 
confirm the exponent j3 = 3/4 using the numerically com- 
puted dynamical form factor. We remind that since in this 
wavevector regime F{q,t) becomes independent of A'^, the 
Rouse model implies the scaling 



^^ = f(x) with x = q^rt 
Fiq,0) ^^' ^ 

being the scaling variable and 



qn 



Fig. 42 Reduced dynamical coherent form factor F(q,t)/F{q,0) as a 
function of reduced time x = q^Ft for several wavevectors q as indi- 
cated. The data has been obtained for BFM melts ((^ = 0.5, £ = 10, i" = 
0.0072) of chain length A' = 1024 with radius of gyra tion R ^ = 32.7. 
Main panel: Although the Rouse model scaling Eq. (147) holds for 
sufficiently small q, deviations are clearly visible for large waveve ctors 
q > 2n/S. The dashed line indicates the numerical solution of Eq. (149) 
for the intermediate wavevector regime. Inset: The double-logarithmic 
representation of ^(x) = 1 —F{q,0)/F{q,0) reveals for large wavevec- 
tors and small times (x <g 1) an exponent /3 = 3/4 (bold line). 



At short times the data for all s collapses on the same power- 
law slope with exponent ft) = 5/4 (bold line) as in Eq. ( 1140b . 
Confirming Eq. ( 112b and Eq. ( 1138b we obtain the scaling 



a^ T 



= ^fl^iy 



(147) 



(148) 



a rescaled monomer mobility |4|. Using the directly mea- 
sured mobility W we have F = 0.0072 for the BFM data 
presented in Fig. |42] As can be seen from the main panel, 
a satisfactory data collapse is indeed obtained if sufficiently 
large subchains are probed {q < 27r/8). This is consistent 
with the fact that our systems are to leading order of Rouse- 
type, especially if large times are probed. The scaling func- 
tion /(x) indicated by the dashed line has been obtained by 
integrating numerically the Rouse model prediction H 



f{x) = / daexp 



-a-jc'/^/z 



(«x-i/2y 



with 



CM) 



Cn(0 



^(wty 



-id+2)/4 



2 r 
forf<7;, (143) '<"'^ = ^l 



dx 



co^{xu) 



i.e. the colored forces acting on subchains add up indepen- 
dently to the colored forces acting on the total chain. 

6.3.4 Dynamic coherent form factor 

Neither the MSD nor the VCF of the (sub)chain CM can be 
measured directly in a real experiment. We discuss now how 
the exponent j3 =2 — ft) = 3/4 may be tested experimentally 
by means of the dynamical coherent form factor ||4]|5l 
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Considering that there is no adjustable parameter this fit is 
satisfactory. Expansion of Eq. ( 11491) for small x yieldo 



f{x) = \-x + x 



3/2 



Ay/2 



forx< 1. 
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F{q,t)=N{px{q,t)p,{q,Q)) 



N 



L (exp(/^-(r„(f)-r„(0))) 



(144) 
(145) 



n,m=l 



where p\{q,t) = L^=iexp(— /g • r„) stands for the Fourier 
transform of the density of a labeled chain. For t = the 
dynamical form factor F(^, f) reduces to the static form fac- 
tor F(^) discussed in Sec. 15.71 In the Guinier regime F{q,t) 
probes the overall chain motion which allows to determine 
quite generally /jn(0 using the expansion lU 



y 



1 



Fiq.t) 



1 — exp 



h^{t)q'^\ h^{t)q^ 



(146) 



F{q,Q) 

at least if sufficiently precise small angle data are available 
II78II79L Thus, y - f/^ for ? < t*. Obviously, Eq. (fT46] l also 
holds for our numerical data (not shown). 



Thus, according to the Rouse model one expects the rescaled 
dynamical form factor y{x) = 1 — F{q,t)/F(q,0) traced in 
the inset of Fig.|42]to increase (to leading order) as y (jc) ^ x" 
with j3 = 1 (dashed line). Instead we find for larger wavevec- 
tors a power-law increase with exponent j3 = 3/4 (bold line) 
in agreement with the scaling relation 

y^q\{t)^qV^ys^q^y' (151) 

for f ^ f* <C Ts [Eq. dMlb lPnln other words, an experimen- 
tal observation of a power law slope of the reduced dynami- 
cal structure factor with j3 w 3/4 would thus confirm scale- 
free CM displacement correlations with ft) = 2 — /3 w 5/4. 

'** Note that the approximation /(x) = exp{—2^/x/K) given in |4| 
is unfortunately wrong. Effective mobilities W determined using this 
formula are therefore also incorrect |52|. A better approximation tak- 
ing into account the fluctuations around the saddle point is given by 

/(x) Ri (;rV4)'/*exp(-2v^). 
'^ That with decreasing wavevector the data approaches systemati- 
cally the Rouse prediction /3 = 1 is expected (i) since for the same 
value X = q'^Ft a smaller wavevector q corresponds to a larger time t 
and (//) since for t 2> t* ~ s^N'^ the colored noise becomes masked by 
the white forces, Eq. ( 1141b . 
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Fig. 43 Robustness of the scaling of CN(f) 1901 . Data obtained for 
A' = 1024 using L26-moves is shown for (a) various overlap penal- 
ties e at volume fraction (j) = 0.5 and (b) various volume fractions ij) at 
overlap penalty e = 10. For sufficiently large times (but still f ^ Tn) 
all data appears to collapse on the same asymptotic power-law expo- 
nent a = 5/4 (bold line) for incompressible polymer solutions. The 
statistics deteriorates with decreasing e and p and additional physics 
is visible at short times. 



6.4 Compressibility effects 

Up to now we have focused on one working point at over- 
lap penalty £ = 10 and volume fraction (/) = 0.5, i.e. all data 
corresponds to the same static properties, especially to the 
same dimensionless compressibility g = 0.32. It is natural 
to ask how the observed scaling changes with the compress- 
ibility of the solution. Our key scaling relation Eq. (1133b cor- 
responding to an incompressible packing of blobs does in 
fact only depend implicitly on g. This suggest that Eq. (1140b 
for times f ^ 7n should also remain valid if one uses for 
rescaling of the axes the effective bond length b, the mobil- 
ity W and the diffusion constant Dn measured independently 
for the given operational parameters e and (j) . This is in fact 
borne out for the asymptotic behavior of the rescaled VCFs 
displayed in Fig. l43l a) for different overlap penalties e at 
(j) — 0.5 and in Fig. l43l b) for different volume fraction (j) 
at e = 10. (We only display chains with A^ = 1024 to avoid 
the colloidal regime for t ^ T^.) The values used for the 
rescaling of the axes are listed in Table [T] for the variation 
of the overlap penalty and in Table |2] for different volume 
fractions. The scaling collapse for long times presented in 
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Fig. 44 Dynamical correlations for 3D BFM melts at volume fraction 
<j> = 0.5 obtained using L06-moves: (a) VCF Cn(0 for e = 10 focusing 
on larger chains with A' > 256. Using W = 0.002 and Dn = 0.001 /N 
we rescale the data as in Fig. l40l a). This allows to collapse the data 
on t he sam e power-law slope (bold line) as obtained using L26-moves, 
Eq. {140). (b) Topology induced additional correlations for chains of 
length A' = 8192. Main panel: MSD /!n(0 vs. time t comparing L26- 
moves for e = 10 (spheres) to L06-moves for different overlap penal- 
ties e. Topological constraints become important for L06-moves with 
e 2> 10 and, as expected from reptation theory for times larger than the 
entanglement time Te ~ lO' 14], the data bends towards a power-law 
exponent 1/2 (dashed line). Inset: The reduced VCF —C-t^{t)N con- 
firms that for all system classes the shoit-time dynamics is described 
by the power-law exponent ft) = 5/4 (bold line). 



panel (b) demonstrates explicitly that Cn(0 '^ 1/b^p if one 
scales out the additional variation of the monomer mobility 
W{(j)). We remind that this density dependence stems orig- 
inally from the factor p^/p in Eq. ( 1133b due to the corre- 
lation hole forces setting the repulsion between chains and 
subchains which drive the dynamical correlations. That the 
data deviates for short time from Eq. ( 1140b is to be expected 
qualitatively since the incompressibility constraint is only 
felt by the chains if the dynamics is probed on a scale cor- 
responding to the static screening length. As we have seen 
in Sec. 15.41 a similar crossover is observed for static proper- 
ties such as the angular correlation function II86L At present 
we are still lacking a detailed description for the short time 
dynamical behavior and how to match it with the correlated 
motion of incompressible blobs at long times. 
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6.5 L06-moves and effects of topological constraints 

Up to now we have deliberately tuned our model to avoid 
topological constraints. Obviously, these constraints are ex- 
pected to matter for the dynamics of real 3D polymer melts 
I3JI13. It is thus of interest to see how the presented picture 
changes if topology is again switched on using the topology 
conserving L06-moves of the classical BFM |50"5T,'56l. 

Obviously, for small penalties the dynamics remains of 
Rouse-type as can be seen from Fig. l44T a) for the same static 
conditions ((/) = 0.5,e = 10) as in Sec. 16.31 One determines 
readily a local mobility W = 0.002 and a self-diffusion co- 
efficient Dn = 0.001 /N, i.e. the dynamics is slightly slower 
than for the larger L26-moves. If using these parameters 
Cn(^) is rescaled as shown in Fig. l40t a). all data sets collapse 
perfectly on the same slope as for L26-moves (bold line). 
This is consistent with the idea that the scaling Eq. (1140b 
only depends implicitly on the specific MC moves used. 

Preliminary data comparing L26-moves and L06-moves 
for chains of length A'^ = 8192 and higher overlap penalties 
is presented in Fig. l44T b). With increasing £ the crossing 
of the chains gets more improbable for L06-moves and the 
topological constraints become more relevant. This can be 
seen for £ = 30 and even more for £ = 100. The latter data 
set approaches the slope 1/2 (dashed line) expected from 
reptation theory for times larger than the entanglement time 
Tg |4|. The vertical arrow indicates a value for T^ obtained 
from an analysis of the monomer displacements h(t) |56|. 
Apart from the much larger chain length A'^ = 8192 used, the 
data for £ = 100 is consistent with the results obtained using 
the classical BFM (£ = oo) 151115611 . As in Fig. [32] the bold 
line represents the exponent j5 = 0.8 II51II79I . Superficially, 
it does a better job for systems with conserved topology 
due to the broad crossover to the entangled regime. How- 
ever, that this empirical exponent is by no means deep as 
revealed by the VCF Cn(0 plotted in the inset. For short 
times all systems with and without topology conservation 
are well-described by the same exponent ft) = 5/4 (bold line) 
in agreement with the proposed deviations from the Rouse 
model, Eq. ( 11341) . Note that the last points given for £ = 100 
decay slightly more rapidly. Unfortunately, the length of the 
analyzed time series does currently not allow us to verify 
whether our data become consistent with the decay 



'^~<""-^(l) ii 



-3/2 



for t > t; 



(152) 



1/2 

expected for reptating chains with d^ ^ T^' being the tube 
diameter |3|. Much longer time series, just as for the L26- 
moves with £ = 10 we have focused on, are currently under 
production to clarify this issue. The numerical demonstra- 
tion is challenging, since the difference between the expo- 
nents, 3/2 — 5/4 = 1/4, is rather small and several orders 
of magnitude in time are needed to discriminate the power 
laws. In any case it is thus due to the dynamical correlations 
first seen in the BFM simulations of Paul et al. II51II79I that 
the crossover between Rouse and reptation regimes becomes 
broader and more difficult to describe than suggested by the 



standard Rouse-reptation theory |[3]|3 which does not take 
into account the (static and dynamical) correlations imposed 
by the incompressibility constraint. 



6.6 Perturbation calculation predictions 

6. 6. 1 Colored collective forces 

Focusing now on the short-time correlations (f ^ Tn) in d- 
dimensional polymer melts without topological constraints, 
we present now a linear response calculation to describe the 
coupling between the degrees of freedom of a tagged test 
chain and the degrees of freedom of the collective bath char- 
acterized, respectively, by the dynamical form factor /^(^,?) 
and the dynamical structure factor S{q,t) at equilibrium. The 
underlying physics is that the displacement of a test chain at 
t = creates a free energy perturbation of the bath which 
decays by collective diffusion but survives at intermediate 
times. The collective force f (t) associated with the pertur- 
bation of the molecular field pushes the test chain towards its 
original position at / = 0~ causing thus the anomalous dif- 
fusion of the CM motion demonstrated numerically above. 
We remind first that the forces acting on the CM ^nIO of ^n 
a priori Rouse chain may be written 14] 



/vc 



dt 



(153) 



i.e. in addition to the random white force /(f) of the stan- 
dard Rouse model we have to account for the force / (f) 
from the molecular field. The friction coefficient ^ may be 
obtained using Eq. ( 1148b from the effective monomer mo- 



bility W. According to the relation Eq. ( 1176b given in Ap- 
pendix |A] the VCF Cn(0 of the chain CM is given by 



(154) 



CN(0 = -(^(/,(0-/.(0))forf>0, 

i.e. our task is to compute the correlation of the collective 
force /(f) due to the molecular field surrounding the ref- 
erence chain. The negative sign corresponds to the fact that 
the collective forces push the reference chain back towards 
its original position, just as for the overdamped colloids. 

6.6.2 Dynamical random phase approximation 



As already discussed in Sec. 16.3.41 the degrees of freedom 
of the test chain are encoded by the dynamical form factor 
F{q,t), Eq. ( 11451) . As seen from the main panel of Fig. l42l 
the scaling F{q,t) = F{q,0)f{x) with x = q^Ft is nicely 
obeyed for not too large wavevectors, although deviations 
from the Rouse prediction are visible for small x. Since these 
deviations are small, it is justified in the spirit of a perturba- 



tion calculation to assume that Eq. ( 11491) holds for all wave- 
vectors q ^ 1/^N and times f <C 7n. 

The degrees of freedom of the bath are described by the 
dynamical collective structure factor 

S{q,t) = {p{q,t)p{q,0))/{pV) (155) 
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Fig. 45 Scaling of the dynamical collective structure factor S{q,t) 
characterizing the degrees of freedom of the bath surrounding the test 
chain. Data obtained from BFM simulati ons (o pen symbols) are com- 
pared to the theoretical prediction, Eq. <158t . indicated by the bold 
line (9T). In agreement with theory we find a good scaling collapse if 
y{x) = S{q,t)F{q)/S^(q) is plotted as a function of the reduced time 
X = q* rt wit h r obtained independently from the monomer mobility 
W, Eq. J148t . Deviations from the Rouse model prediction are visible 
for small x, i.e. for short times f <g; f*, in agreement with Fig. [42] 



withp(^,f) =L"'=T^^P {^^1 'ln(O) being the Fourier trans- 
form of the collective monomer density. For f = Eq. ( 1155b 
reduces to the static structure factor S{q) analyzed in Sec. 14.41 
The dynamical structure factor cannot be devised using the 
Rouse theory alone, for S{q,t) encodes the collective be- 
havior of the chains, a cooperativity destroyed by the as- 
sumptions underlying the Rouse approach. Fortunately, to 
deal with this question the "dynamical Random Phase Ap- 
proximation" (dRR\) is a valuable tool which provides sen- 
sible predictions for S{q,t) |112 113,91 1. Similarly to the 
static RPA discussed in Sec. 14.41 the dRPA is based on a 
self-consistent closure of a mean-field theory: upon a per- 
turbation by an external field, the response of a polymer is 
assumed to be the combination of a direct response to the 
field plus a response mediated by the environment of other 
chains, which are also subjected to the external field. This 
autocoherent reasoning leads to a prediction for S{q,z), the 
Laplace transform of S{q,t), given by 19111 



zS{q,z)-S{q) = 



zF{q,z)-F(q) 



\+pvF{q)[F{q)-zP{q,z)\ 



(156) 



with V = \/pg being the effective monomer excluded vol- 
ume and F{q,z) the Laplace transform of the dynamical form 
factor F{q,t). For z = 0, i.e. by summing over all times 
t, this equation yields of course the static RPA S{q)^^ = 
F{q)~^ + vp discussed in Sec. 14.41 Using this limit one may 
rewrite Eq. (1156b which yields the dRPA formula 



F{q,z)/Fiq) 



s{q,z) ^ 

%) l+pvF{q)[l-zF{q,z)/F{q)] 



Focusing on incompressible systems (g<^l), large wavevec- 
tors (1//?N < ^ < 1/^) and short times t < Tn, Eq. ( 1157b 
can be further simplified. As shown in II91I 



S{q,t) S{q) f{z) 
— r-^ ~ — r-rS(x) with s(z) = ^ — 



(158) 



where f{z) and S{z) stand, respectively, for the Laplace trans- 
form of the scaling functions f{x) and s{x). Assuming for 
f{x) the Rouse model prediction, Eq. ( |149b , this yields for 
y{x) = S{q,t)F{q)/S^{q) the bold line presented in Fig.l45l 
Taking apart deviations for small x, i.e. for small times, we 
obtain a nice scaling collapse for the BFM data presented in 
Fig.|45] especially considering that there is no free fit param- 
eter available. We note finally that since for the wavevectors 
used S{q)Ki g= l/(vp), it follows from Eq. ( 11581) that 



{NpfS{q,t)F{q,t) 



{aq) 



J{x)s{x) 



(159) 



for sufficiently low wavevectors at fixed compressibility or 
for sufficiently small compressibilities at fixed wavevector. 

6.6.3 Lowest-order perturbation 

A collective density fluctuation 5p corresponds to a molec- 
ular field v5p driving the monomers of the test chain with 
a force —TvV5p. The lowest order force correlation cor- 
responds to the following mechanism: the test chain inter- 
acts with the collective field at time and again at time f. 
Using the densities p\{q,t) and p{q,t) of the test chain and 
the bath in Fourier space defined above, the collective force 
f (t) may be rewritten as 



r d"q 
L^t) = Tv j T^{iq)pl{q,t)p{q,t). 



(160) 



The correlation (f {t) ■ f (0) ) for ? > becomes thus after 
lowest-order factorization 



^'^ .2.2 



vV (Pl(^,0)pi(^,0) {p{q,0)piq,t)). 



(161) 



V .1 {2nY 

Using the dynamical form and structure factors, Eq. ( 11451) 
and Eq. ( 1155b . this yields 

f^it)-f^iO)) =T^Np I ^v^q^Fiq,t)Siq,t). (162) 

Using Eq. ( fT54] l and Eq. ( [159] ) the VCF becomes 



Cn(0 = 



d^'q 
Nl^^p J {Inf 



q^s{x)f{x). 



(163) 



Since f{x) and S{x) decrease rapidly for x ^ 1, the integral 
is well-behaved. After substituting the friction coefficients i^ 
and r by W, Eq. (11481) . this may be rewritten as 



(157) Cn(0 



~Cd- 



If'p 



{Wt)-'^Wi\h(0 = {d + 2)/A. 



(164) 
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Hence, Eq. (11641) confirms finally the scaling announced in 
Eq. ( 1133b and Eq. ( 1134b . The coefficient cj is given by 

c, = V^ {Adl^f^ I ^x^f{x)s{x) (165) 

which may be computed numerically [911. For 3D melts 
(d = 3) we find C3 w 0.42. This is of the same order but 
slightly lower as the empirical coefficient c w 1 estimated 
from Fig.|40] Where this small difference comes from is cur- 
rently still an unresolved question. One possibility is that the 
single loop approximation presented here — especially as- 
suming Eq. (1149b and Eq. ( 1158b for small x — is not suffi- 
cient and higher order terms should be considered. Another 
(physically more appealing) option is that only about half 
of the free energy is stored in the longitudinal composition 
fluctuations and that the transverse stress due to the shear- 
ing of the bath at constant density generates an additional 
contribution to the correlations of same magnitude Ii91ii93j . 



with dt being the (arbitrary) time window used to define 
the displacements. Since subchains repel each other due to 
the incompressibility constraint, a tagged subchain is pulled 
back to its original position by the subchain dipole field. 
Since for times dt -^t <^T^ the relevant dipole field decays 
much slower than for colloids (a = 1/2 — ?> 1/4), the correla- 
tions are much more pronounced. That our scaling relations 
for / <C 7n do not depend explicitly on the compressibility 
of the solution, as stated by Eq. ( 11331) , has been checked by 
the systematic variation of excluded volume penalty e and 
volume fraction (j) [90]. As shown for melts with topological 
constraints the early-time behavior, CN(f) ^ —N~^t~^''^, is 
preserved before entanglement effects set in which leads to 
a broad crossover between Rouse and reptation regimes. As 
shown in Sec. 16.61 our scaling approach is consistent with 
a standard linear response calculation. En passant we have 
verified the so-called "dynamical Random Phase Approxi- 
mation" for the collective dynamical response of the bath 
surrounding the reference chain (Fig.|45]) |!9T1. 



6.7 Sunmiary 

Focusing on a variant of the BFM algorithm with topology 
non-conserving moves we investigated the effective dynam- 
ical forces acting on the chains and subchains imposed by 
the incompressibility constraint. Sampling chain lengths up 
to A'^ = 8 1 92 allowed us to carefully check the A'^-dependence 
of deviations with respect to the Rouse model. Such devia- 
tions are visible from the short-time scaling of the CM MSD 
h^{t) in agreement with the literature, however, a more pre- 
cise characterization can be achieved by means of the CM 
VCF Cn(0 w d^hn{t)/2 which allows to probe directly the 
correlated random forces, Eq. ( 1154b . How such a VCF can 
be computed within a MC scheme has been first illustrated 
for dense BFM beads (Fig. |36] | confirming the negative al- 
gebraic decay expected for overdamped colloids, Eq. ( 11321) . 
The observed exponent (0 = {d + 2)a with a = 1/2 can be 
understood by the coupling of a tagged colloid to the gra- 
dient of the collective density field decaying in time. As 
shown in Fig. l40l b). the same exponents are relevant for 
polymer chains for large times (/ ^ 7n) where the chains 
behave as effective colloids. More importantly, the short- 
time deviations for /7N(f) have been traced back to the neg- 
ative analytic decay Cn(0 "^ —N'^t'"" for t <^ 7n with an 
exponent O) w 5/4. That Cn{t) decays inversely with mass 
shows that the process is local, i.e. the displacement corre- 
lations of subchains of arc-length i add up independently 
(Fig. 1411 ). Assuming according to the postulated scaling re- 
lation Eq. (1133b the chain relaxation time Tn to be the only 
characteristic time scale, both asymptotic regimes can be 
brought to a successful data collapse. The a)-exponents for 
short times proposed in Eq. ( 11341) are consistent with the 
crossover scaling, Eq. ( 11331) , and the locality of the corre- 
lations (Cn(0 ^ 1/AO which implies CO = {d + 2)/4. This 
scaling can be understood by the generalization of the above- 
mentioned correlation experienced by overdamped colloids 
to the displacement field of subchains of length .s ^ 5f '/^ 



7 Conclusion 

7.1 Summary 

Until very recently it has been generally assumed that all 
long-range static and dynamical correlations are negligible 
in dense 3D solutions of flexible homopolymers beyond the 
excluded volume screening length ^ ^^ g^'^ which character- 
izes the decay of the density fluctuations |j2||3jl4j,8J. For static 
properties this general screening assumption leads to Flory's 
ideality hypothesis (Sec. 11.3b |2 4| stating that the chains 
must obey Gaussian chain statistics. For the equilibrium dy- 
namics it implies Rouse model dynamics (Sec. 11.41 ) if in ad- 
dition momentum conservation ("hydrodynamic screening") 
and topological constraints [4| may be neglected or are de- 
liberately switched off as one can readily do in a computer 
experiment using local MC hopping moves (Sec. l3.3.Tb . 

That some deviations from Flory's ideality hypothesis 
must exist has been suggested in earlier work by de Gennes 
ll3lfTT4l. Obukhov L69J, Duplantier UJ8J, Schafer et al. \J} 
ItTII or Semenov and Johner ll94l . For instance, it has been 
known for a long time (3] that due to the incompressibility 
constraint an entropic penalty m^ w Pn/P ^^t by the correla- 
tion hole self-density p^ w ^/^n has to be paid for bringing 
two chains of length N close to each other (Fig. [3]). In the 
large-A' limit this effect must become negligible in d = 3, 
however. That the dynamics of polymer melts without rel- 
evant topological interactions may deviate from the Rouse 
model has been pointed out in various theoretical, computa- 
tional and experimental studies by Schweizer 1 175], Guenza 
et al. 1 99 1 or Paul et al. I5W7\ as summarized in Ref. 179|. 
However, since topological and hydrodynamical effects have 
not been systematically separated from the incompressibil- 
ity constraint |99 1 — with the notable exception of the BFM 
simulations of Shaffer II54II55I — it has been difficult to pin- 
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point the precise physical origin of the key finding observed 
for the MSD h^{t) of the chain CM, Eq. (O. 

Summarizing recent theoretical and computational re- 
sults we have argued in the present contribution that due 
to the incompressibility constraint both static 180,81,821831 
[84ll85ll86.87 .88 891 and dynamical ||90ll9Tll correlations of 
the composition fluctuations, i.e. of the density and displace- 
ment fields of marked chains and subchains, must arise for 
distances r ^ ^ and for corresponding arc-lengths s ^ g. 
Since the applied constraint is scale-free, the resulting cor- 
relations are also scale-free — taken apart the usual upper 
and lower cutoffs set by the chain and the monomer size 
— and do not depend explicitly on the dimensionless com- 
pressibility g of the solution (Sec. 12.31) . The central point 
stressed by us [811 is that not only chains but also subchains 
of arbitrary length i must repel each other due to the incom- 
pressibility constraint and this with a penalty m* w s/pR^ 
being for s <C A'^ much larger than the penalty m^ between 
the chains (Fig.lHl. Since subchains repel each other and this 
with a strength decreasing with sind = 3, the chains become 
weakly swollen as discussed in Sec.|5] Our study shows that 
a polymer in dense solutions should not be viewed as one 
soft sphere (or ellipsoid) [123 124 99 1, but as a hierarchy of 
nested segmental correlation holes of all sizes aligned and 
correlated along the chain backbone rFig.[3lb)lFlThe effec- 
tive interaction between subchains has also dynamical con- 
sequences since a subchain displaced at f = causing thus 
a perturbation of the subchain density field (Fig. [35l l will 
be slightly pushed back by the bath to its original position 
(Sec.ili. 

In this review we have focused on numerical results ob- 
tained by means of a BFM variant with finite monomer in- 
teractions without topological constraints (Fig.|6]l. Since this 
"soft BFM" is fully ergodic (in contrast to the classical BFM) 
and very efficient due to its implementation as a Potts spin 
model 1861 . it may be an interesting alternative to various 
popular coarse-grained simulation approaches with self-con- 
sistently calculated effective pair interactions Ill23lll24ll99l 
[180 166 1 . Our study has not been limited to monodisperse 
polymers but we have also investigated (essentially) Flory- 
distributed EP 181 1183 |[87l [881. Since in these systems chains 
break and recombine constantly (Fig. [9]l equilibration and 
sampling become much faster than for their monodisperse 
counterparts. This is computationally of interest since the 
numerical demonstration of the various theoretical predic- 
tions requires the sampling of large chain lengths to avoid 
additional chain end effects. 

The deviations from the screening assumption are in- 
deed scale-free as made manifest by the numerical obser- 
vation of the analytic decay of various correlation functions. 
As announced by Eq. (Illb . it has been shown in Sec.[5]that 
the bond-bond correlation function Pi{s) ^ d^R^ decays as 
P\{.s) ^ +\/s'^l'^ for s<€.N 7mA d = 3. That these correla- 
tions are long-ranged is made explicit by the power-law de- 
cay of the bond-bond correlation function P\{r) ^ +l/r^ 

*" Similar deviations from Flory's ideality h ypothesis have been re- 
ported recently for polymer gels and networks I139I179I . 
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Fig. 46 Polymer melts confined to ultrathin films of width H = 4 
for volume fraction ^ = 0.5, projected 2D density p2 = Hp = 0.25, 
overlap penalty e = 10 and several chain lengths A' as indicated. In- 
set: Bond-bond correlation function Pi{s) compared to the power- 
law exponent a = d/2 = 1 predicted by Eq. (TT). Main panel: VCF 
C-f^{t) using a similar representation as in Fig. l40l"a) . The exponent 
CO = {d + 2)/2 = 2 expected for large times, Eq. l ll32t . is represented 
by the dashed line, the short-time exponent O) = {d + 2)/4- = I pre- 
dicted by Eq. n33t by the bold line. A perfect data collapse is observed 
for?<s:rN, i.e. Cn(0~1/A'- 



as a function of the monomer distance r (Sec. 15.4.21) . As a 
second key result of this study, the displacement correlation 
function Cj^{t) w d^hj^{t) is shown to reveal for short times 
f <C 7n a negative algebraic decay according to NC-^{t) ^ 
— I/r''+^'''* as predicted by scaling arguments and pertur- 
bation theory (Sec. [6]). As stated by Eq. ( II2I) . this analytic 
decay holds more generally for the displacement correla- 
tion function Cs(f) of subchains of length s <N and times 
t <^TsKi s^/W, Fig. [41] which shows that the correlations of 
different subchains add up independently to the correlations 
of the total chain ("locality"). 



Note that earlier computational studies have focused on 
(second and higher) moments of the generalized displace- 
ment field under consideration such as the mean-squared 
size /?s of subchains of arc-length 5 |[73l or the CM MSD 
hj^it) L51 1 and have thus only probed indirectly the respec- 
tive colored forces, Eq. (1176b , corresponding to the devia- 
tions from the general screening assumptionIf3 Computing 
numerically correlation functions such as Pi{s) or Cn(0> 
rather than twice their integral, allows to probe directly the 
colored forces, Eq. ( 1176b . without having to subtract first the 
white noise and the local physics which contribute both to R^ 
and /in(0- ^s a consequence, this allows to demonstrate that 
the deviations from the screening assumption are not due to 
(non-universal) physics at the lower cutoff but, in fact, are 
present for all arc-length s (Fig.[25]i and all times t (Fig.[40li. 
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Strictly ID and 2D self-avoiding polymer melts (e = oo) 
are known to become compact and segregated with g* w 1 
and are thus not Gaussian beyond the monomer scale. While 
the ID case is trivial, it should be noted that for strictly 2D 
melts compactness means V = l/d = 1/2 |3|. This does 
of course not imply ideal chain behavior since other crit- 
ical exponents characterizing the chain conformations dif- 
fer from those of non-interacting Gaussian chains as shown 
in the pioneering work by Duplanier 1 1781 . For instance, 
since 02 = 3/4 II178I this implies that the compact 2D chains 
have a fractal perimeter of dimension d^ = d — 02 =5/4 
[94 1 . This theoretical prediction has been verified numeri- 
cally II119II120II118I from the scaling of the chain perimeter 
(determined by "box counting" [191) and of the form factor 
shown to be consistent with Eq. (1166b . 

Self-avoiding melts in strictly one and two dimensions 
correspond to rather specific universality classes II94II192I 
|193ll96l . Since systems of strictly one monomer layer thick- 
ness at high volume fractions remain experimentally a chal- 

lenge 1 184,188], it is of interest to relax the non-crossing 

I178ll94ll96lll91lll92lll93l . Especially, the perturbations to chainconstraint by either allowing some finite monomer overlap 



7.2 Related questions and outlook 

7.2.1 Confined polymer melts 

In the presented numerical work we have focused on 3D 
polymer bulks (d = 3). It has been shown that even the devi- 
ations with respect to the general screening assumption are 
well-described by means of one-loop perturbation calcula- 
tions (Sec. |2l Sec. 16.6b and this for arbitrarily long chains 
and subchains. As already noted in Sec. 11.11 the success of 
the mean-field approach is expected due to the scaling of 
the (sub)chain self-density p^ /p ^ s^'"^^, i.e. the number 
of subchains a reference subchain interacts with increases 
as p/p* ^ y/s ^ 1 in (i = 3. Due to the increasing experi- 
mental interest on mechanical and rheological properties of 
nanoscale systems in general HIS II and on polymer melts 
confined to surfaces, thin slits {d — 2) or capillaries {d = 1) 
in particular lll82lll83llT84lfT85lfT86lfT87lll88lll89llT90l one 
is naturally led to question theoretically such mean-field cal- 
culations for systems of reduced (effective) dimension 1 1 14 



dynamics due to geometric constraints remain a challenge of 
significant technological relevance with opportunities rang- 
ing from tribology to biology lfT82l[183][TB5][T87l . 

Note that for ideal chain systems with d < 1/v = 2 the 
self-density p* ultimately exceeds the density p at a charac- 
teristic chain length g* which depends on the specific prob- 
lem considered. For (sub)chains larger than this length g* 
one expects the chains to adopt compact and segregated con- 
formations, i.e. Flory's exponent V is given by the spatial di- 
mension, V = l/d \19\. Since these compact structures arise 
ultimately due to the scale-free incompressibility constraint, 
their surfaces are expected not to be ruled by a finite sur- 
face tension (which would imply a length scale) II94II118I 
but to be also scale-free and described by a fractal surface 
exponent d^ < d |il9j| . In most cases of interest the surface 
exponent becomes d^ 



d— 92 11181 . i.e. it is set by the 
well-known contact exponent 02 describing the scaling of 
the distribution G{r,s) ^ r^^ for short distances r <^ R^^ be- 
tween two monomers n and m = n+s ona chain. (Obviously, 
02 = for Gaussian chains, Eq. (|2]l.) According to the gen- 
eralized Porod law the intramolecular structure factor F [q) 
of such compact objects scales as 



F{q)/N^\/{Rnq) 



2d-ds 



(166) 



in the intermediate range of the wavevector q f5"T94 1 1 81 . If 
ds <2d ~ 2, i.e. d2>2 — d, the standard Kratky representa- 
tion reveals thus a much more pronounced no n-mono tonous 
behavior as for the 3D melts discussed in Sec. l5.7Pn 

*' In line with Appendix [A] we view here the bond vectors /, of a 
chain conformation as displacement vectors with the monomer index / 
playing the role of time. 

*^ Incidentally, melts of non-concatenated rings in rf = 3 have been 
argued to become "marginally compact" with v = l/d = 1/3 and 
ds ^ d = 3, i.e. the form factor is predicted to scale as F{q)q^ ~ 1/? 
[118 195] in the large-TV limit. That the rings should adopt compact 
configurations is expec ted due to the mutual repulsion caused by the 
topological constraints I132I121I133I168I196I197I198I122I199[T951 



penalty e (Sec. 13.4b or by considering chains confined to 
capillaries or slits between parallel walls of finite distance H 
as in Fig.|6tc). Note that g* w {ciH^pY ^ H'^ for melts con- 
fined to a ID capillary of finite width H as already shown by 
Brochard and de Gennes in 1979 



111 141 and g* is known to 
increase even exponentially with H for essentially 2D ultra- 
thin films 194] . Since g* becomes thus rapidly large, pertur- 
bation results, such as the scaling of the bond-bond correla- 
tion function Pi (s) stated in Eq. (Illb . must become relevant 
for finite (sub)chain lengths. A lower limit of validity g^, of 
the (i-dimensional perturbation calculation is set for systems 
with finite overlap penalty by g* w g and (sub)chains con- 
fined to capillaries or slits of finite width should be larger 
than g^ w H^, otherwise the chains do not "feel" the wall 
constraint and behave as in the 3D bulk. 

Let us focus now on effectively 2D films of width H 
ri36:T37l l95ll97l as shown in Fig.gH Note that all properties 
considered here are the 2D projection of the 3D observables. 
The relevant density is, e.g., the projected number density 
p2 = Hp = 0.25 and the indicated bond length b = 2.65 is 
the projection of the effective bond length of the 3D melt. 
Since the width H = 4 used here is much smaller than the 
typical subchain size Rg considered, i.e. 5 ^ g* w 1, these 
systems can be regarded as effectively d = 2^ dimensional 
for all s and A'. (We write c? = 2+ to stress that monomer 
overlap and chain crossings are allowed.) Since also g* ^ s 
one expects the perturbation result for the bond-bond cor- 
relation function Pi{s) ^ l/i™ with CO = d/2 = 1 to yield 

200 118|. Since there is no obvious reason for a finite surface tension, 
the surface must become fractal which may be determined both in ex- 
periment as in a computer simulation by the generalized Porod law 
of the form factor. A marginally compact structure allows to keep all 
monomers evenly exposed to the topological constraints imposed by 
other rings, i.e. all subsegments of the rings are thus ruled in a self- 
similar manner by the same statistics. Such behavior is known for var- 
ious biological systems, such as the lungs of mammals, atte mptin g to 
maximize the surface at constant overall embedding volume 12011 . 
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a reasonable fit to the measured data for all s <^ N. This is 
indeed borne out nicely for the data presented in the inset. 
According to Eq. (|4]l this implies that R^/s and R^/N must 
increase logarithmically with, respectively, i and A^. This can 
be checked directly as seen in Re f. 1 95 1 . The same paper also 
investigates the form factor F{q) for different H. It is shown 
that q^F{q) becomes systematically more non-monotonous 
with decreasing width H, i.e. the more p'^/ p increases. 

We turn now to the dynamics of BFM melts confined 
to thin slits of width H = A sampled using local topology 
non-conserving L26-moves as in Sec. 16.31 197|. It can be 
demonstrated that to leading order the dynamics remains 
consistent with the Rouse model, especially Eq. ©Ifj As 
before, we obtain from the monomer MSD h{t) a local mo- 
bility, W w 0.002. That the self-diffusion coefficient scales 
asDN RiO.OOS/A'^canbeseenby \xdicmgNh^{t) /2dt in log- 
linear coordinates. (Note that W and D^ are consistent with 
Eq. (1137b .) Systematic deviations are again revealed by the 
short-time behavior of the MSDs h-^{t) and /js(f ) of chains 
and subchains (not shown). These deviations are, however, 
much more pronounced than in 3D and are visible over time 
scales up to the total chain relation time T^. A more precise 
characterization of the colored forces acting on the chains is 
again achieved by means of the VCF Cn(0 which is repre- 
sented in the main panel of Fig. |46] The representation cho- 
sen is similar to Fig. HOla) for the bulk case; the factor b^ p 
being replaced by b^p2- The faster decay observed for the 
shorter chains at f ^ 7n is described by the exponent (0 = 
{2 + d)/4 = 2 (dashed line) expected for 2D overdamped 
colloids, Eq. (1132b . As may be seen from the collapse for 
short times (f ^ Tn), the VCF scales again as Cn(0 '^ 1/A^, 
i.e. the colored forces are again localo In agreement with 
Eq. ( 1134b or Eq. (1142b the VCF is observed to decay alge- 
braically with an exponent (0 = {d + 2)/4 = 1 (bold line). 
Since Cn(0 ~ d^h^{t)/2 it follows that h^{t)/t is not con- 
stant but must increase logarithmically for all times t ^ 7^. 
In analogy to Eq. (1141b for the bulk case this thus explains 
the observed deviations from the Rouse model. 

We note finally that qualitatively similar results consis- 
tent with Eq. ( fTTT) and Eq. ( 1134b are also found in prelim- 
inary BFM simulations of polymer melts confined to thin 
capillaries with a square-section of width 7/ = 8 ^ /?n and 
a broad range of overlap penalties e. 



7.2.2 Interchain correlations and anti-Casimir forces 

The presented work has focused on intrachain properties 
such as the bond-bond correlation function Pi {s). The main 
reason for this is that the power-law exponents associated 

^^ Since in thin films the ratio R^/N diverges logarithmically with A' 
1 94 95 154] and since Tn r^ R^/Df^ still holds, there are some (rather 
weak) logarithmic corrections to Eq. ^ as discussed in |97|. 

*"* Note that a different scaling, Cn(?) ~ l/x/ZV, has been reported in 
a study of strictly 2D polymer melts sampled by MD simulation us- 
ing a standard Langevin thermostat 1 177|. The coiTclations of the CM 
motion are argued to become non-local in this special limit due to the 
"constant surface constraint" of the compact and segregated subchains. 



(a) 1 external field 




; 5 ; 





(b) 




EP bath: \i 



i, grand-canonical 
MC moves 



density ? energy ? 



Fig. 47 The presented studies focused on intrachain properties which 
are numerically more readily accessible than the long-range interchain 
coiTelations associated to the "anti-Casimir" forces 1 18 1 15 1. The two 
panels indicate two possibilities to probe such correlations for EP melts 
confined to thin slits: (a) A weak external field conjugated to the bond 
direction is applied to the bonds in a thin slap of width 5: The sec- 
ond Legendre polynomial Piir) characterizing the orientation of the 
bonds with respect to the direction perpendicular to the slap and the 
force f[r) experienced by a colloid placed at a distance r should decay 
analytically with not too large exponents, (b) EP in grand-canonical 
contact with a large bath: The difference of, e.g., the densities of the 
confined system and the bath decays analytically with H. 



to the intrachain static and dynamical deviations due to the 
incompressibility constraint are not too large being thus nu- 
merically still accessible with reasonable computational ef- 
fort. In fact, similar correlations have been predicted the- 
oretically also for interchain properties 118111151 which for 
macroscopic thermodynamic and mechanical properties may 
even be more relevant. Note that for these interchain corre- 
lations it is not only the incompressibility constraint which 
matters but in addition the constraint that closed loops are 
disallowed and must be eliminated from an extended grand- 
canonical ensemble containing both linear chains and rings, 
as shown by Obukhov and Semenov 118111151 . The "throw- 
ing-away" of configurations is argued to generate additional 
entropic forces repelling, e.g., two large colloids immersed 
in a linear-chain polymer melt. Unfortunately, these "anti- 
Casimir forces" correspond to such strong exponents as a 
function of distance r or wavevector q (similar to the stan- 
dard van der Waals forces |20|) that, at present, it has turned 
out to be elusive to verify them numerically although a brave 
attempt has been made by means of off-lattice MC simula- 
tions 1 202 1 . Following this recent work, we sketch in Fig.l47] 
two possibilities which might allow to probe correlations 
due to the no-loop constraint by means of EP melts confined 
to thin slits. The use of EP instead of monodisperse chains 
should speed up the sampling of independent configurations 
while confining the monomers reduces the exponents of the 
various analytic scaling relations predicted. For instance, the 
second Legendre polynomial Piir), measuring the orienta- 
tion of the bonds at a distance r from a small region where an 
external field weakly aligns the bond vectors, should decay 
as Pzix) ^ r~^ with 5 = 2 for t/ = 2+ as illustrated in panel 
(a) of Fig. |47] If successful these simulations may stimulate 
real micro-mechanical experiments using optical tweezers 
aligning bond vectors and repelling thus nearby colloids. 
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7.23 Dynamical properties beyond MC stochastics 

We have checked by means of MD simulations of a bead- 
spring model II91II92I that Cn(0 « -TV'f"^/'* for t < 7n 
holds as long as a large friction constant y is used for the 
Langevin thermostat applied |24|. Concerning the power- 
law amplitude it should be noted, however, that the coef- 
ficient c determined from the data using Eq. ( 11401) is much 
larger for these MD simulations (c/c3 w 20) compared to our 
BFM result {c/c^ w 2). As we have pointed out in Sec. 16.61 
small differences of order one to the predicted value C3 = 
0.42 are expected due to the approximations made in the 
one-loop perturbation theory, especially since the forces as- 
sociated to the shear stress may not be properly accounted 
for 1 91]. However, the much larger ratio c/ci, found for our 
MD simulations points to additional physics related to the 
screening of the monomer momentum in these strongly in- 
teracting viscoelastic systems [92,93 1 13 Since our BFM and 
MD simulations correspond to different model systems, it 
is currently not possible to clarify this issue. In the future 
we plan to address this problem using a generic (off-lattice) 
bead-spring Hamiltonian with soft beads and large springs 
such as, e.g., the model used by Spenley |68| since we re- 
quire a model which is (i) essentially incompressible {g ^ 
1), (ii) without having too much useless local structure (as 
for the Lennard-Jones beads used in |91 ,92,93]), (iii) does 
not conserve topology and (iv) can be used for different dy- 
namical methods such that we can compare the melt dynam- 
ics under the same static conditions. For such a model one 
should then compare the dynamics obtained by means of 



• MD simulation with a strong Langevin thermostat 12411 : 

• local MC moves of the monomers II26II : 

• BD simulation in the sense of a "position Langevin equa- 
tion" where the momentum variables are dropped from 
the equations of motion II22II203L 

The BD method has the advantage that one needs not to in- 
ject strong random forces as for MD to get rid of the mo- 
mentum to get closer to the overdamped motion assumption 
implicit to theory 191]. We expect to obtain for the MC and 
BD simulations similar ratios c/cj, w 2 as for our BFM study. 
Interestingly, much stronger dynamical correlations are 
revealed f92"| in MD simulations using a weak Langevin 
thermostat (y — > 0) or a momentum conserving "Dissipa- 
tive Particle Dynamics" (DPD) thermostat 120411241 . As ex- 
pected from the hydrodynamic screening assumption |f4l|, the 
chains are described to leading order by the Rouse model, 
i.e. Eqs. (1136b , ( 1137b or (1149b hold. However, it is shown 
that some long-range correlations arise due to an intricate 
coupling of hydrodynamics (momentum conservation) and 
the (transverse) viscoelastic response of the bath following 
the displacement of a reference chainj^ As a consequence, 

^^ The beads bounce into each other on time scales much shorter than 
1/7. Kicks due to conservative or random forces can thus propagate 
through the system before the thermostat has time to do its job. 

^* Note the analogy to the static excluded volume screening dis- 
cussed in Sec. [2] Although the chains are essentially described md = 3 
by Gaussian chain statistics, not all long-range correlations vanish. 



even the VCF Cs {t ) of subchains for t <^Ts depends now ex- 
plicitly on the total chain length A'^, i.e. the correlations are 
non-local with Cn(0 ^ 1/v^ as suggested by the Zimm 
model. Since the scaling is still described by Eq. ( 11331) this 
implies an exponent ft) = 3/2 for short times. Hence, 



CN(f) ^ 



I) f <^»/')' 



1 



(167) 



which is the central formula real experimental data should be 
compared with. Although the VCF decays now more rapidly 
in time, it is stronger than the VCF for perfectly overdamped 
melts, since the A'^-dependence in Eq. ( 1167b is weaker. Note 
that at f w 7n both mechanisms become of the same magni- 
tude: Cn(7n) ^ \/N''I'^. Using the same bead-spring model 
mentioned above it would now be of high interest to sample 
the dynamical response using 

• a DPD thermostat following the work by Spenley 1681 . 

Obviously, these issues are clearly outside the realm of stan- 
dard MC based stochastics which shows their limitations. 
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A Displacement correlations 

Let x{n) denote a (real valued) stochastic variable in some vector space 
with n being a discrete or continuous index characterizing this variable. 
Examples are the position r, of the monomers of a given chain as a 
function of the monomer index i (Fig. lir'l or the center-of-mass (CM) 
position £n(0 of "^ chain of length A' as a function of a time t measured 
by the discrete number of Monte Carlo Steps (MCS). We assume trans- 
lational invariance of ensemble averaged properties with respect to n, 
i.e. properties such as the "mean-square displacement" (MSD) 



h{n) = h{ni = tio+n.no) = ((x(«i) -x(no))^) 



(168) 



only depend on the difference n = iii — no > and not on the abso- 
lute indices ng or iii . Taking advantage of the translational invariance, 
the ensemble average (...) is often sampled over all available pairs 
of indices (;j| = no -|-n,«o) for a given n. Examples for MSDs ob- 
tained as "glid ing averages" are the typical segment size Rl presented 
in Sec. [5^ or the MSD h{t) of BFM beads (A^ = 1) discussed in 
Sec. 16.21 In this review we consider essentially Fickian stochastic pro- 
cesses where the associated generalized MSDs are to leading order 
proportional to time n, Eq. (O or Eq. ^. Since the deviations with 
respect to the Fickian behavior are small, it is important to sample di- 
rectly correlation functions corresponding to higher derivatives of h{n) 
with respect to time. 

Since we have to deal with indices n which may be either discrete 
or continuous, let us introduce the (forward) difference operator A„ 
[44 1 defined by A„f(n) = f{n + 5n) — f(n) with /(«) being a function 

*^ This monomer index ; may be eith er a discr ete (as in our study) or 
a continuous variable as, e.g., in Refs. 1 1801 166 1. 
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of our vector space and 8n an arbitrary (typically) small shift of n. Let 
us define the "velocity vector" 



8n 



5n 



x{n + 8n) —x{n) 
8n 



(169) 



associated to x{n) and the "velocity correlation function" (VCF) 

C{n) =C{n\ =110 + 11, no) = (v(«i =no + n) •v(/Jo)) (170) 

which measures the correlations of two velocities v(«o) and v(ni ). C{n) 
automatically vanishes if both vectors are uncorrelated. Note that C{n) 
is a four-point correlation function with respect to n depending in gen- 
eral on the shift 5n of the index. If 5n > n both displacements become 
trivially correlated and we have lim„^oC(«) = h{Sn)/5n^ > 0. To see 
how h{t) and C{n) ar e rela ted for n 2> 5n let us apply twice the differ- 
ence operator to Eq. <168t . This yields 

A„ A,„ {ixini)-x{no)f) = -2{v{nt)-v{no)} Sn" (171) 

where we have used that the difference operator and the averaging pro- 
cedure commute. Since the MSD h(t) can be assumed to be mathemat- 
ically well behaved and using the translational invariance it follows on 
the other hand side that 



An^A„^h{n 



I — no) ~ —^n^{") S"^ for " 2> 5n. 
Altogether this demonstrates that 



C{n) » -d„h{n) for n:s> 8n 



(172) 



(173) 



in agreement with Eq. J4) for the bond-bond correlation function Pi {s) 
and with Eq. {7} for the chain velocity correlation function. Note that 
since h{n) does not depend on 8n, C{n) does not depend on Sn either 
in the limit n 2> 5n. 

The stochastic variables of interest in this review are described by 
a position Langevin equation of form 1 116j 



? 



dx{n) 
dn 



■i;v{n)=Mn) = Mn)+Mn) 



(174) 



with ^ being the "friction constant" and f,{n} the total "force" acting 
on x{n) which may be decomposed in a random white force contribu- 
tion /^ («) and a colored force ./c (" ) stemming from the remaining (non- 
white) interactions of the degree of freedom under consideration with 
the bath and its constraints. Perfectly Fickian behavior would be ob- 
tained if only the white force contribution were present (fc{n) = 0)F1 
For x{n) being the position of a monomer or the CM of a chain as a 
function of time t, Eq . <174t corresponds to an ov erdamped motion as 
discussed in Sec. l6.2lf or MC bea ds or i n Sec. l6.3l for Rouse-like chains 
in the melt, Eq. l ll53t . From Eq. l ll74t it follows for the displacement 
correlation function that 



C(«) = (v(n).v(0)) = r'(/.W-/.(0)). 

The total force correlation function {ft(n) 
now into four contributions: 



(175) 
• /t(0)) at n > decomposes 



(a) (/r(«) • /r(0)) =0 since /^(n) is a white force, 

(b) (/r(«) • /c (0)) =0 since future white forces cannot be anticipated, 

(c) (/c(«)-/c(0))and 

(d) {fc{n) ■ f,{0)} = -2(/c(«) •/c(0)) due to the odd «-symmetry of 
the velocity v{n) = -v{-n) t255i FI 

^* While fr{n) and fc(n) are both of zero mean, one may add in addi- 
tion an external force /e which allows to determine the friction constant 
from the drift velocity: ^ (v(n)) = f^. 

^' Since v{n) = — v(— n) this implies ft{n) = —ft{—n) for the to- 
tal force which in turn yields (Mn) ■ fdO)} = - (/t(-«) •/c(0)). 
The identity (d) is then obtained by substituting {ft{n) ■ fc{0}} = 
(/cW-/c(0)) (due to (b)) and (/,(-«) ./,(0)) = (/cW-/t(0)) = 
(/cW-./i-(0)) + (/cW-/c(0)). 



Summing up over all contributions and using Eq. n75t it follows 
C(«) = -r'(/c(n)-/c(0)) forn>0, (176) 

i.e. the displacement correlation function probes directly the correla- 
tions of the colored force fc{n). This general result is used in Sec. 16.61 
for the motion of the chain CM with time t but does also apply for the 
monomer position r, as a function of the monomer index / where the 
bond-bond correlation function Pj {s) measures the colored forces act- 
ing on the chain. It shoul d be stressed, however, that in this case the 
forces noted in Eq. {174) do not correspond directly to the standard 
forces acting on the monomers in real time. 



B Static properties 

B.l Moments and generating function 

Higher moments of the segmental size distribution G{r,s) can be sys- 
tematically obtained from its Fourier transform G(q,s) = ^[G(r,s)] 
which is in this context sometimes called the "generating function" 
|fTT6|. For ideal Gaussian chains the generating function is given by 
Go(g,i) = exp{—s{aq)^) where we have used a^ = b^ /Id to simplify 
the notation. Moments of the size distribution are given by proper 
derivatives of G{q,s) taken at g = 0. For example. 



(177) 



{r^")={-lYA''G{q,s)l^^o 

with A being the Laplace operator with respect to the wavevector q. A 
moment of order 2p is, hence, linked to only one coefficient A2p in the 
systematic expansion, G(q,s) = Y^„=o^2p'f''' , of G{q,s) around q = Q. 
For our example this implies 

{r^>') = (-l'n2p+l)\A2p (178) 

in general and more specifically for a Gaussian distribution (zi^'')n = 
' ''\ ' si'a^i'. The non-Gaussian parameters read, hence. 



6Pp\ {r}P) 
{2p+l)\JPY 
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(179) 



which implies (by construction) a^ = for a Gaussian distribution. As 
various moments of the same global order 2p are linked to the same Ai/, 
they differ by a multiplicative constant independent of the details of 
the (isotropic) distribution G{q,s). For example, (r^) = 6|A2|, (r*) = 

I2OA4, (x^) = (y^) = 2IA2I, (xV) = 8A4 with X and y de noting the 
spatial components of the segment vector r. Using Eq. ( |179t for p = 2 
it follows that 



Kx\ 



(x2)(y2) 



Aa , , 

= l-2-^ = a2{s) 



(180) 



i.e. the properties a2{s) and K^yis) discussed in Figs.|29landl30lmust 
be identical in general provided that G{q,s) is isotropic and can be 
expanded in q^. 



B.2 Deviations of the segmental size distribution 

We turn now to specific properties of G{q,s) computed for formally 
infinite polymer chains in the melt. In practice, these results are also 
relevant for small segments in large chains. A' 2> .s 2> 1, and, espe- 
cially, for segments located far from the chain ends. These chains are 
nearly Gaussian and the generating function can be written asG(q,s) = 
Go(q,s) + 5G(q,s) where according to Eq. l[56t we have a small peitur- 
bation SG{q,s) = — {UG}o + {U}o (G)q due to the effective interaction 
potential \{q) given by Eq. {66). In this paragraph we use b^ = 6a^ for 
the bond length of the Gaussian reference chain of the perturbation 
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q-k 
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Fig. 48 Interaction diagrams used in reciprocal space for tiie calcula- 
tion of 8G{q,t) in the scale free limit. There exist three nonzero contri- 
butions to first-order perturbation, t he fir st involving two points inside 
the segment (first two lines of Eq. {T82)), the secon d one point inside 
and one outside the segment (third line of Eq. {T82))and the third one 
point on either side of the segment (last line of Eq. ( I182t ). Momentum 
q flows from one correlated point to the other. Integrals are performed 
over the momentum k. Dotted lines denote the effective interactions 
v(^) given by Eq. <66t . bold lines the propagators which carry each a 
momentum q or q — k as, indicated. 



calculation and b~ for the measured effective bond length, Eq. \\1\ . To 
compute the different integrals it is more convenient to work in Fourier- 
Laplace space {q,t) with t being the Laplace variable conjugate to s: 



I 



8G{q,t)= / ds8G(q,s)e- 



(181) 



As illustrated in Fig.|48l there are three contributions to this perturba- 
tion: one due to interactions between two monomers inside the segment 
(left panel), one due to interactions between an internal monomer and 
an external one (middle panel) and one due to interactions between two 
external monomers located on opposite sides (right panel). In analogy 
to the derivation of the form factor described in Ref. 1831 this yields: 
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where we have set Q = {aoq)^ +t. The graph given in the left panel 
of Fig. |48] corresponds to the first two lines, the middle panel to the 
third line and the right panel to the last one. Seeking for the mo- 
ments we expand 8G{q,t) around q = 0. Having in mind chain strands 
counting many monomers (s 2> 1), we need only to retain the most 
singular terms for r ^^ 0. Defining the two dimensionless constants 
c=i37i^f^alp)-^ = ^24/7t^/blp ande = v<^/3m;* = Uw^/Tth'}^ this 
expansion can be written as 



8G{q,t) = 



1 -r(2) 

1! f2 

2 r(3) 

2! f3 

3r(4) 

3! f4 



e (cioq) 
e (ao'?)'* 



{aaqf 



1 ^(3/2) , ,2 
c (aoq) - 



1! f3/2 

1 16 r(5/2) 

2!!" f5/2 

1 216r(7/2) 
3! Is" f7/2 



(183) 



c {aoqf 



{amf 



where we have used Euler's Gamma function F(a) |441. The first lead- 
ing term at each order in q^ — being proportional to the coefficient e — 
ensures the renormalization of the effective bond length feo — > b. The 
next term scaling with the coefficient c corresponds to the leading fi- 
nite strand size correction. Performing the inverse Laplace transforma- 
tion r{a)/t" — ^ 5"^' and adding the Gaussian reference distribution 
Go{q,s) this yields the A2p-coefficients for the expansion of G{q,s) 



around q = 0: 
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More generally, one finds 
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:■! """"' \' ' "" 2(2p+l)! ^ 
From this result and using Eq. {178) one immediately verifies that 

^- ^ 6Pp\ ^ ' \ 2{2p+\)\ ^fs \b J ! 



where we have defined 



12 v<^ 



= bk\\ 



12 



Gzo 



(187) 



for the effective bond length with Gzo = ^/vp/b^p. We remind that bg 
refers here to the bond length of the Gaussian reference chain while b 
is the effective bond length for asymptotically long chains. Since our 
key interest is not as for Edwards to predict b t4;43'J but to describe 
the deviations for finite s from the Gaussian limit fors^-N which 
sets the relevant reference length scale. Setting thus bp = b this leads 
to Eq. <76t with b being a fit parameter Using Eq. l |179t one justifies 
similarly Eq. <78t for the Gaussianity parameter Up. 

The moments Eq. {T86) completely determine the segmental dis- 
tribution G{r, s) which is indicated in Eq. l l73t . While at least in prin- 
ciple this may be done directly by inverse Fourier-Laplace transforma- 
tion of the correction 8G{q,t) to the generating function it is helpful to 
simplify further Eq. il82t . We observe first that 8G{q,t) does diverge 
for strictly incompressible systems (v -^ 00) and one must keep v finite 
in the effective potential whenever necessary to ensure convergence 
(actually everywhere but in the diagram corresponding to the interac- 
tion between two external monomers). Since we are not interested in 
the wave vectors larger than l/E, we expand 8G{q,t) f or 1^ — > which 
leads to the much simpler expression 
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with Q = (aoqY +t. The first term diverges as -/v for diverging v. It 
renormalizes the effective bond length in the zero order term which is 
indicated in the first line of Eq. <73t . The next two terms scale both 
as v". Subsequent terms must all vanish for diverging v and can be 
discarded. It is then easy to perform an inverse Fourier-Laplace trans- 
formation of the two relevant v^-terms. This yields 



^ni \ n I \ ^ 3x/Jr / 2 3x x^ 
5G(.,.) = Go(x,..)-^^^-- + --- 



(189) 



with X = rja^'J's = \/6«. This is consistent with the expression given 
in Eq. ((73). 
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B.3 Intramolecular form factor F{q) 

The deviation 5F {q) = F (q) — Fo{q) oi the intramolecular form factor 
from the Gaussian reference Fo{q) can be readily obtained from the 
deviation 5G{q,t) of the subchain size distribution in Fourier-Laplace 
space. Since G{q,s) = <^exp(— ig-r)) we have for asymptotically long 
chains 

SF{q) = 2 fds 5G{q,s) = 2 5G{q,t = 0) (190) 

where we focus on the intermediate scale-free wavevector regime which 
does not depend on the c hain length distribution pN. Using Eq. (T| 
derived in Appendix lB.2l one obtains 



5F{q) : 






9 1 1 
4/?^ bq' 



(191) 



where only the third term of Eq. J188t contributesQNote that b = \/6a 
stands now for both the measured effective bond length and the b ond 
length of the Gaussian reference chain. It follows from Eq. ( |191t that 
within first-order perturbation theory 



F{q)=Fo{q) + 5F{q)^Fo{q){l 



3 bq 
Sb^ 



(192) 



which is equivalent to the prediction Eq. (141 made in the Introduction. 
Note that this perturbation result is consistent with the renormalization 
group calculations of semidilute solutions by Schafer |7 71]. This is 
of course expected since semidilute solutions may be considered as 
incompressible melts of (semidilute) blobs [3}. 



B.4 Angular correlations Pi (s) 

We present here the direct perturbation calculation of the bond-bond 
correlation function Pi (s) of Flory-distrib uted lin ear chains in recipro- 
cal space following the discussion in Sec. 12.4.61 As sketched in panel 
(b) of Fig. [5] one only needs to compute the interaction diagram be- 
tween the monomers of the two dangling tails if the two bonds /j and 
I2 are placed outside the subchain of length s connecting the head of 
the first bond / j with the tail of the second bond l^ . To simplify the 
notations we set immediately bo = b and ao = a and do not distinguish 
between the bond / of the computer model connecting the monomers 
and the effective Gaussian bond length b, i.e. we set c„ = (b/l)- = 1. 
Restating Eq. dSOt the relevant interaction diagram /o(.s) readj^ 



PlW=/oW = (-l 
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dq 

j2Ky 



'Go{q,s).s/{q,q)v{q)w{q). 



(193) 



The negative sign in front of the integral is due to the negative sign of 
the general perturbation calculation formula Eq. l l57t . The first factor 
Go{q,s) = exp{—{aq)^s) in the integral over the wavevector q stands 
for the Fourier transformed Gaussian propagator Go(r,s) between the 
two bonds indicated in Fig.[5jb). Using Eq. l [82t the scalar product of 
the bond vectors in real space is represented by the scalar product of 
the wavevectors flowing through the bonds. Note that the momentum 
flows in the same direction as the two bonds and the operator thus reads 
s^[q,q) = —(bq/d)'^. We assume that chains are Flory-distributed and 

'" The first term in Eq. {TSS} must be discarded as before since it 
only renormalizes the effective bond length and would also contribute 
to the reference form factor Fo{q). 

'' The diagram loiq) = {—l)Go{q,s)i^{q,q)\{q)w{q) in reciprocal 
space is very similar to the diagram (/) indicated in Fig.|49]for the cal- 
culation of the bond-bond correlation function Pi (r) as a function of 
distance r. Since for Pi {s) we want to sample at constant s between the 
bond pairs irrespective of the distance between the bonds no momen- 
tum needs to be injected in the diagram and the momentum flowing 
along the diagram is constant everywhere. 



that the effective interaction y{q) between two monomers in the two 
dangling ends is described by Eq. l l65l l. The combinatorics between the 
interacting monomers — corresponding to the second line in Eq. l |79l l 
— leads using Eq. <29t to an additional weight factor 



w{q) = 
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(194) 



for the two Flory-distributed dangling tails. The task is thus to compute 
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where we have used that i,- = o^g/2 = a~ /2\p, Eq. d44t . Rewriting 
the last factor in Eq. (195) as 
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one may rewrite the interaction integral as 
I 1 1 



Pi{s) = - 



[A{x = iis)-A{x = 2s/g)] 



dp a'' s^/^ 1 - fig/ 2 ' 
We remind that the integral A(x) takes the asymptotics 
1 
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and hm A(x) ~ lim 1/x = 0. 



(199) 



limA(x) = ^. ^,.- 

Considering thus the limit of infinite Flory-distributed poly mers ( x = 
lis — )• 0) and incompressible solutions (x = 2s /g — > 00), Eq. {198} re- 
duces to 
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with ft) = d/2 in agreement with Eq. jTT}. Evaluating the integral A{x) 
over the wavevector ind = 3 dimensions one obtains for general x 



Ax) = 
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with erfc(x) being the complementary eiTor function |44l. Setting cp = 
\/2A/7?JS{b^p) in agreement with the power-law amplitude indicated 
in Eq. J9| and using u = s/g for the reduced subchain length it follows 
from Eq. i fTW} that 



P M = 



cp 



?3/2 1-^^/2 



^(l-M«/2) 



4\/2^e2"erfc(\/2l^) -I- 2(/xg)^/2\/^e'''erfc( ^TIJ)] . (202) 



Since implicit to the effective interaction potential Eq. J65) we have 
^ig/2 <g; 1 and since we focus on short subchains with /x.s = s/ (N) ^ 1, 
this further simplifies to 



PYs) = 



cp 
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(203) 



For large chains, ;U — > 0, the last term in the bracket decays rapidly 

as 1/ (N) ' and can be omitted for reasonable (mean) chain lengths. 
This leads to Eq. J83) which we have checked numerically in Fig. [26] 
for monodisperse chains. The first term in the bracket dominates only 
for small u when the structure within a (very) large thermal blob is 
probed. The bracket reduces to A/y/u- {4/y/u- l/w'l^ + ...)» \/u^l'^ 
in the opposite large-M limit as may be see n by expanding the error 
function |44|. Eq. l i203t thus reduces to Eq. l l200t for rf = 3. This thus 
confirms the key relation Eq. l|9} for the power-law decay of the bond- 
bond correlation function for all g on scales larger than the thermal 
blob (s/g-> I). 
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rings (r) 



linear chains (I) finite chain length (*) 



Fig. 49 Sketch of interaction diagrams in reciprocal space used for 
the computation of the bond-bond con-elation function P[ (r) summing 
over all intrachain bond pairs at constant r = |£|. The injected wavevec- 
tor q is conjugated to the distance r. The bold vertical aiTows represent 
the bond vectors in reciprocal space, Eq. <8U . dashed lines the effec- 
tive monomer interactions y{q), and the thin lines the Fourier-Laplace 
transformed Gaussian propagators G(q,t) with t being conjugated to 
the curvilinear distance s. Cyclic rings are described by the diagram 
(r), the behavior of asymptotically long linear chains by the diagrams 
(/). The last diagram (*) describes the finite-size corrections for Pi (r) 
relevant for large distances r 2> r,. 



chains stems from the minus sign implied by the fi rst o rder perturba- 
tion, Eq. <57t . Using Eq. (82) and assuming Eq. (205) the integrals 
simplify considerably 



h{q) = 
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Using the well-known theorem for the (inverse) Fourier transforma- 
tion of convolutions and that ■^{drf{r)] = iqf{q), the inverse Fourier 
transforms are thus 



/iW = 



1 



2p 

jb/df 
2p 



■ lAr] 



{d,.G(r)f 



(210) 
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with G{r) being the probability distribution to find another monomer 
of the chain around a reference monomer as given by Eq. (30). Inter- 
estingly, up to a constant prefactor the integrals /, and I\ are thus equal 
on large scales (r 2> (^). 



B.5 Distance dependence of angular correlations 
B.5.1 Interaction diagrams 

We focus now on the low-wavevector regime where the melt can be 
considered as incompressible (at least in terms of thermal blobs) and 
assume a Flory-distributed bath of a given inverse mean chain length 
li. The corresponding results for infinite chains may be obtained by 
setting /i = 0. The bond-bond correlation function Pi (r) = (/j -lij /fi 
is obtained by averaging over all intrachain pairs of bonds /j and I2 
at a given distance r = \r\ irrespective of their curvilinear distance s. 
Setting thus f = for the Laplace variable conjugated to s, the summed 
up Fourier-Laplace propagator becomes according to Eq. (29) 



GW.Go(M)U = py 



(204) 



Using Eq. (60) and Eq. (63) the effective interaction potential thus 
reads (87J 



vWGW = 
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on large scales k^l/a and k<Sil/^. The momentum q inserted in the 
interaction diagrams shown in Fig. [49] is conjugated to the distance r 
between both bonds. Momentum is a conserved quantity flowing from 
one correlated point to the other. The wavevector associated with the 
first bond at the origin is called q , the wavevector through the second 
bond q . As may be better seen from the sketch included in Fig. [27] 
the first bond is now within the .s-subchain while the second bond is 
chosen to be outside. Putting both bonds within the subchain leads to 
trivial correlations even for Gaussian chains, putting both bonds out- 
side the subchain leads to a vanishing correlation f or (^ — > as already 
remarked in Footnote [38] for the correlation function Ci{r) |89|. The 
first two diagrams in Fig.|49]for closed rings (r) and linear chains (/) 
are given by the convolution integrals 
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We remind that for closed cycles the bond-bond correlation function 
does not vanish even for perfect Gaussian statistics and it is the corre- 
sponding zero order average (^Z) Ri (•^)o which is computed by the 
first integral. The minus sign in front of the second integral for linear 



B.5. 2 Normalization and bond-bond correlations 

Focusing from now on the 3D case the bond-bond correlation function 
of closed rings Pr(r) is obtained from /r(r) after normalization with 
G{r). Using Eq. (sb) this yields 
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with X = y/Jir/2a comparing the distance r to the typical (z-averaged) 
size i?g_z w a^ {N) of Flory-distributed chains. The reason for the nor- 
malization factor G^(r) is that for P^r) both bonds are known to be 
bonds of the same polymer ring while the interaction integral Eq. (206) 
corresponds only to a probability G{r) for both bonds being in the same 
chain times a probability G{r) that this chain is closed. That Pr(r) is 
negative is of course due to the closure constraint which corresponds 
to an entropic spring force bending the second bond back to the origin. 
For large chain lengths the factor (1 + Ixf- can be neglected and the 
bond-bond correlation function Pi(r) becomes a scale- free p ower law. 
For linear chains it follows immediately from Eq. (2TT) that 

being the limit for asymptotically long chains and using the scaling 
function h{x) = (1 -|-2x)^exp(— 2x) for the finite-/! corrections. The 
normalization factor G{r) is due to the fact that for Pi(r) both bonds 
are known to belong to the same chain while the interaction integral 
Eq. (207) corresponds only to a probability G(r) for both bonds to 
be on the same chain. As compared to the closed cycles the correla- 
tion has the opposite sign since the attractive spring of the ring closure 
indicated by G{q ) in Eq. (208) has been replaced by the effective re- 
pulsion indicated by -G{q^)\{q^)G{q^) = -G{q^)/2p in Eq. (209) . 
This repulsion bends the second bond away from the origin increasing 
thus the bond-bond correlation function. 



B.5. 3 Sum rule and geometrical interpretation 

Interestingly, the perturbation result, Eq. (210), for Flory-distributed 
chains may be rewritten as 



^lW + ^^rW=0 



(214) 



57 



where we have used the normahzation factors mentioned above. This 
"sum rule" suggests a geometrical interpretation of the observed re- 
lation between infinite linear chains and closed cycles which may re- 
main valid beyond the one-loop approximation used here|3 The idea 
is that in an hypothetical ideal melt containing both linear chains and 
closed cycles all correlations disappear (on distances much smaller 
than the typical chain sizes) when summed up over the contributions 
of both architectures. The weight (G(r)/2)/p corresponds to the frac- 
tion of bond pairs in closed loops|3 Since the orientational coitc- 
lations in ideal cycles are necessarily long-ranged due the ring clo- 
sure (Eq. I212t . it follows, assuming the sum rule, that the same ap- 
plies to bond pairs of linear chains. Since bonds in closed cycles are 
anti-correlated {Pr('') < 0), they must be aligned (Pi (r) > 0) for linear 
chains. Interestingly, if one turns the argument around ass umin g the 
sum rule (rather than deriving it as we did) this imposes Eq. <205t and, 
hence, the effective intrachain potential \/{k) = {{akf- + ii)/2p for a 
Flory-distributed linear chain bath. 



B.5.4 Finite chain size ejfects for Flory-distributed chains 

We emphasize that the diagram (/) shown in Fig.|49]is not sufficient to 
characterize Pi {r) for larger distances since the last diagram (*) corre- 
sponding to the convolution integral 



r rllr 
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provides, as we shall see, the actu al cu toff of the power law in this 
limit. Using again Eq. {82} and Eq. ( l2Q5t the integral factorizes 
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where we have introduced in the last line the convenient dimensionless 
constant 



c* : 



{b/a) 



18p6 






(218) 



in which we dump local physics at large wavevector k. Before eval- 
uating the angular correlations in real space it is important to clarify 
the phy sics described by the diagram. The underlined second factor 
in Eq. J2T6) characterizes the alignment of the bond vectors of the 
monomers nj and nj — 1 at a fixed distance r of the monomers «| and 
n2 = ni+sa& shown by the sketch included in Fig.|27] Obviously, even 
for Gaussian chains these two bonds become more and more aligned if 
the distance r = \r\ gets larger than fo'/^, i.e. when the chain segment 
becomes stretched. For perfectly Gaussian chains the bonds /[ and [-, 
aln\ and 112 would still remain uncorrelated, however, since the second 

^^ A similar sum rule is obtained for the bond-bond coixelation 
function Pi (s) as a function of arc-length s. One verifies readily that 
Ii{s) = -Ir(s)/2p and thus 
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being the density of chains of 



with P,{s) 

length A' > .s in a Flory-distributed bath which would return after A' 

steps to the origin. This is the most elegant way to demonstrate that 

Pl(.s)~l/.v''^ 

'^ Note that G{r) is the density of the monomers of both strands the 
reference monomer is connected to. We know for linear chains as for 
cycles that both bonds are connected by a first strand. The probability 
for both bonds to be in a closed loop is given by the density G{r)/2 of 
the second strand. The factor 1/2 is thus needed to avoid counting the 
same ring twice. 



bond is outside the chain segment on which we have imposed the dis- 
tance constraint. As indicated by the dashed line in the diagram, it is 
then due to the effective interaction between the monomers within the 
stretched segment (n < nj) and the monomers outside (n > ni) that the 
bonds at ^2 — 1 and ni get aligned and then in turn the two bonds at ni 
and «2- We note that, strictly speaking, c* depends on the mean chain 
length (A'), since G(k) is a function of ;U. However, one checks readily 
that this effect can be neglected for reasonable mean chain lengths. We 
also note that the constant c^ infinite, since the UV divergence which 
formally arises for large k (where c* ~ k) may be regularized by local 
and, hence, model dependent phvsicsP^ We determine C j, num erically 
from our simulations of self-assembled linear EP in Sec. 15.4.21 

Assuming a finite and chain length in depen dent coefficient c* in 
Eq. J217t and inserting the propagator Eq. l l204t we obtain by inverse 
Fourier transformation 



I,{r)=c,{^G{r)-5{r)) 



(219) 



for the interaction integral in real space. Normalizing /^^(r) as before 
with G(r) and summing over both diagrams for linear chains this yields 



Pl(r) =P^{r)h{x)+c^,^l 



(220) 



for r 2> 1^ > 0. Comparing both terms in Eq. i220t one verifies that a 
crossover occurs at r, Ri fe (A') ' in agreement with Eq. l ll24t stated in 
the main text. The bond-bond correlation function P| (r) of an incom- 
pressible solution of Flory distributed polymers becomes thus constant 
for r » r,. This remarkable result is essentially due to the polydis- 
persity. This allows to find for all distances r pairs of bonds /i and [2 
stemming from segments which are slightly stretched by an energy of 
order /x <g; 1 and which are, hence, slightly shorter than a unstretched 
segment of length s ^ [r/b)^. Since there are more shorter chains and 
chain segments this just compensates the decay of the weight due to 
the weak stretching. Although the number of such slightly stretched 
segments decays strongly with distance, their relative effect with re- 
spect to the typical unstretched segments, e** — 1 » /x, remains con- 
stant for all r. It is for this r eason that the chemical potential appears 
in the second term of Eq. <220t . Please note that bond pairs from 
strongly stretched segments (corxesponding to an energy much larger 
than /x) are, however, still exponentially suppressed and can be ne- 
glected. As explained in detail in Ref. |87 1, this is qualitatively differ- 
ent for monodisperse chains where strongly stretched chain segments 
contribute increasingly to the average for large distances. 



B.6 Non-extensivity of the chemical potential 

We take again as a reference for the perturbation calculation a melt of 
Gaussian chains where the bond length b is set by the effective bond 
length of asymptotically long chains. Averages performed over this un- 
perturbed reference system are labeled by an index 0. The task is now 
to compute the ratio Q{n)/Qo{ii) of the perturbed to the unperturbed 
partition function of the test chain of length n plugged into the bath of 
A'-chains, 



1 



Q{n 



Qo{n) ^ ^° 

" f 

("«)o = £(«-■«) / drGo{r,s)v{r) 



(221) 



with the perturbation potential u„ being the sum of the effective mono- 
mer interactions v(r) of all pairs of monomers of the test chain. The 

''* For soft melts with weak bare excluded volume v = 1/gp, i.e. £, 2> 
b, it can be shown that an upper cutoff wavevector kc~ \/t, regularizes 
the integral over k. The coefficient becomes c,^ = \/2Kpb^t,. Using 
b = 3.244 and t, = 0.5 for p = 0.5/8 this gives c^ » 0.48. This is 
not that far off the best fit value c* » 0.14 considering that for small 
i, <b one expects c,^ to be rather determined by the model-depending 
stiffness between adjacent bonds. 
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factor n — s in Eq. <221t counts the number of equivalent monomer 
pairs separated by an arc-length s. The deviation 5^„ from Flory's hy- 
pothesis is then given by the cont ribution to (m„)q which is non-linear 
in n. The calculation of Eq. J22U in d dimensions is most readily per- 
formed in Fourier-Laplace space with q being the wavevector conju- 
gated to the monomer distance r and l the Laplace variable conjugated 
to the chain length n. The Laplace transformed averaged perturbation 
potential reads 



References 
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Jn=0 



dn{u„)f^e 



I 



(2^ 



Go{q,t)w{q)w{t) 



(222) 



6. 



where Go{q,l) represents the Fourier-Laplace transformed Gaussian 
propagator Go{r,s) as given in Eq. l [28t and the weight factor w{t) = 
l/t^ accounts ag ain fo r the combinatorics between the interaction monomers -j 
— just as in Eq. l[T93). 

The effective interaction potential \{q) being given by E£j59) g 

one realizes that the naive perturbation calculation using Eq. l l222t is 
formally diverging at high wavevectors in three dimensions (becoming g 

only regular below d = 2). This is due to monomer self-interactions 
which must be first subtracted. Using Eq. <64t instead of Eq. l |59t even 
makes things worse due to an additional divergency associated with 
the self-interactions of the blobs whose size was set to zero (g -^ 0). 
However, since we are not interested in (possibly diverging) contribu- 
tions linear in the length of the test chain or independent of it, we can 
freely subtract linear terms (i.e., terms ~ \/t^ in Laplace space) or con- 
stant terms (i.e., terms ~ \/t) to regularize and to simplify ;(,. Such a 
transformation leads to 



Ut = 
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d'^q 1 2F^^ - {aqf 



{iKY'fi [aq)^+t 2(vp-hF(,-') 



(223) 



where ". . ." stands for the linear and constant contributions we do not 
compute. This conver ges n ow for incompressible melts in d < 2 dimen- 
sions. Applying Eq. <223t to incompressible Flory-distributed melts, 
i.e. assuming Eq. <40K this yields 



I ll-t r d''q 
2p~P~ J (InY 



2p 



G{q,t) 



(2j 

l/t)G{r = 0,t)- 



(224) 



2p 

where we have read Eq. l l224t as an inverse Fourier transform taken 
at r = 0. Remembering that a factor 1/f in f-space stands for an inte- 
gral /q ds in ;7-space, the inverse Laplace transform of u, can be ex- 
pressed in terms of integrals of the return probability G(r = 0,s) = 
{Ansa^Y''l'^. We obtain, hence. 



, in n-space 



5Ain = 
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where 8^^ stands for the non-extensive contribution to {u„)q. Note that 
the first term in the brackets scales as the correlation hole in d dimen- 
sion. Its marginal dimension is d = 2. The second term characterizes 
the effective two-body interaction of the test chain with it self. A s one 
expects 1 3 1, its marginal dimension is rf = 4. Although Eq. J226t is for- 
mally obtained for d < 2 it applies to high er dimensions by analytic 
continuation. In three dimensions Eq. ([226) becomes 



5tla = 



1 



1 



(4;r)3/2 pfl3 



-1/2. 



2/x/j 



1/2 



(227) 



which demonstrates finally the non-extensive correc tion to the ideal 
polymer chain chemical potential announced in Eq. l llOSt and repre- 
sented by the solid lines in Fig. fTsl The chemical potential for monodis- 
perse chains is obtained from Eg. ilOSt within the Fade approximation, 
Eq. l [37t . where fi is replaced by 2/N. This result is indicated by the 
dash-dotted lines in Fig. [Ts] A more precise calculation for monodis- 
perse chains using the full Debye function is given in Ref. 1881 . 
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